[gdal-dev] gdal-dev Digest, Vol 124, Issue 17

Thu Sep 11 19:54:59 PDT 2014

Hi,

The gdalinfo in my post is the gdalinfo of the subdataset 37. I have also
tried gdal_translate but when I tried to use gdalwarp using the
gdal_translate output, it has error. I really need to use gdalwarp so that
I could project the image into epsg:4326(+prj=longlat). If I'm going to use
the gdal_translate I am ask to indicate the gcp which is not good in my
case because I have many hdf files in my local directory and maybe they
have different gcp and I need to batch process them later using the
gdalwarp with the +proj=longlat command. What could be the solution?

On Fri, Sep 12, 2014 at 3:00 AM, <gdal-dev-request at lists.osgeo.org> wrote:

> Send gdal-dev mailing list submissions to
>         gdal-dev at lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://lists.osgeo.org/mailman/listinfo/gdal-dev
> or, via email, send a message with subject or body 'help' to
>         gdal-dev-request at lists.osgeo.org
>
> You can reach the person managing the list at
>         gdal-dev-owner at lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gdal-dev digest..."
>
>
> Today's Topics:
>
>    1. Re: S57 and area reconstruction (Antti Castr?n)
>    2. Re: Extract subdataset from netCDF in a given directory
>       (Andre Joost)
>    3. Re: ERROR 1: Too many points (440 out of 441) failed to
>       transform, unable to compute output bounds. (Andre Joost)
>    4. Re: Extract subdataset from netCDF in a given directory
>       (Andre Joost)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 11 Sep 2014 14:10:42 +0300
> From: Antti Castr?n <antti.castren at iki.fi>
> To: vwf <vwf at vulkor.net>
> Cc: gdal dev <gdal-dev at lists.osgeo.org>
> Subject: Re: [gdal-dev] S57 and area reconstruction
> Message-ID:
>         <
> CAPcvoO+scqr+CCO8JbYss5SnAAGWLGGa22GX+KMk-xd0bPKVZQ at mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> Hello Frits,
>
> It is possible to (re)merge the objects across cell boundaries if they
> share common Feature Object IDentifier (FOID), which is a combination
> of AGEN, FIND and FINS subfields. Some ECDISes can hide the boundary
> depending on masking. Usage of FOIDs vary, so it might be possible for
> your datasets - if you are lucky.
>
> See Clause 3.1 of Edition 2.0 (November 2000) of the ENC Product
> Specification (S-57 Appendix B.1) and ENC Encoding Bulletin Number 33.
>
> Cheers,
>
>   Antti
>
>
> 2014-09-11 8:33 GMT+03:00 vwf <vwf at vulkor.net>:
> > Hello,
> >
> > My question is about S57, not gdal, but one of you may be able to help
> > me.
> >
> > When an area (like "prcare") is at the edge of a cell (map) there will
> > be a vertex on the edge of the cel. With several cells the full area can
> > have several lines going through the area. Is there a trick/technique to
> > reconstruct the outer edge or remove the extra lines?
> >
> > Thank you,
> > Frits
> >
> > _______________________________________________
> > gdal-dev mailing list
> > gdal-dev at lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
> --
> Antti Castr?n
> +358505997538
> antti.castren at iki.fi
>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 11 Sep 2014 16:39:11 +0200
> From: Andre Joost <andre+joost at nurfuerspam.de>
> To: gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] Extract subdataset from netCDF in a given
>         directory
> Message-ID: <luscmh$ps8$2 at ger.gmane.org>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Am 10.09.2014 20:02, schrieb VictoriaH:
> > Hello, I am using gdal_translate to extract a subdataset from a netCDF
> file
> > to a geotiff. If I run the command from the same directory where the
> netCDF
> > file is located, there is no problem:
> >
> > Filename: PM25_2001_02_average.nc
> >
> > Ex. gdal_translate -of GTiff HDF5:"PM25_2001_02_average.nc"://Data
> > /var/www/pm25/netCDF/output.tif'
> >
> > However, if I try to run gdal-translate and point to a netCDF file in a
> > different directory, I get an error:
> >
> > Ex. gdal_translate -of GTiff
> > /netCDF/files/HDF5:"PM25_2001_02_average.nc"://Data
> > /tiff/files/pm25/netCDF/output.tif'
> >
>
> I am not sure about your syntax at ://Data. This worked for me:
>
> gdal_translate -a_srs EPSG:4326 -a_ullr -180 -90 180 90 -of "Gtiff"
> HDF4_SDS:UNKOWN:"LISOTD_HRFC_V2.3.2013.hdf":0 test.tif
>
> with all files in the same directory. You have to run gdalinfo first to
> get the exact name of the subdataset including the number after the :
> The output filename needs to be separated by a blank.
>
> HTH
> Andr? Joost
>
>
>
>
> ------------------------------
>
> Message: 3
> Date: Thu, 11 Sep 2014 16:46:05 +0200
> From: Andre Joost <andre+joost at nurfuerspam.de>
> To: gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] ERROR 1: Too many points (440 out of 441)
>         failed to transform, unable to compute output bounds.
> Message-ID: <luscmh$ps8$1 at ger.gmane.org>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Am 11.09.2014 10:45, schrieb Love:
> > I know I'm not the first to post this topic and encountered this problem.
> > But I've search and read some posts related to this and tried their
> > solutions but no luck. I really can't figure out what is causing this
> > problem. How am I going to solve it? Badly need help!
> >
> > I have tried gdalwarp -wo SOURCE_EXTRA=100 -wo SAMPLE_GRID=YES and still
> no
> > luck.
> >
> > Here are the info:
> >
> > gdalwarp -t_srs "+proj=longlat +datum=WGS84"
> > HDF4_SDS:UNKNOWN:"A2014042042000.L2_LAC.SeAHABS.hdf":37 hdf.tif
>
> maybe source crs is missing.
>
>
> > *GDALINFO:*
> >
> > Driver: HDF4Image/HDF4 Dataset
> > Files: A2014042042000.L2_LAC.SeAHABS.hdf
> > Size is 563, 880
> > Coordinate System is `'
>
> That is for the whole dataset. Try gdalinfo on the subdataset 37 you want.
>
> I was successful with extracting the subdatset using gdal_translate, and
> work on with that:
>
> gdal_translate -a_srs EPSG:4326 -a_ullr -180 -90 180 90 -of "Gtiff"
> HDF4_SDS:UNKOWN:"LISOTD_HRFC_V2.3.2013.hdf":0 test.tif
>
>
> In your case, a_ullr might be different.
>
>
> HTH,
> Andr? Joost
>
>
>
> ------------------------------
>
> Message: 4
> Date: Thu, 11 Sep 2014 16:58:17 +0200
> From: Andre Joost <andre+joost at nurfuerspam.de>
> To: gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] Extract subdataset from netCDF in a given
>         directory
> Message-ID: <lusdbc$4t1$1 at ger.gmane.org>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Am 10.09.2014 20:02, schrieb VictoriaH:
> > Hello, I am using gdal_translate to extract a subdataset from a netCDF
> file
> > to a geotiff. If I run the command from the same directory where the
> netCDF
> > file is located, there is no problem:
> >
> > Filename: PM25_2001_02_average.nc
> >
> > Ex. gdal_translate -of GTiff HDF5:"PM25_2001_02_average.nc"://Data
> > /var/www/pm25/netCDF/output.tif'
> >
> > However, if I try to run gdal-translate and point to a netCDF file in a
> > different directory, I get an error:
> >
> > Ex. gdal_translate -of GTiff
> > /netCDF/files/HDF5:"PM25_2001_02_average.nc"://Data
> > /tiff/files/pm25/netCDF/output.tif'
> >
>
> According to http://www.gdal.org/frmt_hdf5.html, the syntax should be:
>
> HDF5:file_name:subdataset
>
> So it has little use to put the path before the HDF5:. I'm not sure if
> GDAL accepts the path before the file_name. Running the command in the
> same folder should be a safe thing. The output file might go to another
> folder.
>
> HTH,
> Andr? Joost
>
>
>
> ------------------------------
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
> End of gdal-dev Digest, Vol 124, Issue 17
> *****************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20140912/e9d06e33/attachment.html>