Frank,<br>
<br>Netcdf now support UTM and LCC. <br>
TM projection works fine with WCS data from NWGISS <span class="q">server </span>.<br>
<br>
I have checked for uniform grid spacing in projected values at three places in the image:<br>
left, middle and right <br>
If I establish that spacing is uniform accross the map, I enable the GeoTransform, <br>
otherwise it stays disable. <br>
<br>
But I need help with LCC metadata found in the test file RUC.nc. <br>
I used osr.SetLCC1SP with STD_PARALLEL as dfScale. (Probably wrong...) <br>
I just
don't know what to do with these values for Lamber Conformal Conic :<br>
<br>
standard_parallel=25<br>
longitude_of_central_meridian=-95<br>
latitude_of_projection_origin=25<br>
GRIB_earth_shape=spherical<br>
GRIB_earth_shape_code=0<br>
<br>
You can see that False_Easting an False_Northing are missing so, I set them to 0 which is probably wrong too.<br>
<br>
Finallly, for latitude/longitude data, I use WKT "WGS84". <br>
(I also check for uniform grid spacing as before ...)<br>
<br>
<br>
<br>
Regards,<br>
Denis<br>
<br><br><div><span class="gmail_quote">2005/8/26, Frank Warmerdam <<a href="mailto:fwarmerdam@gmail.com">fwarmerdam@gmail.com</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Denis,<br><br>The GMU WCS NetCDF file now reports:<br><br>Coordinate System is `'<br>Origin = (94.692708,6.738066)<br>Pixel Size = (0.00271231,-0.00268426)<br><br>Corner Coordinates:<br>Upper Left ( 94.6927079, 6.7380663
)<br>Lower Left ( 94.6927079, 4.8295543)<br>Upper Right ( 96.8978175, 6.7380663)<br>Lower Right ( 96.8978175, 4.8295543)<br>Center ( 95.7952627, 5.7838103)<br><br>But, I am not sure how this could be. The file is apparently
<br>projected in transverse mercator, but those are lat/long values.<br>The file presumably is not projected in lat/long so the lat/long<br>pixel size will not be consistent across the image. It is important<br>we only set the geotransform based on the lat/lon (or projected values)
<br>if we can establish that the pixel spacing is uniform across the<br>dimension maps.<br><br>In this file, we should be picking up the TM values from the xdim and<br>ydim arrays.<br><br>... phone call from you ...<br><br>
Well, I think we have covered this mostly on the phone now, but i'll<br>send what I had in the email so far.<br><br>Best regards,<br>--<br>---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam">http://pobox.com/~warmerdam
</a><br>and watch the world go round - Rush | Geospatial Programmer for Rent<br></blockquote></div><br>