topic Re: spatial join in SAS Programming

spatial join

hanfei28 — Tue, 30 Mar 2021 08:40:25 GMT

Hi, I have two files. One file contains LocationID and associated latitude and longitude for each LocationID. Another file is Census Tiger file (tabblock10.shp) in Shapefile format. File tabblock10 contains blockID and associated latitudes and longitudes (polygon which defines the block). I want to match these two files to determine the blockID for each LocationID (find block each location falls into based on the shape file). How to do such a spatial join in SAS? Thanks.

Re: spatial join

PeterClemmensen — Tue, 30 Mar 2021 09:19:56 GMT

Please provide sample data and your desired result. Makes it much easier to provide a usable code answer.

Re: spatial join

PaalNavestad — Tue, 30 Mar 2021 10:36:43 GMT

This procedure does the trick I believe

SAS Help Center: Syntax: PROC GINSIDE PROC GINSIDE Statement

You must import the shapefile using the appropriate shape importer,

Re: spatial join

sbxkoenk — Tue, 30 Mar 2021 10:42:22 GMT

Here's an example.

The GINSIDE procedure is used to find out if the location (customer) is inside or outside the polygon.

Example 2: Mapping and Annotating Values from the GINSIDE Procedure

https://go.documentation.sas.com/?cdcId=pgmsascdc&cdcVersion=9.4_3.5&docsetId=grmapref&docsetTarget=p1ig1hwfnyblgln17ensk8zous3k.htm&locale=en

(Example 1 might be simpler.)

The above page has extensive info on the how and why of each individual statement in the program.

To import your TIGER file: go to

https://go.documentation.sas.com/?docsetId=helpcenterwlcm&docsetTarget=home.htm&docsetVersion=1.0&locale=en

and search for TIGER in the search field (click on magnifying glass).

Cheers,

Koen

Re: spatial join

hanfei28 — Tue, 30 Mar 2021 14:37:37 GMT

sample_data.csv

tl_2020_24_tabblock10.cpg

tl_2020_24_tabblock10.dbf

tl_2020_24_tabblock10.prj

tl_2020_24_tabblock10.shp

tl_2020_24_tabblock10.shp.ea.iso.xml

tl_2020_24_tabblock10.shp.iso.xml

tl_2020_24_tabblock10.shx

The desired result is to find the BlockID the set of latitude and longitude
(in the sample data) falls into.