Tengo dos tablas espaciales en PostgreSQL/PostGIS llamadas departamentos y distritos. Cada tabla contiene su geometría respectiva.

Yo quiero pasar un punto y que me diga de que departamento es con su respectivo distrito. Siempre voy a querer que me salga el departamento, sabiendo que algunos distritos no están cargados. Como podría hacer con funciones espaciales.

Son dos tablas

    CREATE TABLE `departamentos` (
      `departamento_id` int(2) NOT NULL AUTO_INCREMENT,
      `departamento_nombre` varchar(250) DEFAULT NULL,
      `departamento_capital` varchar(250) DEFAULT NULL,
      `geom` geometry NOT NULL,
      PRIMARY KEY (`departamento_id`),
      UNIQUE KEY `i_departamentoids` (`departamento_id`),
      SPATIAL KEY `i_geomidx` (`geom`)
    ) ENGINE=MyISAM AUTO_INCREMENT=18 DEFAULT CHARSET=utf8;
    
    
    CREATE TABLE `distritos` (
      `distrito_id` int(3) NOT NULL AUTO_INCREMENT,
      `distrito_nombre` varchar(250) DEFAULT NULL,
      `geom` geometry NOT NULL,
      PRIMARY KEY (`distrito_id`),
      UNIQUE KEY `i_distritoids` (`distrito_id`),
      SPATIAL KEY `i_geomidx` (`geom`)
    ) ENGINE=MyISAM AUTO_INCREMENT=103 DEFAULT CHARSET=utf8;

Muchas gracias