SQL Server Spatial Tools

Closed Issue: using GeographyCollectionAggregate on set of points that exceed hemispehere limit raise Exception [3737]

bsdz — Thu, 16 May 2013 10:24:36 GMT

When using GeographyCollectionAggregate or GeographyUnionAggregate on set of points that exceed the hemisphere limit an Exception is raised. What would be better is if NULL is returned as this would match the behaviour of STUnion and would make the aggregate classes more usable.

For example: -

declare @x geography = 'MULTIPOINT ((0 1), (90 1), (180 1), (270 1))';
declare @y geography = 'POINT (45 1)';

select @x.STUnion(@y)

This returns NULL because by unioning @y to @x the new geography exceeds the hemisphere limit.

However if using, say GeographyCollectionAggregate, like so: -

declare @t table (
name varchar(max),
position geography
);

insert into @t (name, position) values ('object 1', 'POINT (0 1)' );
insert into @t (name, position) values ('object 1', 'POINT (90 1)' );
insert into @t (name, position) values ('object 1', 'POINT (180 1)' );
insert into @t (name, position) values ('object 1', 'POINT (270 1)' );
/* including this point causes hemisphere limit Exception */
insert into @t (name, position) values ('object 1', 'POINT (45 1)' );

select name, dbo.GeographyCollectionAggregate(position) 'locus'
from @t
group by name

The inclusion of the POINT(45 1) forces the 'locus' value to exceed the hemisphere limit but instead of returning NULL a GLArgumentException is raised.

Created Issue: Implement a method to export each geography in a resultset as a KML placemark fragment [10777]

Bigjimslade — Tue, 29 Jan 2013 05:09:12 GMT

Often when working with more than one row of data in a sql resultset, you want to return a fragment of KML, not a entire well formed kml document.

I would like to suggest that the spatial tools library implement an additional export method that just return a element as xml and optionally takes a parameter of xml to include as an node.

Patch Uploaded: #13750

Bigjimslade — Tue, 29 Jan 2013 04:47:28 GMT

Bigjimslade has uploaded a patch.

Description:
This patch corrects a missing namespace in the KMLFunctions.cs that prevents registration of the KML functions in SQL Server. I also changed the methods to return proper xml as suggested
in #8717

Finally, the register and unregister scripts have been updated to include the kml functions.

New Post: Project Status

Bigjimslade — Mon, 28 Jan 2013 18:57:02 GMT

ok great I was trying to use the KML Features and it looks like there is a minoor issue with using it I think i have a fix would it be possible to get added as a developer?

Thanks.

New Post: Project Status

edwink — Mon, 28 Jan 2013 18:54:05 GMT

We’ve just finished rolling up most of the features in this project into SQL Server 2012 and are now doing the same for the Windows Azure SQL Database (formerly SQL Azure). This project will remain active on CodePlex while we determine what features to add next after we finish the SQL and Azure synchronizations.

From: Bigjimslade [email removed]
Sent: Monday, January 28, 2013 10:44 AM
To: Ed Katibah
Subject: Project Status [sqlspatialtools:431141]

From: Bigjimslade

What is the status of this project? is anyone still doing active development?

New Post: Project Status

Bigjimslade — Mon, 28 Jan 2013 18:44:01 GMT

What is the status of this project? is anyone still doing active development?

New Post: Showing two spatial columns in the spatial results together

Viking42 — Fri, 30 Nov 2012 07:50:22 GMT

works like a charm!

Thank you very much

New Post: Showing two spatial columns in the spatial results together

edwink — Wed, 28 Nov 2012 18:16:02 GMT

You could try something like this:

select top 100 geography::Point(ddLatitude, ddLongitude, 4326) as point, geom into #foo

from Area, Position

where geography::Point(ddLatitude, ddLongitude, 4326).STIntersects(Area.geom) = 0

select point.STBuffer(1000) from #foo

union all

select geom from #foo

The general approach would be to populate a temporary table with the output from the select (into #foo). Then you can logically merge the two spatial columns in #foo using the union all command.

Note the simplified point constructor, Point(). When you are using the geography type and the Point() method in argument list is LATITUDE, LONGITUDE, SRID.

Please let us know if this answers your question.

Thanks,

-Ed

From: Viking42 [email removed]
Sent: Wednesday, November 28, 2012 8:44 AM
To: Ed Katibah
Subject: Showing two spatial columns in the spatial results together [sqlspatialtools:404916]

From: Viking42

Hi,

I am new to this and I have some geo points that does not intersects with my geo structures and I want to plot these and see where these points are in order to debug the solution.

I run the following query in MS SQL management studio:

select

top 100

geography::STPointFromText('POINT(' + CAST(ddLongitude AS VARCHAR) + ' ' + CAST(ddLatitude AS VARCHAR) + ')', 4326)

,geom

from Area, Position

where

geography

::STPointFromText('POINT(' + CAST(ddLongitude AS VARCHAR) + ' ' + CAST(ddLatitude AS VARCHAR) + ')', 4326).STIntersects(Area.geom) = 0

I can switch between the two columns in the Spatial Results window but I want to see the points and my shapes together. Is there a way to do this?

Thank you for your replies

New Post: Showing two spatial columns in the spatial results together

Viking42 — Wed, 28 Nov 2012 16:43:33 GMT

Hi,

I am new to this and I have some geo points that does not intersects with my geo structures and I want to plot these and see where these points are in order to debug the solution.

I run the following query in MS SQL management studio:

select

top 100

geography::STPointFromText('POINT(' + CAST(ddLongitude AS VARCHAR) + ' ' + CAST(ddLatitude AS VARCHAR) + ')', 4326)

,geom

from Area, Position

where

geography

::STPointFromText('POINT(' + CAST(ddLongitude AS VARCHAR) + ' ' + CAST(ddLatitude AS VARCHAR) + ')', 4326).STIntersects(Area.geom) = 0

I can switch between the two columns in the Spatial Results window but I want to see the points and my shapes together. Is there a way to do this?

Thank you for your replies

Source code checked in, #80132

Project Collection Service Accounts — Mon, 01 Oct 2012 21:49:39 GMT

Upgrade: New Version of LabDefaultTemplate.xaml. To upgrade your build definitions, please visit the following link: http://go.microsoft.com/fwlink/?LinkId=254563

Source code checked in, #80131

Project Collection Service Accounts — Mon, 01 Oct 2012 21:43:10 GMT

Checked in by server upgrade

Patch Uploaded: #11602

kaveh096 — Wed, 29 Feb 2012 03:36:32 GMT

kaveh096 has uploaded a patch.

Description:
Continuing from the discussion here (http://sqlspatialtools.codeplex.com/discussions/69864) I noticed that the TranverseMercator implementation results in slightly different numbers than other GIS tools like ESRI ArcGIS. So I decided to create another class for my own needs (TranverseMercatorMeter). I was wondering if you guys are interested in this so I decided to submit a patch. Let me know how else I can help in this regard.

New Post: Project to UTM - long0 values

kaveh096 — Wed, 29 Feb 2012 01:59:37 GMT

Don't know if this issue is still valid... I'm simply posting for the sake of future reference. I had some WSG 1984 geographies which needed to be projected into UTM on run-time using TranverseMercator. So I used this SqlSpatialTools library like adair but the numbers did not exactly matched after compering to ESRI ArcGIS outputs. So I found this page (http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html) and based on the javascript implementation, coded a C# class within the SQLSpatialTools framework to take care of it. Now, does anyone knows how I can submit this new class (patch files) to the project admins?

New Post: C# library for converting geom field of SQLServer into set of google latlng

Mustafa646 — Fri, 03 Feb 2012 12:55:55 GMT

Is there any C# library which can convert the geometry field into set of google latlngs? I have stored a shapefile into SQLserver table. From Sqlserver i want to read the geometry field data and want to convert this geom data into set of google latlngs in my C# application. Finally i want to save these google latlngs into another table in SQLServer.

New Post: NuGet package?

jalbert — Wed, 25 Jan 2012 21:56:12 GMT

A NuGet package would be very useful to us. I think NuGet has enough momentum behind it by now...

New Post: Conversion of State Plane to WGS84

alembic6 — Mon, 23 Jan 2012 15:26:30 GMT

Ed -

Thank you very much for the comprehensive answer and your suggestions for further avenues of pursuit! I very much appreciate your time and effort. I'll look into FME Microsoft SQL Server Data Loader Edition.

Warmest regards,

Joel

New Post: Conversion of State Plane to WGS84

edwink — Sat, 21 Jan 2012 00:41:13 GMT

@alembic6

As you have discovered, SQL Server Spatial does not offer any coordinate projection/transformation methods as part of the core offering. There are several reasons for this which I won’t delve into here. The coordinate system projection code provided in the CodePlex SQL Server Spatial Tools Project is designed as working examples of how this might be done but will not, in themselves, solve this particular problem for you. There are several alternatives. First, you could, using the CodePlex sample code as a framework, create new CLR-based T-SQL functions to do this task. The algorithms are well known and published on the Internet in a variety of places. Or you could use the Proj.4 library available from OSGeo along with some CLR-code to create some new T-SQL functions. Or you could purchase the Feature Manipulation Engine (FME) from Safe Software. If you are not familiar with FME but are interested in this approach, I would suggest that you avail yourself of the US$500 package, FME Microsoft SQL Server Data Loader Edition (please check with Safe Software for current pricing). All FME editions allow you to easily and efficiently project coordinates and perform datum transformations. Lastly, you could check with Blue Marble Geographics and see what they have to offer. I have not worked with their product but there has been some discussion with them on this subject.

Let us know if this information assists you or if you need addition guidance.

Thanks,

-Ed

From: alembic6 [email removed]
Sent: Friday, January 20, 2012 4:07 PM
To: Ed Katibah
Subject: Conversion of State Plane to WGS84 [sqlspatialtools:286893]

From: alembic6

I am trying to find the correct T-SQL code to convert X,Y State Plane coordinates (NAD83) to WGS84 coordinates. I have successfully installed SQL Spatial Tools and I have searched these discussion topics for assistance, but I am unable to assist myself. Can anyone point me in the right direction? Which function(s) should I use and how? I have a large table with X,Y values stored as integers. Thank you in advance for your assistance!

New Post: Conversion of State Plane to WGS84

alembic6 — Sat, 21 Jan 2012 00:07:02 GMT

New Post: wgs84 lat/longs to web mercator projection

kraghavk — Sun, 08 Jan 2012 08:09:44 GMT

Hi,

I am storing data as Geometry. However, the Project method seems to accept only Geography type. So, I changed the above T-SQL a bit to suit to my needs. Below is a sample script performing re-projection to one Polygon.

DECLARE @geom Geometry;
SET @geom = Geometry::STGeomFromText('POLYGON ((77.9672499974569 13.3208888689677, 77.968444442749 13.320666662852, 77.9683610916138 13.3203333218892, 77.9671388626099 13.3205555597941, 77.9672499974569 13.3208888689677))', 4326)

-- Correct the ring orientation if it is incorrect
-- SqlGeography, unlike SqlGeometry data type is very unforgiving
-- This will force a validation of the Spatial object and correct the ring orientation internally
-- For more info see: http://blogs.msdn.com/b/edkatibah/archive/2008/08/19/working-with-invalid-data-and-the-sql-server-2008-geography-data-type-part-1b.aspx
SET @geom = @geom.STUnion(@geom.STStartPoint());

DECLARE @geography Geography
SET @geography = Geography::STGeomFromText(@geom.STAsText(),@geom.STSrid)

DECLARE @mercator Projection;
SET @mercator = Projection::Mercator(0)

SELECT 
	@geom.STAsText() AS WGS84,
	@mercator.Project(@geography).STAsText() AS WebMercator

However, the values I get are totally off! Here is what I get:

Original Geometry: POLYGON ((77.968361091625411 13.320333321884391, 77.968444442754844 13.32066666285391, 77.967249997454928 13.320888868969632, 77.967138862615684 13.320555559796048, 77.968361091625411 13.320333321884391))

Re-projected Geometry: POLYGON ((1.3608046134327036 0.23460666206021522, 1.3608060681843475 0.23461264080687308, 1.3607852211811289 0.2346166262623624, 1.3607832815122705 0.23461064808050117, 1.3608046134327036 0.23460666206021522))

I am not able to figure out what I am doing wrong. Any help will be great.

Thanks in advance.

Regards,
Raghu

New Post: Conversion from ("non-vacuous") geography to geometry

StJohnKettle — Thu, 10 Nov 2011 01:19:53 GMT

There seems to be an asymmetry in the spatial tools provided. Given a geometry with an SRID we can convert to a geography with the same SRID using MakeValidGeographyFromGeometry What about the reverse process? There is a symmetric pair of functions for conversion when the items are "vacuous", which I >think< means that they have no SRID. There is also an asymmetry in the spatial aggregations: we have GeographyUnionAggregate but no GeometryUnionAggregate. Is this forthcoming or it is redundant?