search

esrivas17 / gacos_correction_for_GMTSAR

Atmospheric GACOS correction adapted for GMTSAR
17 stars 6 forks source link

GACOS ERROR #5

Open CUMT530849827 opened 1 year ago

CUMT530849827 commented 1 year ago

Hello, author, I would like to ask about the error I reported during the operation of GACOS correction. The following is the specific log of the error error1: rdtrack [WARNING]: Some input points were outside the grid domain(s). grdtrack [WARNING]: Some input points were outside the grid domain(s). grdmath [ERROR]: Operation "SUB" requires 2 operands grdmath [ERROR]: Operation "SUB" requires 2 operands

error2: surface [INFORMATION]: Writing grid to file lla.grd grdtrack [WARNING]: Some input points were outside the grid domain(s). grdtrack [WARNING]: Some input points were outside the grid domain(s). grdmath [ERROR]: Operation "SUB" requires 2 operands grdmath [ERROR]: Operation "SUB" requires 2 operands Problems with the reference point. Is the reference point within the AOI? GACOS correction done

image

image image ![Uploading image.png…]()

CUMT530849827 commented 1 year ago

ref_point 101.9389 37.4718

esrivas17 commented 1 year ago

Are your GACOS grids slightly larger in longitude with respect to the DEM? Maybe trying with a larger DEM?

Elcappo2019 commented 10 months ago

@esrivas17 , I faced the same problem, the limits of dem is more than the GACOS grid file !!! but still the same error , Can you help me to resolve it ? error1 demVsdate1_ztd

esrivas17 commented 10 months ago

Hiii. According to the log, it seems the problem comes from the reference point. Could you check it is in the right format and that it is within the AOI?

Elcappo2019 commented 10 months ago

ref_point

esrivas17 commented 10 months ago

Ok. It seems that you are processing 3 subswaths. Originally, the scripts were not meant to work with more than 2 subswaths. I would need to add this feature to the code. I think that it is crashing when reprojecting GACOS to radar coordinates.

esrivas17 commented 10 months ago

Also, there are some inaccuracies when reprojecting the grid from geographic coordinates to radar. Meaning that when the limits are too close to each other (between Gacos and Dem) when reprojecting some points might fall outside of the grid. You could try with a new Dem that is slightly longer in the Y minimum boundary

esrivas17 commented 10 months ago

Also, there are some inaccuracies when reprojecting the grid from geographic coordinates to radar. Meaning that when the limits are too close to each other (between Gacos and Dem) when reprojecting some points might fall outside of the grid. You could try with a new Dem that is slightly longer in the Y minimum boundary

You can better try first to run line 73 from operation.csh and to see if it runs successfully or not

Elcappo2019 commented 10 months ago

Dear @erik-17 , I run the line 73 in the operation.csh in the top terminal. and i think is Ok (see the bottom terminal" . run173

esrivas17 commented 10 months ago

Dear @erik-17 , I run the line 73 in the operation.csh in the top terminal. and i think is Ok (see the bottom terminal" . run173

Right. Looks fine. One other thing. Could you check the reference point is within phasefilt.grd ? It looks like the Gacos grid dimensions are fine.

Also, what script are you using? Single_gacos_correction.csh?

Elcappo2019 commented 10 months ago

I confused, the limits of phasefilt_ll.grd appeared as in the first attached terminal goog i think the problem is the study area is located in the west of (Greenwich line) with " minus coordinates " and i put the reference point as " -8.077282 31.504392 " but the limit of " phasefilt_ll.grd are " x_min: 350.638888889 x_max: 353.588888889 " and " y_min: 30.2958333338 y_max: 32.3250000005" So, How to fix it ?

Elcappo2019 commented 10 months ago

I used " operation.csh " script .

Elcappo2019 commented 10 months ago

When i make a new phasefilt.grd file and make sure the coordinates with X in minus signs fitting the coordinate of ref_point. in addition, during this case i used "Single_gacos_correction.csh " tt Still this error !!!! where is the problem ?

esrivas17 commented 10 months ago

When i make a new phasefilt.grd file and make sure the coordinates with X in minus signs fitting the coordinate of ref_point. in addition, during this case i used "Single_gacos_correction.csh " tt Still this error !!!! where is the problem ?

Mmm everything looks fine. Could you check again the info of resample_zpddm.grd? And ensure it has the same dimensions and incs of phasefilt_ll.grd? If everything is fine you can check what output you get by running gmt grdtrack using your reference point to phasefilt_ll.grd and resample_zpddm.grd It's indeed strange. I can't see an obvious mistake. But I have spotted some misleading information in the script