Re: [gdal-dev] How do I add a projection to proj 8?

2024-04-13 Thread Stephen Woodbridge via gdal-dev
On 4/13/2024 1:26 PM, Javier Jimenez Shaw wrote: On Sat, 13 Apr 2024 at 17:35, Stephen Woodbridge via gdal-dev wrote: Thanks, this is NOT the standard Web Mercator projection. I am aware of EPSG:90013 and EPSG:3857. This projection is used with HYCOM data that I have extracted

Re: [gdal-dev] How do I add a projection to proj 8?

2024-04-13 Thread Stephen Woodbridge via gdal-dev
something like EPSG:90014? Thanks,   -Steve On Sat, 13 Apr 2024 at 06:17, Sebastiaan Couwenberg via gdal-dev wrote: On 4/12/24 11:24 PM, Stephen Woodbridge via gdal-dev wrote: > and was able to access it in gdal, mapserver, postgis, etc with > "EPSG:900914"

[gdal-dev] How do I add a projection to proj 8?

2024-04-12 Thread Stephen Woodbridge via gdal-dev
Hi all, I've been gone for a while, but got called back to update a site I built and need to move from proj4 to proj 8 on Ubuntu 22.04. In the past I just added the following to /usr/share/proj/epsg # HYCOM Mercator projection <900914> +proj=merc +a=6371001 +b=6371001 +lat_ts=0.0 +lon_0=0.0

Re: [gdal-dev] Considering drivers removal ?

2021-01-27 Thread Stephen Woodbridge
I think a 4th option is a hybrid approach of moving to a more modular plug-in architecture that allows the core more flexibility to evolve at the same time by moving to a more plug-in driver allows for more independent development, testing and release lets the community participation. This does

Re: [gdal-dev] Driver maintenance - long-term solution ?

2021-01-15 Thread Stephen Woodbridge
So reading through this thread, my cynical side agrees with the "this project is dying from lack of funding, ..." approach, but I'm not sure that works for the long run as people get tired of hear the "sky is falling" and ignore it over time. A different approach could be something along the

Re: [gdal-dev] Considering drivers removal ?

2021-01-12 Thread Stephen Woodbridge
On 1/12/2021 5:36 PM, Even Rouault wrote: On mardi 12 janvier 2021 12:56:13 CET Frank Warmerdam wrote: On Tue, Jan 12, 2021 at 12:38 PM Howard Butler wrote: The only question that matters here is "Who is going to maintain it?" and if the answer to that is "no one", it should be removed. There

Re: [gdal-dev] Considering drivers removal ?

2021-01-10 Thread Stephen Woodbridge
Even, This makes a lot of sense to me. How would you handle this in Python? Would it make sense to create a GDAL-removed repository and move stuff into it just so it is available if someone wants it. This would not be supported or updated by GDAL just making it available if someone

Re: [gdal-dev] Contour Line Thinning

2021-01-01 Thread Stephen Woodbridge
On 1/1/2021 9:47 PM, Richard Greenwood wrote: On Fri, Jan 1, 2021 at 2:36 PM Stephen Woodbridge mailto:stephenwoodbridg...@gmail.com>> wrote: Hi all, I'm contouring bathemetry data using gdal_contour and it works really great. The problem I have is that when depth fal

[gdal-dev] Contour Line Thinning

2021-01-01 Thread Stephen Woodbridge
Hi all, I'm contouring bathemetry data using gdal_contour and it works really great. The problem I have is that when depth falls off rapidly like at the continental shelf or into a canyon, I get too many contour lines that all bunch up. If I change the contour step size to fix this, then the

Re: [gdal-dev] Problem writing with ODS and XLSX drivers

2020-09-12 Thread Stephen Woodbridge
, Stephen Woodbridge wrote: Try exporting to CSV and check the that the records do not have any embeded or chars. One simple way to do this is to compare the record count of the shp file to the line count of the CSV file. Try importing the CSV file into EXCEL  or Libreoffice and see

Re: [gdal-dev] Problem writing with ODS and XLSX drivers

2020-09-12 Thread Stephen Woodbridge
Try exporting to CSV and check the that the records do not have any embeded or chars. One simple way to do this is to compare the record count of the shp file to the line count of the CSV file. Try importing the CSV file into EXCEL  or Libreoffice and see if it reports a problem. This may

Re: [gdal-dev] I have a problem getting higher resolution in a VRT file

2020-08-20 Thread Stephen Woodbridge
downsampling-geotiff-with-a-low-pass-filter-td5385890.html. In this gis.SE question cubic spline was considered good https://gis.stackexchange.com/questions/30627/smoothing-reinterpolating-raster-with-gdal. -Jukka- -Alkuperäinen viesti- Lähettäjä: Stephen Woodbridge Lähetetty: torstai 20.

[gdal-dev] I have a problem getting higher resolution in a VRT file

2020-08-19 Thread Stephen Woodbridge
Hi all, I've been puzzling over a problem with trying to get a higher resolution images using a VRT file so low res data looks smoother in mapserver. This is my source tif image: gdalinfo  /maps/wms/data/HYCOM/HYCOM_today_mlt.tif Driver: GTiff/GeoTIFF Files:

Re: [gdal-dev] Dealing with "default" TOWGS84 values

2020-03-31 Thread Stephen Woodbridge
Even, Would it make sense to add an option to GDAL like -ignore-bound-towgs84 or -use-bound-towgs84 depending on what the default behavior is set to? This would allow the user to choose which is best for their data. -Steve W On 3/31/2020 10:21 AM, Even Rouault wrote: Hi, I wanted to get

[gdal-dev] Solved: Re: Having a problem getting VRT ComplexSource to Scale

2020-01-23 Thread Stephen Woodbridge
Never mind! I just checked the histogram and the data is getting scaled, I forgot that I was also scaling it in my mapfile so the results  where getting masked. -Steve On 1/23/2020 8:20 PM, Stephen Woodbridge wrote: Hi all, I'm having trouble getting my VRT file to scale the source data

[gdal-dev] Having a problem getting VRT ComplexSource to Scale

2020-01-23 Thread Stephen Woodbridge
Hi all, I'm having trouble getting my VRT file to scale the source data into my ColorTable. I've been looking at this https://gdal.org/drivers/raster/vrt.html The source GeoTiff has: Band 1 Block=4500x1 Type=Float32, ColorInterp=Gray   Minimum=1.082, Maximum=46.322, Mean=34.181, StdDev=2.049

Re: [gdal-dev] ogr2ogr PDF

2019-11-14 Thread Stephen Woodbridge
Paul, Your issue might be one of scale. You have specified units in px try changing it to pt. (units (g, px, pt, mm, cm, in)) See: https://gdal.org/user/ogr_feature_style.html#ogr-feature-style You can also try named pens like: Here is the current list of OGR pen ids (this could grow over

Re: [gdal-dev] Question on gdal_contour and limiting the range of contours generated.

2019-10-31 Thread Stephen Woodbridge
Thanks I have not seen the raster common option file thing. I’ll look at it. Sent from my iPhone > On Oct 31, 2019, at 2:11 PM, jratike80 > wrote: > > Hi, > > One more suggestion, have you tried the -fl option "Name one or more “fixed > levels” to extract"? Requires some writing but only

Re: [gdal-dev] Question on gdal_contour and limiting the range of contours generated.

2019-10-31 Thread Stephen Woodbridge
Jukka, Thanks! After trying shapefile and exceeding 2GB, I switched to postgres and deleted the elevation as you suggested. I'm working with a global relief model the includes both elevation and bathymetry depth over the whole world at about 15 arc-sec/pixel. I generates a lot of contour

[gdal-dev] Question on gdal_contour and limiting the range of contours generated.

2019-10-30 Thread Stephen Woodbridge
Hi, I want to use gdal_contour to generate shapefile contours which is pretty straight forward except I'm working with coastal relief data that has both elevation and bathymetry. I only want contours for the bathymetry. I tried setting the -off -5 and -i -5 with the idea that this would

[gdal-dev] Problems with masking while create a color-relief from a DEM

2019-10-24 Thread Stephen Woodbridge
Hi, I have a collection of region netCDF DEM files that overlap and have a nodata=-nan. I have been able to convert them to color-relief using gdaldem but somewhere in the processing chain I seem to be loosing the nodata/mask/alpha info and the nodata in the color-relief is opaque, ie: where

Re: [gdal-dev] Question on converting netCDF to GTiff

2019-10-23 Thread Stephen Woodbridge
On 10/23/2019 1:54 PM, Even Rouault wrote: On mercredi 23 octobre 2019 13:41:45 CEST Stephen Woodbridge wrote: Hi, I have a netCDF grid of Float32 values that I want to convert to GTiff with Int16 values. The min/max values will support this without scaling, but I'm not sure what will happen

[gdal-dev] Question on converting netCDF to GTiff

2019-10-23 Thread Stephen Woodbridge
Hi, I have a netCDF grid of Float32 values that I want to convert to GTiff with Int16 values. The min/max values will support this without scaling, but I'm not sure what will happen to the NODATA value = 9.96920996838686905e+36. Ideally, I would like to set it to something like 32000 but

Re: [gdal-dev] vrt guidance: 31 band classified raster

2019-10-03 Thread Stephen Woodbridge
Matt, Have you tried create a cascading VRT, where the first VRT composites all the GTiff files into a single virtual raster and then reference the composite in another VRT that assigns colortable. I think this should work nicely. -Steve W On 10/3/2019 7:01 PM, matt.wil...@gov.yk.ca wrote:

Re: [gdal-dev] Is there an easy way to clip an image to the realdata?

2019-06-27 Thread Stephen Woodbridge
how to tackle your use case. The full script can be found here: https://gist.github.com/LarsSchy/9ecb31eb964dd83820c139b2f2769a7c Have fun Lars Schylberg 20 juni 2019 kl. 22:06, "Stephen Woodbridge" skrev: Hi, I'm working with VIIRS L3U images of L2P (level 2 patches). The L3U

Re: [gdal-dev] adding features to a layer/dataset...

2019-06-26 Thread Stephen Woodbridge
Shayne, To solve your problem, given the constraints the Even explained, you probably need to do something like: * program that is writing to shapefile    - creates a socket or pipe    - open the file for exclusive access    - after it updates the shapefile    - close the file    - writes to

Re: [gdal-dev] Help with band pixel functions

2019-06-24 Thread Stephen Woodbridge
Hi Adriana, Here are a few ideas: import numpy as np ima = gdal.Open(filea, 0) banda = ima.GetRasterBand(1) dataa = banda.ReadAsArray() print type(dataa), dataa.shape, isinstance(dataa, np.ma.MaskedArray) imb = gdal.Open(fileb, 0) bandb = imb.GetRasterBand(1) datab = bandb.ReadAsArray()

Re: [gdal-dev] How do I get lat long extents from this Ogrinfo output on a .shp?

2019-06-20 Thread Stephen Woodbridge
Hi Paul, Its not just the extents but all the data in the file is in a UTM projection. You can reproject the file to WGS84 and then ogrinfo will report in long-lat. ogr2ogr -t_srs EPSG:4326 OldStreamsPolyline-wgs84.shp OldStreamsPolyline.shp ogrinfo -so OldStreamsPolyline-wgs84.shp

[gdal-dev] Is there an easy way to clip an image to the realdata?

2019-06-20 Thread Stephen Woodbridge
Hi, I'm working with VIIRS L3U images of L2P (level 2 patches). The L3U data is gridded to the +-180 x +-90 but the patch only fills a small percentage of that. I'm compositing a the patch for a day using a vrt. Performance is an issue because all the patches cover the whole globe instead

[gdal-dev] Stuck with python gdal.Grid throwing an error

2019-06-14 Thread Stephen Woodbridge
Hi, I can't figure out what I'm doing wrong but gdal.Grid() keeps throwing an error and I've run out of things to try and work around the issue. On Ubuntu using these packages: gdal-bin:amd64/xenial 2.2.2+dfsg-1~xenial1 uptodate gdal-data:all/xenial 2.2.2+dfsg-1~xenial1 uptodate

Re: [gdal-dev] Question on python gdal.Rasterize()

2019-06-14 Thread Stephen Woodbridge
OK, After some more research, I found the rasterize.py autotest file which was a huge help and I have code working that burns the tracks into an image. But its not exactly what I hoping for so I'll give the gdal.Grid() a try. Thanks,   -Steve On 6/14/2019 1:58 PM, Stephen Woodbridge wrote

[gdal-dev] Question on python gdal.Rasterize()

2019-06-14 Thread Stephen Woodbridge
Hi all, My goal is to take satellite track data and create a gtiff file using Python. The satellite data is in a NetCDF file, which I can read in Python and has variables lat, lon, ssha. There are a continue stream of the NetCDF data over time so I plan to just keep loading them as they

Re: [gdal-dev] Question of file ordering in VRT

2019-03-07 Thread Stephen Woodbridge
Thanks Even and Ivan, Good idea to build to build the list and use that for gdalbuildvrt. -Steve On 3/7/2019 11:42 AM, Even Rouault wrote: On jeudi 7 mars 2019 11:30:04 CET Stephen Woodbridge wrote: Hi all, When I have multiple overlapping files in VRT which pixel takes president? I'm

[gdal-dev] Question of file ordering in VRT

2019-03-07 Thread Stephen Woodbridge
Hi all, When I have multiple overlapping files in VRT which pixel takes president? I'm assuming that the later files in the VRT would be the one presented. For example, I have multiple satellite images with names like MMDDHHMM-.tif so when I use gdalbuiltvrt MMDD.vrt

[gdal-dev] NOAA satellite swath NetCDF file to georeferenced Gtiff?

2019-03-05 Thread Stephen Woodbridge
Hi All, Can someone please give me a hint on how to convert a NOAA satellite swath NetCDF file to georeferenced Gtiff? I found these:   https://trac.osgeo.org/gdal/ticket/4513   https://trac.osgeo.org/gdal/wiki/rfc4_geolocate which implies it cannot be done directly with gdal yet. Is there

Re: [gdal-dev] gdalwarp ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)

2018-07-01 Thread Stephen Woodbridge
Corps of Engineers On Jul 1, 2018, at 10:07 AM, Stephen Woodbridge wrote: Hi all, I'm trying to convert a utm jp2 file to a wgs84 tif file and getting "ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)" but the src file has only 3 RBG bands. gdalwarp -

Re: [gdal-dev] gdalwarp ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)

2018-07-01 Thread Stephen Woodbridge
of Engineers On Jul 1, 2018, at 10:07 AM, Stephen Woodbridge wrote: Hi all, I'm trying to convert a utm jp2 file to a wgs84 tif file and getting "ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)" but the src file has only 3 RBG bands. gdalwarp -t_srs EPSG:4326 -d

[gdal-dev] gdalwarp ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)

2018-07-01 Thread Stephen Woodbridge
Hi all, I'm trying to convert a utm jp2 file to a wgs84 tif file and getting "ERROR 6: PHOTOMETRIC=YCBCR requires a source raster with only 3 bands (RGB)" but the src file has only 3 RBG bands. gdalwarp -t_srs EPSG:4326 -dstalpha -r bilinear -of GTiff -co BIGTIFF=YES -co TILED=YES -co

Re: [gdal-dev] Having trouble working with jp2 files

2018-06-20 Thread Stephen Woodbridge
20 juin 2018 12:33:40 CEST Stephen Woodbridge wrote: Even, Thanks! I did built gdal 2.2.4 with openjpeg 2.1.2 on a docker container and that solves the simple case of: gdalinfo -stats /data/850011.jp2 but it still fails with: root@436a09289084:/usr/local# gdal_translate -of GTiff -co TILED=YES

Re: [gdal-dev] Having trouble working with jp2 files

2018-06-20 Thread Stephen Woodbridge
iff -co TILED=YES /data/850011.jp2 /data/850011.tif Input file size is 29364, 25856 0ERROR 1: psImage->comps[0].data == NULL ERROR 1: /data/850011.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 On 6/20/2018 10:27 AM, E

Re: [gdal-dev] Having trouble working with jp2 files

2018-06-19 Thread Stephen Woodbridge
And got the same errors. Even, I sent you a link to download the file from my dropbox. -Steve On 6/19/2018 5:11 PM, Even Rouault wrote: On mardi 19 juin 2018 16:57:30 CEST Stephen Woodbridge wrote: Hi all, I'm having a problem readding jp2 files. My goal is to be able to display them through

[gdal-dev] Having trouble working with jp2 files

2018-06-19 Thread Stephen Woodbridge
Hi all, I'm having a problem readding jp2 files. My goal is to be able to display them through mapserver, but I'm have problems with gdalinfo and gdal_translate. I set up a docker image with the latest gdal2 in it like the following.  I have 1074 jp2 files and all that I have tried are

[gdal-dev] First attempt to use raster - need help

2018-01-13 Thread Stephen Woodbridge
Hi all, I'm trying to use the raster functions for the first time. I have a simple black and white image with some polygons drawn on it. and I want to extract the polygons as geometry. $ gdalinfo -hist 12322240_lines.png Driver: PNG/Portable Network Graphics Files: 12322240_lines.png   

Re: [gdal-dev] ogr2ogr PROMOTE_TO_MULTI with linestrings?

2017-05-16 Thread Stephen Woodbridge
of the Census roads data has a mix of Linestring and MultiLinestring features, which is why I used PROMOTE_TO_MULTI in the first place. -Steve On 5/16/2017 2:58 AM, jratike80 wrote: Stephen Woodbridge wrote On 5/15/2017 6:15 PM, Stephen Woodbridge wrote: Hi, I'm trying to load both poly

Re: [gdal-dev] ogr2ogr PROMOTE_TO_MULTI with linestrings?

2017-05-15 Thread Stephen Woodbridge
On 5/15/2017 6:15 PM, Stephen Woodbridge wrote: Hi, I'm trying to load both polygons and linestrings and would like PROMOTE_TO_MULTI to work for both, but it appears to define Multipolygon type and does not work for linestrings/multilinestrings. Maybe it would make more sense to have

[gdal-dev] ogr2ogr PROMOTE_TO_MULTI with linestrings?

2017-05-15 Thread Stephen Woodbridge
Hi, I'm trying to load both polygons and linestrings and would like PROMOTE_TO_MULTI to work for both, but it appears to define Multipolygon type and does not work for linestrings/multilinestrings. Maybe it would make more sense to have: PROMOTE_TO_MULTI - work for either polygon or

Re: [gdal-dev] A couple of new issues with gdal

2017-04-14 Thread Stephen Woodbridge
of imagery. For example California has something like 10-11,000 DOQQs. Anyway, let me know if you need more info. Thanks, -Steve On 4/14/2017 2:06 PM, Even Rouault wrote: On vendredi 14 avril 2017 13:39:12 CEST Stephen Woodbridge wrote: > Hi Even, > > GDAL 2.1.0, released 2

[gdal-dev] A couple of new issues with gdal

2017-04-14 Thread Stephen Woodbridge
Hi Even, GDAL 2.1.0, released 2016/04/25 I'm pretty happy with the new workflow so far, but have te following minor issues. The color interp is still problematic with gdalinfo, even though the *.aux.xml file looks correct, except it is missing the Alpha band. I generated the final tif

Re: [gdal-dev] What is the best strategy for handling multiple banded tiffs

2017-04-13 Thread Stephen Woodbridge
erp=Blue Mask Flags: PER_DATASET Band 4 Block=256x256 Type=Byte, ColorInterp=Gray Mask Flags: PER_DATASET Band 5 Block=256x256 Type=Byte, ColorInterp=Alpha Mask Flags: PER_DATASET And the .aux.xml is the same as above. On 4/12/2017 5:40 PM, Stephen Woodbridge wrote: I'm working wit

Re: [gdal-dev] What is the best strategy for handling multiple banded tiffs

2017-04-12 Thread Stephen Woodbridge
to do that? Or is the best I can do is to create separate RGBA file, an IR file and a sobel file and bring them together via a VRT file? Thanks, -Steve On 4/12/2017 2:01 PM, Even Rouault wrote: On mercredi 12 avril 2017 13:23:06 CEST Stephen Woodbridge wrote: > Hi all, > > I'm

[gdal-dev] What is the best strategy for handling multiple banded tiffs

2017-04-12 Thread Stephen Woodbridge
Hi all, I'm reworking my code dealing with NAIP imagery that has R, G, B, IR bands, and I generate a mask band when reprojecting it. I also have the option to generate another computed band based on a sobel operator. In the past, I separated them R, G, B, mask into jpeg ycbcr compressed

Re: [gdal-dev] ogr2ogr intersect very slow at command line with large sqlite input

2017-04-10 Thread Stephen Woodbridge
Chris, I'm not suggesting that you change your workflow yet, only that you try some experiments. I've worked with both Spatialite and postgis, but I have much more experience with postgis. While they both have many of same functions, the interactions with the indexes is much more automatic

Re: [gdal-dev] ogr2ogr intersect very slow at command line with large sqlite input

2017-04-10 Thread Stephen Woodbridge
On 4/10/2017 5:51 PM, CTL101 wrote: I've removed the intersection part after realising that I didn't need it as the intersection geometry is already calculated in the input (slow hand clap me!) So after investigating a bit more I have: ogr2ogr -gt unlimited -nlt PROMOTE_TO_MULTI25D --config

Re: [gdal-dev] I can't get ogr2ogr to overwrite an existing table

2017-04-08 Thread Stephen Woodbridge
On 4/8/2017 12:29 PM, Even Rouault wrote: On vendredi 7 avril 2017 15:59:49 CEST Stephen Woodbridge wrote: > Hi all, > > I'm having trouble figuring out how to overwrite an existing table with > ogr2ogr. I have tried various combos of arguments without success. > &g

[gdal-dev] I can't get ogr2ogr to overwrite an existing table

2017-04-07 Thread Stephen Woodbridge
Hi all, I'm having trouble figuring out how to overwrite an existing table with ogr2ogr. I have tried various combos of arguments without success. What am I doing wrong? ogr2ogr -t_srs EPSG:4326 -nln roads -nlt PROMOTE_TO_MULTI -f PostgreSQL -overwrite -lco OVERWRITE=YES -lco PRECISION=NO

Re: [gdal-dev] Is it possible to add new fields to a shapefile and populate then?

2017-04-03 Thread Stephen Woodbridge
On 4/3/2017 11:37 PM, jratike80 wrote: Stephen Woodbridge wrote Hi All, I have a shapefile(s) and I want to read the features, generate some metrics about each feature and then add them to the that feature. I'm using python and one obvious way to to in effect create a new shapefile

[gdal-dev] Is it possible to add new fields to a shapefile and populate then?

2017-04-03 Thread Stephen Woodbridge
Hi All, I have a shapefile(s) and I want to read the features, generate some metrics about each feature and then add them to the that feature. I'm using python and one obvious way to to in effect create a new shapefile with the columns I need to add, then copy the existing shapefile to the

Re: [gdal-dev] Getting black lines between doqqs

2017-03-06 Thread Stephen Woodbridge
Thank you Even! that was the trick and it works great now. On 3/6/2017 1:35 PM, Even Rouault wrote: On lundi 6 mars 2017 13:20:50 CET Stephen Woodbridge wrote: Hi, I have a collection of doqqs that look like the two below. I trying to extract a small image the crosses between

[gdal-dev] Getting black lines between doqqs

2017-03-06 Thread Stephen Woodbridge
Hi, I have a collection of doqqs that look like the two below. I trying to extract a small image the crosses between the two doqqs, but I'm getting a black line between them. When I serve them via mapserver, I'm not getting this. gdalbuildvrt -addalpha tmp-tc-13110.vrt

[gdal-dev] Help working with hdf5 images

2017-02-19 Thread Stephen Woodbridge
Hi all, I'm trying to create GTiff files from IMARS Modis hdf5 images. I did something like this about 3-4 years ago when they were using hdf4 datasets but things have changed since then. So I have two questions: 1. How do I extract the sst data (subdataset 3) into a GTiff and create a mask

Re: [gdal-dev] Ubuntu GDAL 1.10.1 support for ECW formats

2017-02-16 Thread Stephen Woodbridge
On 2/16/2017 3:04 PM, Sebastiaan Couwenberg wrote: On 02/16/2017 08:46 PM, Stephen Woodbridge wrote: Any ideas on this? ECW support is not enabled in the Debian package, that's why your client doesn't have it. You probably have a custom build installed in /usr/local or elsewhere. Yes, I

[gdal-dev] Ubuntu GDAL 1.10.1 support for ECW formats

2017-02-16 Thread Stephen Woodbridge
Hi all, On my system: Ubuntu 14.04.5 LTS I have installed: gdal-bin:amd64/trusty 1.10.1+dfsg-5ubuntu1 upgradeable to 2.1.0+dfsg-1~trusty2 libgdal-dev:amd64/trusty 1.10.1+dfsg-5ubuntu1 upgradeable to 2.1.0+dfsg-1~trusty2 libgdal1h:amd64/trusty 1.10.1+dfsg-5ubuntu1 uptodate $ gdalinfo

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

2017-02-13 Thread Stephen Woodbridge
: Stephen Woodbridge wrote Even, Thanks for the quick feedback. I'll will give these suggestions and try today and let you know how it goes. Thanks, -Steve Hi Steve, They are probably good suggestions but folks did not bother to think what is your ultimate target. It is not to improve you

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

2017-02-12 Thread Stephen Woodbridge
On 2/12/2017 9:54 AM, Even Rouault wrote: On samedi 11 février 2017 18:18:35 CET Stephen Woodbridge wrote: Hi All, I need your wisdom. I'm downloading NAIP DOQQs in GTiff format and I have a processing chain something like the following: gdalwarp -t_srs EPSG:4326 -dstalpha -r

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

2017-02-11 Thread Stephen Woodbridge
Hi All, I need your wisdom. I'm downloading NAIP DOQQs in GTiff format and I have a processing chain something like the following: gdalwarp -t_srs EPSG:4326 -dstalpha -r "bilinear -multi -co TILED=YES -dstnodata '0 0 0' srctiff tmpfile nearblack -nb 15 -q tmpfile gdal_translate -co

[gdal-dev] Question on writing a GTiff using python

2016-06-18 Thread Stephen Woodbridge
Hi all, I'm writing a GTiff using GDAL/Python and it is mostly working. The source data is coming from a NETCDF file that is georeferenced in Longitude from 74.1600037 to 434.160 degrees If I use this I get: Size is 4500, 1782 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984",

Re: [gdal-dev] Questions on working with HYCOM data in netcdf format

2016-06-06 Thread Stephen Woodbridge
#file-mapfile-for-hycom-ocean-currents-map Thanks for the help! -Steve On 6/6/2016 2:29 PM, Stephen Woodbridge wrote: On 6/6/2016 1:25 PM, Kurt Schwehr wrote: +gdal-dev Some hackish test code I had laying around from last year... https://gist.github.com/schwehr/01f6604afc7757ea0a676f0eb28be582

Re: [gdal-dev] Questions on working with HYCOM data in netcdf format

2016-06-06 Thread Stephen Woodbridge
heavily used mapserver, but you could probably do the same there. https://www.flickr.com/photos/schwehr/26895009754/ Very cool picture! This really helps to get me started. Thanks, -Steve On Mon, Jun 6, 2016 at 7:51 AM, Stephen Woodbridge <wood...@swoodbridge.com <mailto:wood...@swoodbrid

Re: [gdal-dev] Questions on working with HYCOM data in netcdf format

2016-06-05 Thread Stephen Woodbridge
Hi Even, Kurt, Thank you for the suggestions. I'm reading through stuff and searching for solutions. I need access to the 2d layer which is only available via the GLBa0.08 dataset and not in the GLBu0.08 dataset which is easier to deal with. I'll start with the GLBu0.08 dataset for the other

[gdal-dev] Questions on working with HYCOM data in netcdf format

2016-06-03 Thread Stephen Woodbridge
Hi All, I want to convert some HYCOM data into GTiff format and I'm a little lost on how to get started, whether I can do what I need just using a vrt or if I'm going to have to write code to process the data. To start with I have placed gdalinfo on one file here:

Re: [gdal-dev] Reading DNC nautical charts

2015-11-22 Thread Stephen Woodbridge
It looks like the DNC you are looking at are in Vector Product Format, so google "gdal vpf" and you will get some links that might help with that. -Steve On 11/22/2015 4:06 AM, Hans Rijsdijk wrote: Hi, Here is a link:

Re: [gdal-dev] Reading DNC nautical charts

2015-11-22 Thread Stephen Woodbridge
There are two version of Nautical charts from what I can see: 1. DNC from www.nga.mil that you pointed at in VPF format 2. NOAA ENC charts Here is a link explaining the difference: http://www.nauticalcharts.noaa.gov/mcd/learn_diffENC_DNC.html The ENC charts are available in S57 format and

Re: [gdal-dev] Adopt RFC48: Geographical networks support

2015-07-31 Thread Stephen Woodbridge
On 7/31/2015 2:51 PM, Dmitry Baryshnikov wrote: Hi everybody, The motion of RFC48 has been adopted with support from PSC members JukkaR, TamasS and EvenR. The code merged in trunk now (r29585). Let me know if issues arise. Also the RFC

[gdal-dev] How is NODATA interperted with a color table?

2015-06-24 Thread Stephen Woodbridge
I have a geotiff with one band of byte with a color table, but it does not have NODATA set in it. How would I set it NODATA to be to be entry 0 in the color table. $ gdalinfo A2015173174000.L0_LAC.L2_OC.tif Driver: GTiff/GeoTIFF Files: A2015173174000.L0_LAC.L2_OC.tif Size is 2433, 1727

Re: [gdal-dev] gdal_calc.py looses the color lookup table

2015-06-24 Thread Stephen Woodbridge
On 6/24/2015 8:59 AM, Even Rouault wrote: Le mercredi 24 juin 2015 03:49:15, Stephen Woodbridge a écrit : Hi All, I have a geotiff file like: [] ~/work/oceandata/test$ gdalinfo A2015173174000.L0_LAC.L2_OC.tif -noct Driver: GTiff/GeoTIFF Files: A2015173174000.L0_LAC.L2_OC.tif

Re: [gdal-dev] ogr2ogr: Stuck trying to append lots of gml files to postgis tables

2015-03-22 Thread Stephen Woodbridge
On 3/22/2015 3:32 PM, Even Rouault wrote: Le dimanche 22 mars 2015 20:18:12, Stephen Woodbridge a écrit : Hi all, I have a directory of gml files. I can load any of them into postgis without a problem. What I'm having a problem with is appending each file to the existing tables. I'm using

Re: [gdal-dev] ogr2ogr: Stuck trying to append lots of gml files to postgis tables

2015-03-22 Thread Stephen Woodbridge
, there's also the lco option of -lco OVERWRITE=YES Which also overwrites an existing layer, but it will retain any VIEWS in the postgresql database that are built using the layer you're appending or updating. Donovan On Sun, Mar 22, 2015 at 1:33 PM, Stephen Woodbridge wood...@swoodbridge.com

Re: [gdal-dev] a postgis table to another

2014-11-16 Thread Stephen Woodbridge
On 11/16/2014 8:01 AM, Ahmet Temiz wrote: Hello Is it possible to export (CREATE )a postgis table to another in different projection? regards There are a few ways that I can think of to do this: 1. add another column to the table with the project you want. Then export using that column,

[gdal-dev] Question on scaling between tif and vrt

2014-05-19 Thread Stephen Woodbridge
Hi all, I have a .tif file with -ot Int16 and an associated .vrt file that is applying a color palette and is type Byte. How does if decide on how to scale the pixel values? Does it compute the range and then just scale and offset that values to the new range? Below are the gdalinfo for

Re: [gdal-dev] Question on scaling between tif and vrt

2014-05-19 Thread Stephen Woodbridge
only the low or high order byte as the value? 3. something else? Thank you for your patience with all my questions. -Steve On 5/19/2014 5:48 PM, Even Rouault wrote: Le lundi 19 mai 2014 23:42:06, Stephen Woodbridge a écrit : Hi all, I have a .tif file with -ot Int16 and an associated .vrt

Re: [gdal-dev] How to map Int32 or Float32 GTiff to color map?

2014-05-15 Thread Stephen Woodbridge
On 5/15/2014 1:56 PM, Even Rouault wrote: Le jeudi 15 mai 2014 16:25:12, Stephen Woodbridge a écrit : Hi all, I'm trying to create a VRT file for a GTiff that is Int32 or Float32 to add a ColorTable but I seem to be missing a key piece of information. How are the pixel values mapped

Re: [gdal-dev] GSoC proposal looking for mentors and suggestions

2014-03-20 Thread Stephen Woodbridge
I believe that OSGeo is expecting every approved student to have a mentor and a co-mentor this year. I have been a mentor for the last 5+ years for pgRouting. We have had 2 students most of those years and two mentors with each mentor being the co-mentor for the other project. This has worked

[gdal-dev] Fwd: Re: Problem with results on two different versions of gdal_wrap

2014-03-09 Thread Stephen Woodbridge
Sorry meant to send this to the list. Original Message Subject: Re: [gdal-dev] Problem with results on two different versions of gdal_wrap Date: Sun, 09 Mar 2014 15:22:23 -0400 From: Stephen Woodbridge wood...@swoodbridge.com To: Even Rouault even.roua...@mines-paris.org

Re: [gdal-dev] Problem with results on two different versions of gdal_wrap

2014-03-09 Thread Stephen Woodbridge
On 3/8/2014 5:22 AM, Even Rouault wrote: This is the problem step: gdalwarp -rcs -ts 8800 6600 -s_srs EPSG:32662 -t_srs EPSG:4326 temp.tif target.tif gdalinfo -mm -stats target.tif is showing that the range of values in the image are dramatically different on the two servers! summary old:

[gdal-dev] Problem with results on two different versions of gdal_wrap

2014-03-07 Thread Stephen Woodbridge
Hi All, First off, sorry this is really long with all the gdalinfo pasted into it, but I thought is might be more useful to have real information to look at. I have been running down a problem caused by upgrading a server to a newer version of gdal. Here are the details: old server: GDAL

Re: [gdal-dev] GSoC 2014

2014-02-17 Thread Stephen Woodbridge
Mikhail, This is a very interesting idea. You might want to add specific support pgRouting which already supports building graphs, creating routing topology and solving various graph problem using postgis for geometry and tables for building and linking the topology. You might also want to

Re: [gdal-dev] gdalwarp EPSG:32662 problem

2014-02-08 Thread Stephen Woodbridge
On 2/8/2014 2:55 AM, Jean-Claude Repetto wrote: On 08/02/2014 07:11, Stephen Woodbridge wrote: I've run into a problem that I think is related to a change in the EPSG definition for EPSG:32662 The symptom is that the image after gdalwarp reprojection is about 20km south of where it should

Re: [gdal-dev] How to fix gdalwarp warning?

2014-02-08 Thread Stephen Woodbridge
On 2/8/2014 6:44 PM, Even Rouault wrote: Le vendredi 07 février 2014 16:07:46, Jukka Rahkonen a écrit : Stephen Woodbridge woodbri at swoodbridge.com writes: How does one do Converting the dataset prior to 24/32 bit is advised.? For example with http://www.gdal.org/pct2rgb.html

[gdal-dev] How to fix gdalwarp warning?

2014-02-07 Thread Stephen Woodbridge
How does one do Converting the dataset prior to 24/32 bit is advised.? Warning: Input file /maps/images.tmp/844219fb7ca86093d82b14376b56b875/modis-chlo ra-201402071850.tif has a color table, which will likely lead to bad results whe n using a resampling method other than nearest neighbour.

[gdal-dev] gdalwarp EPSG:32662 problem

2014-02-07 Thread Stephen Woodbridge
Hi, Sorry for cross posting. This seems like a proj4 issue but might be gdal related. I've run into a problem that I think is related to a change in the EPSG definition for EPSG:32662 The symptom is that the image after gdalwarp reprojection is about 20km south of where it should be.

Re: [gdal-dev] Large shapefile issues

2013-11-25 Thread Stephen Woodbridge
Unless something has changed, I have never been able to work with dbf file over 2GB using shapelib. -Steve W On 11/25/2013 5:52 PM, Even Rouault wrote: Le lundi 25 novembre 2013 11:42:23, CARMAN, Darren a écrit : Hi List I notice on the OGR formats page for ESRI Shapefile the following is

Re: [gdal-dev] Slow convertion from OSM to PG with -skipfailures

2013-05-27 Thread Stephen Woodbridge
Even, Would it make sense to create and load a stored procedure that you use to load the feature with, then you can trap exceptions and ignore them. I would think that this would be much faster. Well the stored procedure is run in a transaction bout you could do something like: begin;

Re: [gdal-dev] Slow convertion from OSM to PG with -skipfailures

2013-05-27 Thread Stephen Woodbridge
On 5/27/2013 3:06 PM, Even Rouault wrote: Le lundi 27 mai 2013 20:52:57, Stephen Woodbridge a écrit : Even, Would it make sense to create and load a stored procedure that you use to load the feature with, then you can trap exceptions and ignore them. I would think that this would be much

Re: [gdal-dev] PostGIS date field output

2013-02-28 Thread Stephen Woodbridge
On 2/28/2013 4:40 PM, Kralidis,Tom [Ontario] wrote: Hi: we are using 1.9.2 via MapServer WFS to serve out PostGIS data in GeoJSON, which works very well -- great feature! I was having issues with GeoJSON output on date fields when hour is 10. Digging deeper: The column in question is of

[gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
Hi Even, Frank pointed me to your blog entry on the GDAL virtual file system at: http://erouault.blogspot.com/2012/05/new-gdal-virtual-file-system-to-read.html This is most excellent. I am trying to get this to work using the /vsicurl/ option to connect to a php page that generates CSV

Re: [gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
Just a followup to this, I have downloaded gdal-1.9.2 and get the same results. So no joy there. -Steve On 2/9/2013 4:28 PM, Stephen Woodbridge wrote: Hi Even, Frank pointed me to your blog entry on the GDAL virtual file system at: http://erouault.blogspot.com/2012/05/new-gdal-virtual-file

Re: [gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
On 2/9/2013 6:03 PM, Jukka Rahkonen wrote: Stephen Woodbridge woodbri at swoodbridge.com writes: Just a followup to this, I have downloaded gdal-1.9.2 and get the same results. So no joy there. -Steve Perhaps vsicurl_streaming that is also mentioned in the same blog is something worth

Re: [gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
On 2/9/2013 6:41 PM, Even Rouault wrote: Le dimanche 10 février 2013 00:23:37, Stephen Woodbridge a écrit : On 2/9/2013 6:03 PM, Jukka Rahkonen wrote: Stephen Woodbridge woodbri at swoodbridge.com writes: Just a followup to this, I have downloaded gdal-1.9.2 and get the same results. So

Re: [gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
On 2/9/2013 6:41 PM, Even Rouault wrote: Le dimanche 10 février 2013 00:23:37, Stephen Woodbridge a écrit : On 2/9/2013 6:03 PM, Jukka Rahkonen wrote: Stephen Woodbridge woodbri at swoodbridge.com writes: Just a followup to this, I have downloaded gdal-1.9.2 and get the same results. So

Re: [gdal-dev] Virtual Filesystem Question(s)

2013-02-09 Thread Stephen Woodbridge
On 2/9/2013 8:24 PM, Stephen Woodbridge wrote: On 2/9/2013 6:41 PM, Even Rouault wrote: Le dimanche 10 février 2013 00:23:37, Stephen Woodbridge a écrit : On 2/9/2013 6:03 PM, Jukka Rahkonen wrote: Stephen Woodbridge woodbri at swoodbridge.com writes: Just a followup to this, I have

  1   2   >