spatial join

[hanfei28]

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.

[PeterClemmensen]

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

[hanfei28]

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.

[PaalNavestad · Accepted solution]

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,

[sbxkoenk]

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=...

(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&lo...

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

Cheers,

Koen