X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 15 May 2003 04:16 AM by  anon
idl scripting
 0 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
15 May 2003 04:16 AM
    Hi all, I'm working on a batch script to produce warped sst images of a specific region from avhrr files, the calibration part works fine but I'm having problems with the warping part, I need to convert to a mercator projection but am unsure of the parameters that need to be entered (leading to my script failing) this is the part that is failing, where am I going wrong? Secondly is it possible to warp from input geometry, would that be a better way of going about the final warping and georeferancing, if it is how should I approach it. Thanks Neil Ackermann ;Warp the calibrated data to Sample Mercator ; ; Set the keywords ; Query file to get the size envi_file_query, fid, ns=ns, nl=nl, nb=nb ;set the dimensions to be warped dims = [-1l, 0, ns-1, 0, nl-1] ;set which band to be warped pos = [0] out_name = tmp3_name units = envi_proj_units_translate('Meters') ; proj = envi_proj_create(/geographic, unit=units) ; Parameters are a, b, lat0, lon0, x0, y0 Params = [6378160.0, 6356774.7, 0.000000, 18.000000, 500000., 10000000.] datum = 'cape' name = 'West Coast Mercator' proj = envi_proj_create(type=20, name=name, datum=datum, params=params) ;pixel size of 1100m by 1100m ps = [1100,1100] ; ; Call the doit ; envi_doit, 'envi_avhrr_warp_doit', fid=cl_fid, $ avhrr_fid=fid, dims=dims, pos=pos, $ out_name=out_name, method=8, zero_edge, $ background=0, grid=[51,500], proj=proj, $ pixel_size=ps ; ,r_fid=r_fid
    You are not authorized to post a reply.