Re: [gdal-dev] Problem with black edges to DOQQs using JPEG in Tiff compression

2017-02-12 Thread Brian Case
steve and all after reading my reply I noticed a major error "first, if your data is non lossy you need to nearblack -nb 0 -near 15" this should read "nearblack -nb 0 -near 0" Brian On Sun, 2017-02-12 at 13:16 -0800, Brian Case wrote: > Steve, > > first,

Re: [gdal-dev] Problem with black edges to DOQQs using JPEG in Tiff compression

2017-02-12 Thread Brian Case
Steve, first, if your data is non lossy you need to nearblack -nb 0 -near 15 however your DOQQ's dont have a collar so there is no reason to nearblack gdalwarp -dstalpha will crate an alpha band to mask out the collar after gdalwarp rotates the image. no nead to nearblack the output of gdal warp

Re: [gdal-dev] how to create just the msk file from a rgba vrt file

2014-02-06 Thread Brian Case
Duarte gdaltranslate -of outfrmt -b 1 -b 2 -b 3 -mask 4 infile outfile brian you can output to vrt's just dont try to mosaic masked or aphad files with vrts brian On Thu, 2014-02-06 at 13:20 +, Duarte Carreira wrote: > Hi there. > > > > Continuing with trying to cut processing times,

Re: [gdal-dev] Fast Pixel Access

2014-02-03 Thread Brian Case
-Jukka tileindex, mapserver, and the gdal wms driver On Mon, 2014-02-03 at 17:20 +, Jukka Rahkonen wrote: > Luke Roth gmail.com> writes: > > > > > Another thing that might speed up access is setting the config > option GDAL_DISABLE_READDIR_ON_OPEN = TRUE, either as an environment > varia

Re: [gdal-dev] Fast Pixel Access

2014-01-31 Thread Brian Case
evenr what about the use of a tileindex? seems an intersection with a set of polys first would be quick brian On Fri, 2014-01-31 at 19:30 +0100, Even Rouault wrote: > Le vendredi 31 janvier 2014 17:15:53, David Baker (Geoscience) a écrit : > > Dev's, > > > > I have a set of 55,501 bil file

Re: [gdal-dev] Can gdal_rasterize utilize OGR_STYLE?

2014-01-30 Thread Brian Case
Jukka you could use mapserver to create your image. it will make use of ogrstyle shp2img or if that is an issue use mapserver and the wms driver in gdal brian On Thu, 2014-01-30 at 15:51 +, Jukka Rahkonen wrote: > Hi, > > By reading the manual page and after doing some experiments I guess

Re: [gdal-dev] Question about gdaladdo and raster mosaicking

2014-01-23 Thread Brian Case
Ammar you want to turn on internal masks in the shell, gdaladdo will not use them otherwise. also you will need to turn this on if you use mapserver. gdal will need to have this turned on for any further operations. gdaladdo should use the correct image options without specifying them export GD

Re: [gdal-dev] Reading CurveString and CurvePolygon Geometries from FGF Data

2014-01-21 Thread Brian Case
bob, if you decide to go the mask band route you wont be able to use a vrt. you will need to use gdalwarp or mapserver Brian On Tue, 2014-01-21 at 09:04 -0800, Bob Cave wrote: > Even Rouault wrote > > Le vendredi 17 janvier 2014 17:08:38, Bob Cave a écrit : > >> Hello, > >> > >> I am using the

Re: [gdal-dev] gdalwarp producing rgb 0 0 0 for input nodata areas

2013-02-24 Thread Brian Case
nput-nodata-areas-tp5035602.html > Sent from the GDAL - Dev mailing list archive at Nabble.com. > ___ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Brian Case KF7WPK r...@winkey.org ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] gdal_polygonize and polygon edges

2013-02-19 Thread Brian Case
gt; > ---+-- > I set the clouds in motion - turn up | Frank Warmerdam, > warmer...@pobox.com > light and sound - activate the windows | > http://pobox.com/~warmerdam >

Re: [gdal-dev] substitute pixel values

2013-01-09 Thread Brian Case
Station 145 > 51 Commerce Drive > Augusta, ME 04333-0145 > 69o 47' 58.9"W 44o 21' 54.8"N > > > _______ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Brian Case KF7WPK r...@w

Re: [gdal-dev] Multiple pens in OGR feature style

2012-12-10 Thread Brian Case
w a two-color line, so > it makes sense that LIBKML would not allow this. > ___ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Brian Case KF7WPK r...@winkey.org _

Re: [gdal-dev] Unable to open GRIB file

2012-02-14 Thread Brian Case
On Wed, 2012-02-15 at 09:10 +0300, SupportLists wrote: Carlos is it raster or vector data? can you make the file accessible some place? Brian > Hi, > > I am trying to open a grib file (meteorological data) with GDAL 1.8.1 > but i get the following error: > > ./2011.grib is a grib file, but

Re: [gdal-dev] building geotransform with lower left and upper right coordinates

2012-02-13 Thread Brian Case
Brice, After some thought your image might not be bottom up. if thats the case calc the origin the way you did before, but don't convert back to lat lon Brian On Mon, 2012-02-13 at 05:40 -0700, la...@ucar.edu wrote: > Hello, > > I have a raster where I know the projection and the lower left and

Re: [gdal-dev] building geotransform with lower left and upper right coordinates

2012-02-13 Thread Brian Case
Brice, with that projection the origin should be in meters. convert the ll to meters and try setting dy to dy * -1 Brian On Mon, 2012-02-13 at 05:40 -0700, la...@ucar.edu wrote: > Hello, > > I have a raster where I know the projection and the lower left and upper > right coordinates and pixel s

Re: [gdal-dev] Geolocation Arrays ??

2012-02-12 Thread Brian Case
This could fix when the algorithm fails or generates an extent that is > too large. > > Etienne > > On Sun, Feb 12, 2012 at 12:38 PM, Brian Case wrote: > > the forward transform is pretty strait forwards, however the reverse is > > not, it creates back-mask (a geo-r

Re: [gdal-dev] Geolocation Arrays ??

2012-02-12 Thread Brian Case
egards, > Etienne > > > On Tue, Feb 7, 2012 at 8:46 PM, Brian Case wrote: > > Etienne > > with the transformer the way it sits now, in most cases you will need to > > specify -te xmin ymin xmax ymax with gdalwarp to avoid that error > > > > Brian > &

Re: [gdal-dev] Geolocation Arrays ??

2012-02-11 Thread Brian Case
Etienne with the transformer the way it sits now, in most cases you will need to specify -te xmin ymin xmax ymax with gdalwarp to avoid that error Brian On Tue, 2012-02-07 at 19:22 -0200, Etienne Tourigny wrote: > Even, Frank, > > thanks for your answers. I have run into a problem with using >

Re: [gdal-dev] End User Questions :( KML Raster to Shapefile

2012-01-17 Thread Brian Case
Dennis gdal does not have any code for reading raster files in kml/kmz like that. can you make the files you have available so i can look? Brian On Mon, 2012-01-16 at 17:11 -0600, Dennis Burgess wrote: > I have a .KML File and .PNG file or a .TIF along with a .GEO file that > my application exp

Re: [gdal-dev] Warping onto an image with only GCPs with Python

2012-01-13 Thread Brian Case
Anton, you can get the extents of the modis swath file from the metadata. you have to use the the geolocation arrays to create a grid from the swath file. the further left and right from the center of the image the more area the single pixel in the swath file covers. also modis is a pushbroom sen

Re: [gdal-dev] Warping onto an image with only GCPs with Python

2012-01-12 Thread Brian Case
Anton EOS_SWATH has 2 geolocation arrays, you might try the -geoloc switch to gdalwarp, otherwise it creates a handfull of gcps from the arrays, if you get a error about points failing to transform, set the output extent with the -te switch. also you could try the swath2grid application in https:

Re: [gdal-dev] Motion: Authorize Islandwood Funding for Brian Case

2012-01-09 Thread Brian Case
the Islandwood Code Sprint for Brian Case. > > Folks, > > I declare this motion passed with support from Tamas, Even, > Frank, Daniel and Howard. If there is any other GDAL contributors > interested in attending the sprint under similar arrangements let me > know. The

Re: [gdal-dev] Reading OGRStyle information using C++

2012-01-06 Thread Brian Case
Aneesh i believe the only driver that supports layer andf dataset style tables is the libkml driver Brian On Thu, 2012-01-05 at 18:06 +0530, Aneesh Muralidharan wrote: > Chaithanya, > > Ahhh that explains it ! Thanks a lot. By the way where do I find out > these information (like DXF driver do

Re: [gdal-dev] Writing Pixel Functions with VRT

2011-12-21 Thread Brian Case
Yves, You need to be write a custom application to make use of If you are looking for a command line utility solution you should use gdal_calc.py if your end goal is to serve this data via a wms you could try the thredds server http://www.unidata.ucar.edu/projects/THREDDS/tech/TDS.html http:

Re: [gdal-dev] re-gridding to a coarser grid

2011-12-15 Thread Brian Case
Andreas I think thats referring to my suggestion of GMT to do part of the work On Thu, 2011-12-15 at 14:00 +0100, Andreas H. wrote: > Joaquim, > > > nearneighbor is not really a good thing to use. Grid it with surface. > > What do you mean by "grid it with surface"? In my GDAL (1.8.1), I o

Re: [gdal-dev] Which is the way to handle altitude layers within one file?

2011-12-14 Thread Brian Case
Andreas, you can use gdal to read grib, hdf, and netcdf files built this way. take a look at the format pages for those formats. http://www.gdal.org/formats_list.html note: gdal requires all the "bands" or subdatasets to be the same size if you find yourself needing to put multiple files togethe

Re: [gdal-dev] Newbie question about gdal_grid for 3d data

2011-12-14 Thread Brian Case
Andreas, There are no vector drivers for hdf and nc. take a look here http://www.gdal.org/ogr/ogr_formats.html As an alternative to gdal_grid you can use GMT's blockmean, blockmedian, or blockmode. then use GMT's surface http://www.soest.hawaii.edu/gmt5/ On Wed, 2011-12-14 at 15:03 +0100, Andre

Re: [gdal-dev] how to combine overlapping rasters?

2011-12-06 Thread Brian Case
Etienne you can gdal_rasterise the polys into a alpha band, then add your alpha band to your image in a vrt. Brian On Tue, 2011-12-06 at 19:59 -0200, Etienne Tourigny wrote: > > A related question: How can I effectively mask (mark as nodata) raster > pixels inside a polygon? For example, I wo

Re: [gdal-dev] how to combine overlapping rasters?

2011-12-06 Thread Brian Case
Etienne if they overlap your limited to one of 2 options to mosaic gdalwarp or mapserver to get what ones are on top with gdalwarp i would imagine you would just need to play with the order on the command line. Brian On Tue, 2011-12-06 at 19:59 -0200, Etienne Tourigny wrote: > Hi all. > > I

Re: [gdal-dev] Editing KML files using ogr_csharp.dll

2011-12-05 Thread Brian Case
Chris, Have you tried the libkml driver? http://www.gdal.org/ogr/drv_libkml.html Brian On Mon, 2011-12-05 at 14:50 +, Pouliot, Christopher (DNR) wrote: > Hello All, > > > > I have a C# application that uses OGR_CSHARP.DLL to read/write KML > files. I’m able to read from KML files fin

Re: [gdal-dev] Re: Gdal_translate, question from an amateur.

2011-10-23 Thread Brian Case
e. > > > > > > Thanks > > J > > > > On 23 Oct 2011, at 15:26, "Jay L. [via OSGeo.org]" <[hidden > > email]> wrote: > > > >

Re: [gdal-dev] Re: Gdal_translate, question from an amateur.

2011-10-23 Thread Brian Case
joolek, I noticed you used the word "strip". Are you trying to remove portions of the image that are all black or white? On Sun, 2011-10-23 at 03:01 -0700, joolek wrote: > Basically - I have an imagery with RGB 0-255, but I would like to strip it to > 1-254. > So i want to get rid off two extrem

Re: [gdal-dev] gdal_fillnodata.py filling interior holes in RGB images

2011-10-12 Thread Brian Case
Travis, you have to run gdal_fillnodata on each band, then splice the outputs back together. if the nodata areas on the inside do not connect to the outside area you could create an alpha band ### -color 0,0,0 assumes odata value of 0 nearblack -dstalpha -near 0 -nb 0 -setalpha -color 0,0,0 \

Re: [gdal-dev] Re: KML's extendedData not visible

2011-10-11 Thread Brian Case
joss, you should try the libkml driver. http://www.gdal.org/ogr/drv_libkml.html ogr2ogr -f CSV test.csv test.kml cat test.csv Name,description,timestamp,begin,end,altitudeMode,tessellate,extrude,visibility,holeNumber,holeYardage,holePar Club house,,-1,-1,-1,1,234,4 By the lake,,-1,-1,-1,

Re: [gdal-dev] Resquest for comments (RFC 36)

2011-10-03 Thread Brian Case
Ivan Makes me wonder what kind of performance gain you could see with mapserver on a layer with a large number of images to open. Seems to me that a tileindex would be a good application of this. Brian On Mon, 2011-10-03 at 19:02 -0500, Ivan Lucena wrote: > Hi All, > > RFC 36: http://trac.o

Re: [gdal-dev] Re: Can't get gdalwarp to create a raster in Mercator-2SP projection?

2011-09-29 Thread Brian Case
Patrick from http://remotesensing.org/geotiff/proj_list/mercator_2sp.html "The difference between 1SP and 2SP formulations is that the initial parameters are different: in 1SP you start from the longitude of natural origin and the scale factor, while in 2SP you start from the longitude of natural

Re: [gdal-dev] Can't find a way of properly warping into a buffer

2011-09-12 Thread Brian Case
Daniele a simple solution could be warp your image to an in memory raster http://www.gdal.org/frmt_mem.html Brian On Mon, 2011-09-12 at 19:24 +0200, Daniele Tessaro wrote: > Hello everybody, > I'm really sorry if this post will be very long, but I swear I've > spent a lot of time (days!) > in t

Re: [gdal-dev] Re: discussion on improvements to the NetCDF driver and CF-1 convention

2011-09-02 Thread Brian Case
in the wiki entry at > http://trac.osgeo.org/gdal/wiki/NetCDF_Improvements . > > Etienne > > > > > - Original Message - > From: Brian Case > To: Etienne Tourigny >

Re: [gdal-dev] RE: GDALWARP Chunk Problem (Tiling Effect)

2011-09-01 Thread Brian Case
Chris have you tryed other kernels? cubic, etc.. read this too http://www.gdal.org/structGDALWarpOptions.html#0ed77f9917bb96c7a9aabd73d4d06e08 perhaps someone else may have other ideas. Brian On Wed, 2011-08-31 at 23:58 -0700, heng.feng wrote: > Tried the "-et 0.0" option with no luck :( > >

Re: [gdal-dev] GDALWARP Chunk Problem (Tiling Effect)

2011-08-31 Thread Brian Case
Chris, Try passing gdalwarp -et 0.0 I am also suspicious of doing compression with the warp. you might want to do that as a second step Brian On Wed, 2011-08-31 at 23:16 -0700, heng.feng wrote: > Hi, > > I have used the (OSGeo4W) GDALWARP command to warp some height data image > files (NTv2 tr

Re: [gdal-dev] Re: discussion on improvements to the NetCDF driver and CF-1 convention

2011-08-28 Thread Brian Case
any thoughts on a ogr netcdf driver? I have a simple one for madis data that needs work in my sandbox Brian On Sun, 2011-08-28 at 14:56 -0700, Etienne Tourigny wrote: > I have created a new wiki page at > http://trac.osgeo.org/gdal/wiki/NetCDF_Improvements , following Even's > suggestion (Tha

Re: [gdal-dev] pruge images

2011-08-26 Thread Brian Case
? http://www.gdal.org/gdaladdo.html Brian On Fri, 2011-08-26 at 15:53 +0200, Jan Tappenbeck wrote: > > Am 26.08.2011 15:48, schrieb Brian Case: > > Jan > > > > if you already have a wms setup with your data you could use the wms > > driver to read it. gda

Re: [gdal-dev] pruge images

2011-08-26 Thread Brian Case
Jan if you already have a wms setup with your data you could use the wms driver to read it. gdal_translate to a bigtiff. then cut it up into tiles with gdal_translate -srcwin Brian On Fri, 2011-08-26 at 15:17 +0200, Jan Tappenbeck wrote: > > hi! > > i have following problem. openstreetmap c

Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff

2011-07-26 Thread Brian Case
with 150 files? > > --Nikos > > On Tue, Jul 26, 2011 at 9:17 AM, Brian Case wrote: > Nikolaos > > I have warped modis with swath2grid and mosaicked it with > mapserver, I > have not seen that issue. > > also y

Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff

2011-07-26 Thread Brian Case
transform them first and mosaic them > after you might have > white lines. > > --Nikos > > 2011/7/25 Brian Case > anna, Nikolaos, > > after much experimentation I found something that seems to > work. > > gdalwarp

Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff

2011-07-25 Thread Brian Case
anna, Nikolaos, after much experimentation I found something that seems to work. gdalwarp --config GEOL_AS_GCPS FULL \ -tps -geoloc -t_srs EPSG:4326 \ -te -119.88004906118 13.5092548294379 \ -92.6036641645677 34.7070050511582 \ HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011169.1750

Re: [gdal-dev] Draw map in Visual C++

2011-07-22 Thread Brian Case
mark That is a pretty general question so difficult to give you an answer. I would like to point out that gdal will not draw the maps in your application but will let you read the data from many different formats with a single api. Here is the urls to the tutorials. raster data: http://gdal.org

Re: [gdal-dev] [gdal-1.8.0] Get new Coordinates from Geotiff avoiding black area

2011-07-12 Thread Brian Case
Luiz, do you have a url handy for the source code of that plugin? Brian On Tue, 2011-07-12 at 14:53 -0300, Luiz Motta wrote: > Met, > > I made a QGIS's plugin where calculate the 4 corners of polygon with > the real pixels inside. Maybe is useful for you. > The name of plugin is "Image Boundary

Re: [gdal-dev] Unwanted partial transparency when clipping

2011-07-07 Thread Brian Case
Michael, another option is to use gdalwarp to merge images together. Brian On Thu, 2011-07-07 at 09:47 -0700, Michael Corey wrote: > Thanks for your responses, Chaitanya and Eli. The semitransparency isn't > really noticeable until you put the image on top of something else (in > Photoshop or

Re: [gdal-dev] KML 2 SHP with attributes

2011-06-30 Thread Brian Case
Patrick, Can you give an example of what tags your loosing? are you using the kml or the libkml driver? Brian On Thu, 2011-06-30 at 15:30 -0600, Patrick Wild wrote: > I am converting KML to .shp with OGR2ORG, however, due to the KML tags > I am loosing my attributes. Does anyone have any experie

Re: [gdal-dev] jpeg in tiff artifacts

2011-06-16 Thread Brian Case
matt, jpeg is a lossy format. the nodata is ignored because those pixels are no longer the same value as the nodata value. You can add a mask band to the tiff but unfortunately aside from gdalwarp and mapserver I don't think there is any other way to read it. see http://docs.oam.osgeo.org/storag

Re: [gdal-dev] Converting from 16 bit to 8 bit images (best practices?)

2011-06-09 Thread Brian Case
Nikos 143, and 721 Brian On Thu, 2011-06-09 at 16:49 -0700, Nikolaos Hatzopoulos wrote: > for modis truecolor images? > > --Nikos > > On Thu, Jun 9, 2011 at 4:44 PM, Brian Case wrote: > Jonathan > > I don't know about worldview, bu

Re: [gdal-dev] Converting from 16 bit to 8 bit images (best practices?)

2011-06-09 Thread Brian Case
ng applications need 3-band data, usually with > only 8 bits per band. > > > Pushing forward, but if someone knows how to figure out an optimal set > of scale values to use within GDAL, please post the solution. > > > --j > > On Wed, Jun 8, 2011 at 9:45 PM

Re: [gdal-dev] Converting from 16 bit to 8 bit images (best practices?)

2011-06-08 Thread Brian Case
Jonathan you can do scaling on individual bands. gdal_translate -b 1 -scale 0 5000 -ot Byte -of VRT infile outfile1.vrt gdal_translate -b 2 -scale 0 4000 -ot Byte -of VRT infile outfile2.vrt gdal_translate -b 3 -scale 0 3000 -ot Byte -of VRT infile outfile3.vrt gdalbuildvrt -separate output.vrt

Re: [gdal-dev] Adding libKML

2011-05-18 Thread Brian Case
Paul You need libkml 1.3.0. this was only svn but i do not know if they have released since i last looked.\ Brian On Wed, 2011-05-18 at 11:35 +0200, Paul Meems wrote: > Hi list, > > I'm still struggling enabling new (at least for us) support for > certain file formats. > I'm working on Win7 wit

[gdal-dev] gdal internal vrs external libs and gdal-config

2011-05-14 Thread Brian Case
hi, This week I ran into an issue building OTB against a gdal with internal libtiff and libgeotiff. After spending a bit of time trying to figure out why i had a problem I thought it might be useful if gdal-config could report on internal libs and if they are hidden. What came up with is a gdal-

Re: [gdal-dev] gdal reading in grib values in awrong state?

2011-05-05 Thread Brian Case
ome of our own grib2 files that may > at times have unitless values. For some odd reason they are being converted > to Celsius as well. > > Sent from my iPad > > On May 5, 2011, at 5:09 PM, "Brian Case" wrote: > > > Matt, > > > > At a quic

Re: [gdal-dev] gdal reading in grib values in awrong state?

2011-05-05 Thread Brian Case
Matt, At a quick glance I would say pygrib is outputting in kelvin and gdal Celsius. also it looks like there not in the same order. Brian On Thu, 2011-05-05 at 14:57 -0600, maf...@nmsu.edu wrote: > Hi, > i am a little confused. I am using the gdal/python and trying to read > values from a gri

Re: [gdal-dev] Re: cropping a tif with lat/lon

2011-04-20 Thread Brian Case
Dan use nearblack to add a alpha or mask band, and collar the image Brian On Wed, 2011-04-20 at 19:15 -0700, spalt wrote: > Hi Brian- > Thanks so much for your help. That works, though I feel a little silly now > because I am realizing the lat/lon the FAA provides with these is the > corners o