[gdal-dev] Gdal 32 bit libraries -reg
Hi , Please let me know where to download Gdal libraries of 32 bit Thanks Nikhil Sai Parupalli ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Display of S57 data sets
Hi All, I could parse the S57 data set to my application written in writing it in Windows C# using VS2010 but I am stuck without a module to display the nautical charts. Nautical charts are formatted to the NOAA S57 or S60 formats of the original .BSB format. Does anyone know of a DLL or plug in or source code that I can use to display these chart formats. These particular formats are exclusive to the nautical world as such it is difficult to locate code (dll etc) that I can use. And I am yet to find any examples on how to build your own viewer. Thanks Nikhil Sai Parupalli ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Gdal compatibility with windows 8
Hi All. After all the system migration I could use GDAL in my windows 8 system and as of now no problem has encountered Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: gdal-dev-boun...@lists.osgeo.org [gdal-dev-boun...@lists.osgeo.org] on behalf of Joaquim Luis [jl...@ualg.pt] Sent: Thursday, February 07, 2013 7:45 PM To: Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] Gdal compatibility with windows 8 Nikhil, I don't have the python tests running and my experience on a Win8 machine is still limited (it's not my main working computer) but so far I haven't encounter any problem with either VC10 and VC11 GDAL builds (which I built too). Joaquim Hi All, I have used gdal latest version for my application developed on windows 7.. As our organisation is migrating to windows 8 , I want to check whether gdal is supported in windows 8 Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Gdal compatibility with windows 8
Hi All, I have used gdal latest version for my application developed on windows 7.. As our organisation is migrating to windows 8 , I want to check whether gdal is supported in windows 8 Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] FIDN-issue
Hi , In my application I am fetching S57 data set where it is fetching FIDN with negative value which is wrong because FIDN should be in Range: 1 to 2 power 32-1 Is there any instance that any one got the same issue. If the same data source is loaded in Caris Easy view it is showing some other legal value. Please load the sample data with any application in open source it will show the same as mentioned above Thanks and Regards Nikhil Sai Parupalli ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] (no subject)
Hi All, With the some useful tips and code from gdal-dev@lists.osgeo.org I could successfully complete a huge task to validate S58 checks . Next up is developing a software where it creates a S57 data set . Please give me any information pertaining to S57 data set creation Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] map selection -reg
Hi All,
I am working on an application where I have used Dotspatail dll for display of
map.
How do we get selected layer information.
I have tried implementing the same with following code but all is not well
DotSpatial.Controls.IMapLayer laySelected = null;
DotSpatial.Symbology.FeatureLayer fea = null;
//string layerSelected = "";
foreach (DotSpatial.Controls.IMapLayer lay in
uxMap.GetFeatureLayers())
{
fea = (FeatureLayer)lay;
ISelection sel = fea.Selection;
// F = sel.ToFeatureSet();
List ls1 = new List();
ls1 = sel.ToFeatureList();
if (ls1.Count > 0)
{
laySelected = lay;
break;
}
}
Thanks and Regards
Nikhil Sai Parupalli
Note: Do not print this email until and unless it is really required. Save
paper , stay Green
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] How to read catalog file
Hi All, I have developed an application on VB.Net with GDAL as interface to read S57 dataset. Everything is fine , here I have to read some information form S57 catalog file . I have tried using stream reader in VB.net, but unable to fetch the values correctly. Can any one help me out in reading S57 catalog file using GDAL and vb.net . Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] ordering of layers
Hi ,
I am trying to load a S57 dataset to my map control using GDAL.
I have diff types of layers i.e
Polygon, Point, Line and multipoint
Here I want to show polygon in the base and points and lines on top of it my
code is as follows
if (layerType == "POLYGON")
{
FeatureSet fs = new FeatureSet(dt, 1, true,
DotSpatial.Topology.FeatureType.Polygon);
fs.DataTable = dt.Clone();
var c = new
DotSpatial.Topology.Coordinate[arrCoordinatesX.Count - 1];
for (int k = 0; k < arrCoordinatesX.Count - 1;
k++)
{
c[k] = new
DotSpatial.Topology.Coordinate(Convert.ToDouble(arrCoordinatesX[k]),
Convert.ToDouble(arrCoordinatesY[k]));
}
DotSpatial.Topology.Polygon LS = new
DotSpatial.Topology.Polygon(c);
DotSpatial.Data.IFeature currentFeature =
fs.AddFeature(LS);
uxMap.Layers.Add(fs);
}
else if (layerType == "LINESTRING")
{
FeatureSet fs = new FeatureSet(dt, 1, true,
DotSpatial.Topology.FeatureType.Line);
fs.DataTable = dt.Clone();
List
coordinates = new List();
for (int k = 0; k <= arrCoordinatesX.Count - 1;
k++)
{
coordinates.Add(new
DotSpatial.Topology.Coordinate(Convert.ToDouble(arrCoordinatesX[k]),
Convert.ToDouble(arrCoordinatesY[k])));
}
DotSpatial.Topology.LineString LS = new
DotSpatial.Topology.LineString(coordinates);
DotSpatial.Data.IFeature currentFeature =
fs.AddFeature(LS);
uxMap.Layers.Add(fs);
foreach (DataRow row in dt.Rows)
{
}
}
else if (layerType == "POINT")
{
FeatureSet fs = new FeatureSet(dt, 1, true,
DotSpatial.Topology.FeatureType.Point);
fs.DataTable = dt.Clone();
foreach (DataRow row in dt.Rows)
{
for (int k = 0; k <= arrCoordinatesX.Count
- 1; k++)
{
DotSpatial.Topology.Coordinate c = new
DotSpatial.Topology.Coordinate(Convert.ToDouble(arrCoordinatesX[k]),
Convert.ToDouble(arrCoordinatesY[k]));
DotSpatial.Topology.Point p = new
DotSpatial.Topology.Point(c);
DotSpatial.Data.IFeature currentFeature
= fs.AddFeature(p);
//PointCheck(ds);
}
uxMap.Layers.Add(fs);
//uxMap.Layers.OrderByDescending(fs, 1);
}
}
else if (layerType == "MULTIPOINT")
{
FeatureSet fs = new FeatureSet(dt, 1, true,
DotSpatial.Topology.FeatureType.Point);
fs.DataTable = dt.Clone();
foreach (DataRow row in dt.Rows)
{
DotSpatial.Topology.Point p = new
DotSpatial.Topology.Point();
for (int k = 0; k <= arrCoordinatesX.Count
- 1; k++)
{
DotSpatial.Topology.Coordinate c = new
DotSpatial.Topology.Coordinate(Convert.ToDouble(arrCoordinatesX[k]),
Convert.ToDouble(arrCoordinatesY[k]));
p = new DotSpatial.Topology.Point(c);
DotSpatial.Data.IFeature currentFeature
= fs.AddFeature(p);
}
DotSpatial.Topology.MultiPoint mp = new
DotSpatial.Topology.MultiPoint(p);
uxMap.Layers.Add(fs);
}
}
}
dt.Rows.Clear();
}
}
So how to control this layer order so that i can see all the layers in a single
map window
Thanks and Regards
Nikhil Sai Parupalli
Note: Do not print thi
[gdal-dev] dotspatial map control navigation
Hi I am developing an application which has dotspatial map control do display S57 dataset. Here it is displaying fine but I have to write a code where user should be navigated to particular coordinates when clicking and RCID(it has coordinate information). If you have worked on similar type of thing do share me some idea. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Map control in dotspatial
Hi , Here in my application I have a dotspatial map control , when user clicks particular place on map how to get information of that particular point. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] 2 geometries share function
Hey Loskot This code deals with lot of sql stmts. but here in our development we never worked with databases. We are using gdal libraries , is there any function in gdal or any other libraries that serves the purpose Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: Mateusz Loskot [mate...@loskot.net] Sent: Monday, October 15, 2012 3:44 PM To: Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] 2 geometries share function On 15 October 2012 10:46, Nikhil Sai Parupalli wrote: > Hi all, > > What is the function to find two geometries shares common boundary. > Is that geometry.touches(geometry) PostGIS docs provide a fantastic visual aid: http://postgis.refractions.net/documentation/manual-2.0/ST_Touches.html Best regards, -- Mateusz Loskot, http://mateusz.loskot.net ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] 2 geometries share function
Hi all, What is the function to find two geometries shares common boundary. Is that geometry.touches(geometry) Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Data Set Parameter [DSPM]
Hi, Can any one help me out how to get Data Set Parameter [DSPM] field for S57 dataset into application using GDAL. I need to check below things for my application Check that HDAT = 2 [WGS 84]. Check that DUNI = 1 [metres]. Check that PUNI = 1 [metres]. Check that COUN = 1 [latitude/longitude]. Thanks and Regards Nikhil Sai Parupalli ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Open Source GIS Raster Viewer
Hey Kishore, We can use DotSpatial library that provides Map control wjer we can display map and ot has standard ZoomIn,ZoomOut ,PAN etc... I am using the same in my application which is working fine. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: gdal-dev-boun...@lists.osgeo.org [gdal-dev-boun...@lists.osgeo.org] on behalf of kishore Borra [kishorebo...@gmail.com] Sent: Friday, October 12, 2012 6:22 PM To: gdal-dev@lists.osgeo.org Subject: [gdal-dev] Open Source GIS Raster Viewer Hi All, we are implementing support of elevation processing of various raster images in a standalone desktop Viewer using opensource. we are using GDAL library to support performing various operations like hillshade, slope, aspect map & color coded relief. I would like to know your suggestions on what opensource control can be used as a map control so that I can load various type of processed raster images into it. I need to provide the user with the coordinates readout information and standard zoom In, zoom out & zoom By area operations. I should also be able to read the elevation info from the displayed raster and display the same in the status bar. Please suggest. Regards, Kishore. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Data Set Parameter -reg
Hi, Can any one help me out how to get Data Set Parameter [DSPM] fieldfor S57dataset into application using GDAL. Do i need to change any system variables. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] S57 update (was: (no subject))
Hey , This was the check i was referencing to For record updates for feature/vector updates, check that if it is DELETE: the record does not contain further fields, or MODIFY/INSERT: the record contains more information about the update. Could you understand something about this... Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: s duclos [sylvain_duc...@yahoo.com] Sent: Friday, October 12, 2012 12:19 AM To: Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] S57 update (was: (no subject)) Hi Nikhil, Frank as put up a page for this (sorry Frank for highjacking the thread) http://www.gdal.org/ogr/drv_s57.html Also IHO have the specs for updating ENC specifically. (I don't have the link right now) rgds, Sylvain. - Original Message - > From: Nikhil Sai Parupalli > To: s duclos > Cc: "gdal-dev@lists.osgeo.org" > Sent: Thursday, October 11, 2012 2:26:54 PM > Subject: RE: [gdal-dev] S57 update (was: (no subject)) > > Hey > > HOw do we get to know in data set i.e any attribute information varies for > update... > > Thanks and Regards > Nikhil Sai Parupalli > > > > Note: Do not print this email until and unless it is really required. Save > paper > , stay Green > > > From: s duclos [sylvain_duc...@yahoo.com] > Sent: Thursday, October 11, 2012 8:55 PM > To: Nikhil Sai Parupalli; gdal-dev@lists.osgeo.org > Subject: Re: [gdal-dev] S57 update (was: (no subject)) > > Hi Nikhil, > > In the US1GC09M.ZIP you send you can see update as > > US1GC09M.000 <-- base cell > > > US1GC09M.001 - .008 <-- update to base cell > > > > rgds, > > Sylvain. > > > >> >> From: Nikhil Sai Parupalli >> To: "gdal-dev@lists.osgeo.org" >> Sent: Thursday, October 11, 2012 2:28:40 AM >> Subject: [gdal-dev] (no subject) >> >> >> >> Hi All, >> >> >> How do we get to know which S57 dataset has update sin it.. >> >> >> >> >> Thanks and Regards >> Nikhil Sai Parupalli >> >> >> >> Note: Do not print this email until and unless it is really required. Save > paper , stay Green >> >> ___ >> gdal-dev mailing list >> gdal-dev@lists.osgeo.org >> http://lists.osgeo.org/mailman/listinfo/gdal-dev >> >> > ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] S57 update (was: (no subject))
Hey HOw do we get to know in data set i.e any attribute information varies for update... Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: s duclos [sylvain_duc...@yahoo.com] Sent: Thursday, October 11, 2012 8:55 PM To: Nikhil Sai Parupalli; gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] S57 update (was: (no subject)) Hi Nikhil, In the US1GC09M.ZIP you send you can see update as US1GC09M.000 <-- base cell US1GC09M.001 - .008 <-- update to base cell rgds, Sylvain. > > From: Nikhil Sai Parupalli >To: "gdal-dev@lists.osgeo.org" >Sent: Thursday, October 11, 2012 2:28:40 AM >Subject: [gdal-dev] (no subject) > > > >Hi All, > > >How do we get to know which S57 dataset has update sin it.. > > > > >Thanks and Regards >Nikhil Sai Parupalli > > > >Note: Do not print this email until and unless it is really required. Save >paper , stay Green > >___ >gdal-dev mailing list >gdal-dev@lists.osgeo.org >http://lists.osgeo.org/mailman/listinfo/gdal-dev > > ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] (no subject)
Hi All, How do we get to know which S57 dataset has update sin it.. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] S57 dataset boundary issue
Hey Sylvain, Thanks for your inputs, we have successfully implemented the same in our project . I have one more question i.e If I have a polygon how to find out edges of it programmatic do we have any function to get all the edges of polygon(area) Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: Nikhil Sai Parupalli Sent: Wednesday, October 10, 2012 9:07 PM To: s duclos Subject: RE: [gdal-dev] S57 dataset boundary issue Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: s duclos [sylvain_duc...@yahoo.com] Sent: Wednesday, October 10, 2012 8:43 PM To: Nikhil Sai Parupalli; Frank Warmerdam Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] S57 dataset boundary issue Hi Nikhil, > Based on the code which you have sent we could find out CW and CCW, but not > outer/inner The outer ring of OGRGeometryH poly is the ring 0 (zero). In C code the OGR call: guintnRingCount = OGR_G_GetGeometryCount(hGeom); will give you the number of rings in a poly. Say there is 2 rings for this hGeom. Then the outer ring is OGRGeometryH *hGeomRef_outter = OGR_G_GetGeometryRef(hGeom, 0); and the inner ring is OGRGeometryH *hGeomRef_inner = OGR_G_GetGeometryRef(hGeom, 1); > Please have a look at that data set once and let me know. There is nothing in attachment rgds, Sylvain. - Original Message - > From: Nikhil Sai Parupalli > To: s duclos ; Frank Warmerdam > Cc: "gdal-dev@lists.osgeo.org" > Sent: Wednesday, October 10, 2012 10:40:27 AM > Subject: RE: [gdal-dev] S57 dataset boundary issue > > Hey Sylvain, > > Thanks for trying the data which I have sent. > > Could you be helping me how this works like outer boundary and inner boundary. > > > > Based on the code which you have sent we could find out CW and CCW, but not > outer/inner > Please have a look at that data set once and let me know. > Or if you have any sample data with simple representation of outer and inner > boundaries please do share with me. > > > > > Thanks and Regards > Nikhil Sai Parupalli > > > > Note: Do not print this email until and unless it is really required. Save > paper > , stay Green > > ____________ > From: s duclos [sylvain_duc...@yahoo.com] > Sent: Wednesday, October 10, 2012 7:59 PM > To: Nikhil Sai Parupalli; Frank Warmerdam > Cc: gdal-dev@lists.osgeo.org > Subject: Re: [gdal-dev] S57 dataset boundary issue > > Hi Nikhil, > > I loaded US1GC09M//US1GC09M.000 .. it's huge :) > > It's the Gulf of Mexico .. now where should I look ? > > > > About the inner / outer ring. As I said before, I just assume that OGR > send me the outer ring first. And, up to now, my code is working fine. > > > rgds, > > Sylvain. > > > > > - Original Message - >> From: Nikhil Sai Parupalli >> To: s duclos ; Frank Warmerdam > >> Cc: "gdal-dev@lists.osgeo.org" >> Sent: Wednesday, October 10, 2012 9:40:54 AM >> Subject: RE: [gdal-dev] S57 dataset boundary issue >> >> Could you try on this data set please >> >> Thanks and Regards >> Nikhil Sai Parupalli >> >> >> >> Note: Do not print this email until and unless it is really required. Save > paper >> , stay Green >> >> >> From: s duclos [sylvain_duc...@yahoo.com] >> Sent: Wednesday, October 10, 2012 7:02 PM >> To: Nikhil Sai Parupalli; Frank Warmerdam >> Cc: gdal-dev@lists.osgeo.org >> Subject: Re: [gdal-dev] S57 dataset boundary issue >> >> Nikhil, >> >> I loaded your ENC. >> >> I see a number of depth area (DEPARE) and anchor zone area (ACHARE). >> >> But no inner ring! >> >> Is there an inner rings somewhere I should see? >> >> >> >> rgds, >> >> Sylvain. >> >> >> >> - Original Message - >>> From: Nikhil Sai Parupalli > >>> To: s duclos ; Frank Warmerdam >> >>> Cc: "gdal-dev@lists.osgeo.org" > >>> Sent: Wednesday, October 10, 2012 8:00:41 AM >>> Subject: RE: [gdal-dev] S57 dataset boundary issue >>> >>> Duclos, >>> >>> Please see this sample data attached. >>> So as per the code given by you we have tried and got result i.e able > to >> find >>> out area >
Re: [gdal-dev] S57 dataset boundary issue
Hey Sylvain,
Thanks for trying the data which I have sent.
Could you be helping me how this works like outer boundary and inner boundary.
Based on the code which you have sent we could find out CW and CCW, but not
outer/inner
Please have a look at that data set once and let me know.
Or if you have any sample data with simple representation of outer and inner
boundaries please do share with me.
Thanks and Regards
Nikhil Sai Parupalli
Note: Do not print this email until and unless it is really required. Save
paper , stay Green
From: s duclos [sylvain_duc...@yahoo.com]
Sent: Wednesday, October 10, 2012 7:59 PM
To: Nikhil Sai Parupalli; Frank Warmerdam
Cc: gdal-dev@lists.osgeo.org
Subject: Re: [gdal-dev] S57 dataset boundary issue
Hi Nikhil,
I loaded US1GC09M//US1GC09M.000 .. it's huge :)
It's the Gulf of Mexico .. now where should I look ?
About the inner / outer ring. As I said before, I just assume that OGR
send me the outer ring first. And, up to now, my code is working fine.
rgds,
Sylvain.
- Original Message -
> From: Nikhil Sai Parupalli
> To: s duclos ; Frank Warmerdam
> Cc: "gdal-dev@lists.osgeo.org"
> Sent: Wednesday, October 10, 2012 9:40:54 AM
> Subject: RE: [gdal-dev] S57 dataset boundary issue
>
> Could you try on this data set please
>
> Thanks and Regards
> Nikhil Sai Parupalli
>
>
>
> Note: Do not print this email until and unless it is really required. Save
> paper
> , stay Green
>
>
> From: s duclos [sylvain_duc...@yahoo.com]
> Sent: Wednesday, October 10, 2012 7:02 PM
> To: Nikhil Sai Parupalli; Frank Warmerdam
> Cc: gdal-dev@lists.osgeo.org
> Subject: Re: [gdal-dev] S57 dataset boundary issue
>
> Nikhil,
>
> I loaded your ENC.
>
> I see a number of depth area (DEPARE) and anchor zone area (ACHARE).
>
> But no inner ring!
>
> Is there an inner rings somewhere I should see?
>
>
>
> rgds,
>
> Sylvain.
>
>
>
> - Original Message -
>> From: Nikhil Sai Parupalli
>> To: s duclos ; Frank Warmerdam
>
>> Cc: "gdal-dev@lists.osgeo.org"
>> Sent: Wednesday, October 10, 2012 8:00:41 AM
>> Subject: RE: [gdal-dev] S57 dataset boundary issue
>>
>> Duclos,
>>
>> Please see this sample data attached.
>> So as per the code given by you we have tried and got result i.e able to
> find
>> out area
>>
>> some times it shows area<0 and some times area>0 hence we get CW and
> CCW.
>> But how do we classify in layer which is outer boundary and which is inner
>> boundary
>>
>>
>>
>>
>> Geometry testgeo = geomLAKARE.GetBoundary();
>> double area = 0;
>> for (int l = 0; l <
>> testgeo.GetPointCount() - 1; l++)
>> {
>> double x1 = testgeo.GetX(l);
>> double y1 = testgeo.GetY(l);
>> double x2 = testgeo.GetX(l + 1);
>> double y2 = testgeo.GetY(l + 1);
>> area += (x1 * y2) - (x2 * y1);
>> }
>>
>> Thanks and Regards
>> Nikhil Sai Parupalli
>>
>>
>>
>> Note: Do not print this email until and unless it is really required. Save
> paper
>> , stay Green
>>
>>
>> From: s duclos [sylvain_duc...@yahoo.com]
>> Sent: Tuesday, October 09, 2012 11:58 PM
>> To: Frank Warmerdam; Nikhil Sai Parupalli
>> Cc: gdal-dev@lists.osgeo.org
>> Subject: Re: [gdal-dev] S57 dataset boundary issue
>>
>> Frank & Nikhil,
>>
>>
>> Frank asked:
>>
>>> are, I'm not sure what I or others can do about
>>> the performance cost of checking.
>>
>>
>> If I understand the problem, then it cost nothing!
>>
>> You have to check the winding while loading a OGRGeometryH
>> like so:
>> double area = 0;
>> for (unsigned int i=0; i> double x1 = OGR_G_GetX(hRing, i);
>> double y1 = OGR_G_GetY(hRing, i);
>> double x2 = OGR_G_GetX(hRing, i+1);
>> double y2 = OGR_G_GetY(hRing, i+1);
>> area += (x1*y2) - (x2*y1);
>> }
>>
>>
>> If the area < 0 then it
[gdal-dev] S57 dataset boundary issue
HI All, , Thanks for the reply in below data NAME_RCNM (IntegerList) = (7:130,130,130,130,130,130,130) NAME_RCID (IntegerList) = (7:691,391,690,52,1336,1503,1512) ORNT (IntegerList) = (7:2,2,2,2,1,1,1) USAG (IntegerList) = (7:1,1,1,1,1,1,1) MASK (IntegerList) = (7:255,255,255,255,255,255,255) RCID 691 ORNT has 2 and USAG is 1 based on this data we could identify that RCID 691 is of exterior boundary with reverse direction. So if I have 1000 RCID in NAME_RCID which is vise versa for the above condition I need to display that RCIDS list as errors. So it makes 1000 conditional checks hence reducing performance. Please let us know further. Nikhil, I don't know why you don't see it. I see it fine on the "Edge" layer which contains the linestring vector objects. eg. OGRFeature(Edge):18 RCNM (Integer) = 130 RCID (Integer) = 19 RVER (Integer) = 1 RUIN (Integer) = 1 NAME_RCNM_0 (Integer) = 120 NAME_RCID_0 (Integer) = 21 ORNT_0 (Integer) = 255 USAG_0 (Integer) = 255 TOPI_0 (Integer) = 1 MASK_0 (Integer) = 255 NAME_RCNM_1 (Integer) = 120 NAME_RCID_1 (Integer) = 20 ORNT_1 (Integer) = 255 USAG_1 (Integer) = 255 TOPI_1 (Integer) = 2 MASK_1 (Integer) = 255 LINESTRING (60.775 -32.6249972) > I have one more question: > How do we find area outer boundaries are encoded clockwise based > on Orientation or USAG. > here when I retrieve S57 dataset I get multiple ORNT values for a single > object. I am not aware of any definition of outer ring winding direction. I think the ORNT field values you see on polygons are indicating which direction the edges were used in when composing the ring(s) of the polygon. NAME_RCNM (IntegerList) = (7:130,130,130,130,130,130,130) NAME_RCID (IntegerList) = (7:691,391,690,52,1336,1503,1512) ORNT (IntegerList) = (7:2,2,2,2,1,1,1) USAG (IntegerList) = (7:1,1,1,1,1,1,1) MASK (IntegerList) = (7:255,255,255,255,255,255,255) So in the above feature, RCID 691 is used with orientation "2", usage "1" and mask "255". Best regards, Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Topological pointer
Frank, Thanks for the reply in below data NAME_RCNM (IntegerList) = (7:130,130,130,130,130,130,130) NAME_RCID (IntegerList) = (7:691,391,690,52,1336,1503,1512) ORNT (IntegerList) = (7:2,2,2,2,1,1,1) USAG (IntegerList) = (7:1,1,1,1,1,1,1) MASK (IntegerList) = (7:255,255,255,255,255,255,255) RCID 691 ORNT has 2 and USAG is 1 based on this data we could identify that RCID 691 is of exterior boundary with reverse direction. So if I have 1000 RCID in NAME_RCID which is vise versa for the above condition I need to display that RCIDS list as errors. So it makes 1000 conditional checks hence reducing performance. Please let us know further. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: Srikanth Rao Sidduri Sent: Monday, October 08, 2012 11:42 AM To: Frank Warmerdam; Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: RE: [gdal-dev] Topological pointer Frank, Thanks for the reply in below data NAME_RCNM (IntegerList) = (7:130,130,130,130,130,130,130) NAME_RCID (IntegerList) = (7:691,391,690,52,1336,1503,1512) ORNT (IntegerList) = (7:2,2,2,2,1,1,1) USAG (IntegerList) = (7:1,1,1,1,1,1,1) MASK (IntegerList) = (7:255,255,255,255,255,255,255) RCID 691 ORNT has 2 and USAG is 1 based on this data we could identify that RCID 691 is of exterior boundary with reverse direction. So if I have 1000 RCID in NAME_RCID which is vise versa for the above condition I need to display that RCIDS list as errors. So it makes 1000 conditional checks hence reducing performance. Please let us know further.. Thanks & Regards Srikanth Rao Sidduri, From: gdal-dev-boun...@lists.osgeo.org [gdal-dev-boun...@lists.osgeo.org] on behalf of Frank Warmerdam [warmer...@pobox.com] Sent: Monday, October 08, 2012 10:49 AM To: Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] Topological pointer On 12-10-07 10:02 PM, Nikhil Sai Parupalli wrote: > Frank, > > As said I have my system variable as > > "RETURN_PRIMITIVES=ON,RETURN_LINKAGES=ON,LNAM_REFS=ON, > UPDATES=APPLY,SPLIT_MULTIPOINT=ON,ADD_SOUNDG_DEPTH=ON" > > but still I couldn't see TOPI field retrieved in my dataset.. > Will you elaborate and let me know. Nikhil, I don't know why you don't see it. I see it fine on the "Edge" layer which contains the linestring vector objects. eg. OGRFeature(Edge):18 RCNM (Integer) = 130 RCID (Integer) = 19 RVER (Integer) = 1 RUIN (Integer) = 1 NAME_RCNM_0 (Integer) = 120 NAME_RCID_0 (Integer) = 21 ORNT_0 (Integer) = 255 USAG_0 (Integer) = 255 TOPI_0 (Integer) = 1 MASK_0 (Integer) = 255 NAME_RCNM_1 (Integer) = 120 NAME_RCID_1 (Integer) = 20 ORNT_1 (Integer) = 255 USAG_1 (Integer) = 255 TOPI_1 (Integer) = 2 MASK_1 (Integer) = 255 LINESTRING (60.775 -32.6249972) > I have one more question: > How do we find area outer boundaries are encoded clockwise based > on Orientation or USAG. > here when I retrieve S57 dataset I get multiple ORNT values for a single > object. I am not aware of any definition of outer ring winding direction. I think the ORNT field values you see on polygons are indicating which direction the edges were used in when composing the ring(s) of the polygon. NAME_RCNM (IntegerList) = (7:130,130,130,130,130,130,130) NAME_RCID (IntegerList) = (7:691,391,690,52,1336,1503,1512) ORNT (IntegerList) = (7:2,2,2,2,1,1,1) USAG (IntegerList) = (7:1,1,1,1,1,1,1) MASK (IntegerList) = (7:255,255,255,255,255,255,255) So in the above feature, RCID 691 is used with orientation "2", usage "1" and mask "255". Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://home.gdal.org/warmerda and watch the world go round - Rush| Geospatial Software Developer ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Topological pointer
Frank, As said I have my system variable as "RETURN_PRIMITIVES=ON,RETURN_LINKAGES=ON,LNAM_REFS=ON, UPDATES=APPLY,SPLIT_MULTIPOINT=ON,ADD_SOUNDG_DEPTH=ON" but still I couldn't see TOPI field retrieved in my dataset.. Will you elaborate and let me know. I have one more question: How do we find area outer boundaries are encoded clockwise based on Orientation or USAG. here when I retrieve S57 dataset I get multiple ORNT values for a single object. Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green From: fwarmer...@gmail.com [fwarmer...@gmail.com] on behalf of Frank Warmerdam [warmer...@pobox.com] Sent: Saturday, October 06, 2012 2:32 AM To: Nikhil Sai Parupalli Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] Topological pointer Nikhil, >From the code it appears the TOPI field will be represented on the vector >primitives when you enable treating vector primitives as features. To enable >this you need to set the RETURN_PRIMITIVES=ON configuration option. See the >S57 driver page for some background. http://www.gdal.org/ogr/drv_s57.html It looks like the attributes will be named TOPI_0 and TOPI_1 and will be integer fields containing the start and end point id I gather. The S57 topology support is a pretty fuzzy in my head. There is also other topological information (I think) accessable on the vector primitive features. Some likely as hard to use list attributes. Best regards, Frank On Fri, Oct 5, 2012 at 11:06 AM, Nikhil Sai Parupalli mailto:nikhil.parupa...@iictechnologies.com>> wrote: Hi All, We are working on project which renders S57 dataset , I used GDAL for rendering GIS data and DOTSPATIAL for displaying map. Here I could read dataset , but can any one know how do we get Topological Indicator value Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> http://lists.osgeo.org/mailman/listinfo/gdal-dev -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com<mailto:warmer...@pobox.com> light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| Geospatial Software Developer ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Topological pointer
Hi All, We are working on project which renders S57 dataset , I used GDAL for rendering GIS data and DOTSPATIAL for displaying map. Here I could read dataset , but can any one know how do we get Topological Indicator value Thanks and Regards Nikhil Sai Parupalli Note: Do not print this email until and unless it is really required. Save paper , stay Green ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
