*** STEP 1. Specify the folder location that SAS will use to write HTML output files. ***;
ods html path="R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818"                            /* Need to change */
        gpath="R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818"(URL=none);                 /* Need to change */  
ods html newfile=proc; 


*** STEP 2. READ IN A TAB-DELIMITED FILE THAT CONTAINS THE DATABASE FIELDS OR PATIENT ID/TUMOR RECORD NUMBER         ***;
*** THAT YOU WILL NEED AT THE END OF THE PROCESS TO UPDATE YOUR CENTRAL REGISTRY DATABASE.                           ***;
*** IN ADDITION TO THESE, YOU ONLY NEED THE LATITUDE AND LONGITUDE VALUES FROM PREVIOUS HIGH-QUALITY GEOCODED CASES. ***;
*** SAS CODE IS NOT INCLUDED FOR THIS STEP, AS IT WILL DEPEND ON YOUR REGISTRY SOFTWARE.                             ***;
*** LATITUDE AND LONGITUDE FIELDS SHOULD HAVE THOSE EXACT NAMES.                                                     ***;
*** THE CODE FOR READING IN THE TAB-DELIMITED FILE WILL GO HERE.                                                     ***;

proc contents data="R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818\NM_NPI300_lat_long";      /* Need to change */
run;

proc print data="R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818\NM_NPI300_lat_long" noobs n;  /* Need to change */
run;


*** STEP 3. Create the dataset for PIP in SAS. This requires X and Y numeric variables for lat and long. ***;
data forPIP;
   set "R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818\NM_NPI300_lat_long";                  /* Need to change */
   X = longitude + 0;
   Y = latitude  + 0;
run;


*** STEP 4. Retrieve 2020 census block-level tiger line file for your state from: 
https://www.census.gov/geographies/mapping-files/time-series/geo/tiger-line-file.html

Choose the FTP Archive and navigate to the TABBLOCK20 folder.

https://www2.census.gov/geo/tiger/TIGER2024/TABBLOCK20/
***;


*** STEP 5 - APPEND CENSUS 2020 GEOGRAPHY TO PREVIOUSLY GEOCODED POINTS ***;
proc mapimport
    datafile='R:\EPI\CANCER\NAACCR\GIS_Work_Group\PIP_webinar_20250818\tl_2020_35_tabblock20.shp'      /* Need to change */
    out=NM_blk_2020 create_id_;                                                                        /* Need to change state */   
run;
proc contents data=NM_blk_2020;                                                                         /* Need to change state */   
run;
proc ginside
   data=forPIP        /* Need to change - this is the name of the SAS dataset created from reading in your case data as a tab-delimited file */  
   map=NM_blk_2020                                                                              /* Need to change state */      
   out=NM_PIP_2020 dropmapvars includeborder;                                                   /* Need to change state */  
   id GEOID20;
run;
proc contents data=NM_PIP_2020 position;                                                        /* Need to change state */    
run;
proc print data=NM_PIP_2020 noobs n;                                                            /* Need to change state */         
run;


*** NOW YOU JUST NEED TO PULL THE CENSUS BITS OUT OF THE GEOID20 FIELD ***;
*** SAS CODE FOR THIS STEP IS NOT PROVIDED, AS IT WILL DEPEND ON THE CENTRAL REGISTRY SOFTWARE YOU USE. ***;
/*
EXAMPLE:
350130018071007

35 013 001807 1007
ST CTY TRACT  BLOCK (BLOCK GROUP IS THE FIRST DIGIT OF BLOCK, "1" IN THIS CASE)
*/