Median Income of each NYC Taxi Zone

Hi everyone, two posts in one night here, but I just wanted to make anyone's lives easier if they were curious about this or needed this for any type of project.

Basically, as we all know, New York has 256 taxi zones. However, if you're trying to figure out the median income of each taxi zone, it's quite hard to do because there is no data. There is only data for blocks, census tracts, municipalities, etc. Therefore, we need to use weighted averages when overlaying the different shapes of NTAs and Taxi zones. The code I used is below, for anyone interested in running this regression themselves.

// Load Taxi Zones and NTAs from your assets

var taxiZones = ee.FeatureCollection('projects/ee-YOURNAMAE/assets/FILENAME');

var ntas = ee.FeatureCollection('projects/YOURNAME/assets/FILENAME');

// Assuming NTAs have a property 'median_income' which contains the median income data

var ntaWithIncome = ntas.filter(ee.Filter.notNull(['median_income']));

// Function to calculate weighted median income for each taxi zone

var calculateWeightedIncome = function(zone) {

  // Get the area of the taxi zone

  var zoneArea = zone.geometry().area();

  // Intersect the taxi zone with the NTAs

  var intersections = ntaWithIncome.map(function(nta) {

    var intersection = nta.geometry().intersection(zone.geometry(), ee.ErrorMargin(1));

    var intersectionArea = intersection.area();

    // Calculate the proportion of the taxi zone that overlaps with this NTA

    var weight = intersectionArea.divide(zoneArea);

    // Get the median income from the NTA

    var medianIncome = nta.get('median_income');

    // Multiply the income by the weight (proportion of the area)

    var weightedIncome = ee.Number(medianIncome).multiply(weight);

    // Set the weighted income to the taxi zone

    return zone.set('weighted_income', weightedIncome);

  });

  // Sum up the weighted incomes for all intersecting NTAs

  var totalWeightedIncome = intersections.aggregate_sum('weighted_income');

  // Return the taxi zone with the total weighted income and the location_i

  return zone.set('median_income_estimated', totalWeightedIncome);

};

// Apply the function to all taxi zones

var taxiZonesWithIncome = taxiZones.map(calculateWeightedIncome);

Again if anyone is interested in the data directly pop me a message and I'll respond ASAP.