Hey Scott, You could extract that layer from the netCDF, burn it with gdal_rasterize as a geotiff, convert the geoTiff into a GMT grid, and then possibly add that back into the netCDF (somehow).
Sounds like a pain though. Roger -- On Mon, Sep 21, 2009 at 11:49 AM, Scott Lewis <[email protected]> wrote: > Ah, thanks. I must have missed that. Looks like I'll have to find > another way to accomplish what I'm trying to do then! > > Thanks for your help! > > Scott > > Hermann Peifer wrote: > > The driver doesn't support updating netCDF files. gdalinfo --formats > says: > > > > GTiff (rw+): GeoTIFF > > (...) > > netCDF (rw): Network Common Data Format > > > > The "+" indicates update support, which is obviously missing for netCDF > > format. > > > > Hermann > > > > > > -------- Original Message -------- > > Subject: Re:gdal_rasterize with a NetCDF file? > > From: Scott Lewis <[email protected]> > > To: Date: 21/09/2009 20:37 > > > >> Evan, > >> > >> Thanks, that seems to be getting me in the right direction: I'm no > >> longer getting the segfault error, or "ERROR 5" at all. > >> > >> However, now it's giving me a different error. The first line has the: > >> 0...10...20...30...40...50...60...70...80...90...100 - done. > >> message, like I would expect from when I do this with GeoTIFF files. > >> But following that I am getting a whole bunch of this error, repeated > >> many times: > >> > >> ERROR 6: WriteBlock() not supported for this dataset. > >> > >> > >> Doing a search online, I'm getting the impression that this might mean > >> that gdal_rasterize doesn't support writing to a NetCDF file? Is this > >> correct, or does the problem lie elsewhere with my calling of the > >> utility. > >> > >> Here's the command line I'm using now: > >> > >> /usr/local/FWTools-2.0.6/bin_safe/gdal_rasterize -i -burn 0 -l > >> "OGRGeoJSON" dummy.json > >> 'NETCDF:"Albedo.nc":Albedo_with_1400_Local_Time_of_Measurement' > >> > >> > >> Thanks again for the help! > >> > >> Scott > >> > >> > >> Even Rouault wrote: > >>> Selon Scott Lewis <[email protected]>: > >>> > >>> Scott, > >>> > >>> Anytime you want to access to a subdataset with a gdal utility, you > >>> must use the > >>> value of the relevant SUBDATASET_xxx_NAME metadata item, in that > >>> instance > >>> 'NETCDF:"Albedo.nc":Albedo_with_1400_Local_Time_of_Measurement' > >>> > >>> The segmentation fault however wasn't expected, so I commited fixes > >>> per ticket > >>> http://trac.osgeo.org/gdal/ticket/3146 to error out more nicely. > >>> > >>> Best regards > >>> > >>>> Hi, > >>>> > >>>> I'm still fairly new at using GDAL, but I am having trouble burning a > >>>> vector image onto a NetCDF file. I've had success with using > >>>> gdal_rasterize to burn a vector onto a GeoTIFF file, but when I try to > >>>> burn the same vector onto a NetCDF file, I get an error. > >>>> > >>>> Based on the error, I'm suspecting it's because it doesn't know how to > >>>> find the right "band" in the NetCDF File (I want the > >>>> "Albedo_with_1400_Local_Time_of_Measurement" subdataset), but although > >>>> I've looked at the documentation, I'm not sure how to specify this to > >>>> gdal_rasterize. > >>>> > >>>> Any help would be much appreciated! Below is the info. > >>>> > >>>> > >>>> > >>>> Here is the command I am trying to run: > >>>> > >>>> /usr/local/FWTools-2.0.6/bin_safe/gdal_rasterize -i -burn 0 -l > >>>> "OGRGeoJSON" dummy.json Albedo.nc > >>>> > >>>> > >>>> And the error I get: > >>>> > >>>> ERROR 5: GDALDataset::GetRasterBand(1) - Illegal band # > >>>> > >>>> /usr/local/FWTools-2.0.6/bin_safe/gdal_rasterize: line 9: 26950 > >>>> Segmentation fault $FWTOOLS_HOME/bin/`basename $TARGET` "$@" > >>>> > >>>> > >>>> > >>>> Here is the gdalinfo for Albedo.nc: > >>>> > >>>> Driver: netCDF/Network Common Data Format > >>>> Files: Albedo.nc > >>>> Size is 512, 512 > >>>> Coordinate System is `' > >>>> Metadata: > >>>> NC_GLOBAL#Conventions=CF-1.4 > >>>> NC_GLOBAL#institution=National Snow & Ice Data Center, Boulder, CO > >>>> NC_GLOBAL#title=AVHRR Polar Pathfinder Twice-Daily 5 km EASE-Grid > >>>> Composites Albedo > >>>> NC_GLOBAL#source=See the title and references for this information > >>>> NC_GLOBAL#comment=Not set at this time > >>>> NC_GLOBAL#references=Documentation available at: > >>>> http://nsidc.org/data/docs/daac/nsidc0066_avhrr_5km.gd.html > >>>> NC_GLOBAL#history=Mon, 21 Sep 2009 09:31:24: File created. > >>>> Subdatasets: > >>>> > >>>> > >>> > SUBDATASET_1_NAME=NETCDF:"Albedo.nc":Albedo_with_1400_Local_Time_of_Measurement > >>> > >>>> SUBDATASET_1_DESC=[1x244x242] > >>>> Albedo_with_1400_Local_Time_of_Measurement (8-bit integer) > >>>> SUBDATASET_2_NAME=NETCDF:"Albedo.nc":latitude > >>>> SUBDATASET_2_DESC=[244x242] latitude (64-bit floating-point) > >>>> SUBDATASET_3_NAME=NETCDF:"Albedo.nc":longitude > >>>> SUBDATASET_3_DESC=[244x242] longitude (64-bit floating-point) > >>>> Corner Coordinates: > >>>> Upper Left ( 0.0, 0.0) > >>>> Lower Left ( 0.0, 512.0) > >>>> Upper Right ( 512.0, 0.0) > >>>> Lower Right ( 512.0, 512.0) > >>>> Center ( 256.0, 256.0) > >>>> > >>>> > >>>> And for completeness, the dummy.json file (it's basically a triangle): > >>>> > >>>> { > >>>> "type": "FeatureCollection", > >>>> "features": [{ > >>>> "type": "Feature", > >>>> "properties": { "BASIN_ID": 1.000000, "AREA_SQKM": > >>>> 1739.016470 }, > >>>> "geometry": { > >>>> "type": "Polygon", > >>>> "coordinates": [ [ > >>>> [ -2000000, -2500000 ], > >>>> [ -2500000, -1400000 ], > >>>> [ -1500000, -1400000 ], > >>>> [ -2000000, -2500000 ] > >>>> ] ] > >>>> } > >>>> }] > >>>> } > >>>> > >>>> > >>>> > >>>> Thanks! > >>>> > >>>> Scott Lewis > >>>> National Snow & Ice Data Center > >>>> _______________________________________________ > >>>> gdal-dev mailing list > >>>> [email protected] > >>>> http://lists.osgeo.org/mailman/listinfo/gdal-dev > >>>> > >>> > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
