<div dir="ltr"><br clear="all"><div>Hi all,<br><br>I'm trying to figure out how to convert GOES 16 & 17 weather data into a form readable by GeoTools.  It looks like I'm possibly loosing the sweep information.  Take for example this set of 16 bands:<br><br>gsutil cp gs://gcp-public-data-goes-17/ABI-L1b-RadC/2019/151/21/OR_ABI-L1b-RadC-*_G17_s20191512131197_*.nc .<br><br>if I do:<br><br>gdalwarp --version<br>GDAL 2.2.4, released 2018/03/19<br># e.g. for all C1..C16<br>gdalwarp -t_srs EPSG:4326 NETCDF:OR_ABI-L1b-RadC-M6C16_G17_s20191512131197_e20191512133581_c20191512134014.nc:Rad OR_ABI-L1b-RadC-M6C16_G17_s20191512131197_e20191512133581_c20191512134014_4326.tif<br><br>The imagery once tiled by and displayed by Earth Engine, it shows up in the correct location.<br><br>But I would rather have the data stay in it's native projection and only be transformed when used to which ever projection is needed.  (this is typical of Earth Engine)<br><br>gdal_translate NETCDF:OR_ABI-L1b-RadC-M6C16_G17_s20191512131197_e20191512133581_c20191512134014.nc:Rad OR_ABI-L1b-RadC-M6C16_G17_s20191512131197_e20191512133581_c20191512134014.tif<br><br>Gives the imagery not quite in the correct location.   I briefly looked at gdal_translate from proj at head and gdal at head, but our tiler is an older version of gdal that can read the post "RFC 73: Integration of PROJ6 for WKT2" info.<br><br>Any one have suggestions on how to build a geotiff using gdal 2.2.4 from the netcdf that keeps all the necessary projection info in a form that GeoTools 19 can understand?<br><br>I've attached a screenshot from the Gulf of California where an island is clearly offset between the renderings of the two assets (transformed from the gdal_translate geotiff and the gdalwarp to EPSG:4326 geotiff)<br><br>Thanks!<br>-kurt<br>Engineer on Google Earth Engine<br><br>Earth Engine screen shots:</div><div><a href="https://photos.app.goo.gl/sW1upk1iZ9p62MuY9">https://photos.app.goo.gl/sW1upk1iZ9p62MuY9</a></div><div><a href="https://photos.app.goo.gl/UAGccF2rjPDfenS27">https://photos.app.goo.gl/UAGccF2rjPDfenS27</a><br><br>diff -u <a href="http://netcdf.info">netcdf.info</a> <a href="http://tif.info">tif.info</a><br>--- <a href="http://netcdf.info">netcdf.info</a>     2019-06-26 10:54:22.512939636 -0700<br>+++ <a href="http://tif.info">tif.info</a>   2019-06-26 10:54:39.596731603 -0700<br>@@ -1,22 +1,23 @@<br>-Driver: netCDF/Network Common Data Format<br>-Files: OR_ABI-L1b-RadC-M6C02_G17_s20191512131197_e20191512133570_c20191512133595.nc<br>+Driver: GTiff/GeoTIFF<br>+Files: OR_ABI-L1b-RadC-M6C02_G17_s20191512131197_e20191512133570_c20191512133595.tif<br> Size is 10000, 6000<br> Coordinate System is:<br>-PROJCS["unnamed",<br>-    GEOGCS["unknown",<br>-        DATUM["unknown",<br>+PROJCS["Geostationary_Satellite",<br>+    GEOGCS["GCS_unknown",<br>+        DATUM["D_unknown",<br>             SPHEROID["Spheroid",6378137,298.2572221]],<br>         PRIMEM["Greenwich",0],<br>-        UNIT["degree",0.0174532925199433]],<br>+        UNIT["Degree",0.017453292519943295]],<br>     PROJECTION["Geostationary_Satellite"],<br>     PARAMETER["central_meridian",-137],<br>     PARAMETER["satellite_height",35786023],<br>     PARAMETER["false_easting",0],<br>     PARAMETER["false_northing",0],<br>-    EXTENSION["PROJ4","+proj=geos +lon_0=-137 +h=35786023 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs  +sweep=x"]]<br>+    UNIT["Meter",1]]<br> Origin = (-2505021.493779813405126,1583173.639055640902370)<br> Pixel Size = (501.004328871885207,501.004328871885264)<br> Metadata:<br>+  AREA_OR_POINT=Area<br>   goes_imager_projection#grid_mapping_name=geostationary<br>   goes_imager_projection#inverse_flattening=298.2572221<br>   goes_imager_projection#latitude_of_projection_origin=0<br>@@ -81,23 +82,16 @@<br>   y#scale_factor=-1.4e-05<br>   y#standard_name=projection_y_coordinate<br>   y#units=rad<br>-Geolocation:<br>-  LINE_OFFSET=0<br>-  LINE_STEP=1<br>-  PIXEL_OFFSET=0<br>-  PIXEL_STEP=1<br>-  SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]<br>-  X_BAND=1<br>-  X_DATASET=NETCDF:"OR_ABI-L1b-RadC-M6C02_G17_s20191512131197_e20191512133570_c20191512133595.nc":x<br>-  Y_BAND=1<br>-  Y_DATASET=NETCDF:"OR_ABI-L1b-RadC-M6C02_G17_s20191512131197_e20191512133570_c20191512133595.nc":y<br>+Image Structure Metadata:<br>+  COMPRESSION=DEFLATE<br>+  INTERLEAVE=BAND<br> Corner Coordinates:<br>-Upper Left  (-2505021.494, 1583173.639) (161d34'43.88"W, 14d47'44.39"N)<br>-Lower Left  (-2505021.494, 4589199.612) (175d33' 0.44"E, 53d31'38.26"N)<br>-Upper Right ( 2505021.795, 1583173.639) (112d25'16.11"W, 14d47'44.39"N)<br>-Lower Right ( 2505021.795, 4589199.612) ( 89d33' 0.41"W, 53d31'38.26"N)<br>+Upper Left  (-2505021.494, 1583173.639) (161d33'27.65"W, 14d49'58.11"N)<br>+Lower Left  (-2505021.494, 4589199.612) (175d46'57.22"E, 53d43' 3.62"N)<br>+Upper Right ( 2505021.795, 1583173.639) (112d26'32.33"W, 14d49'58.11"N)<br>+Lower Right ( 2505021.795, 4589199.612) ( 89d46'57.19"W, 53d43' 3.62"N)<br> Center      (       0.151, 3086186.626) (136d59'59.99"W, 29d58' 0.14"N)<br>-Band 1 Block=10000x1 Type=Int16, ColorInterp=Undefined<br>+Band 1 Block=10000x1 Type=Int16, ColorInterp=Gray<br>   NoData Value=4095<br>   Unit Type: W m-2 sr-1 um-1<br>   Offset: -20.2899112701416,   Scale:0.158592373132706<br></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"></div></div>