A method of operating a weather radar network, the method comprising: generating respective radar beams with each of a plurality of N radars (110, 210, 310, 610) disposed at different positions within an environment; determining a plurality of respective measured reflectivities z j m (x,y) of the environment along a respective path of each of the respective radar beams from the generated respective radar beams, wherein each respective measured reflectivity has a respective cross-azimuthal resolution, characterised in that the method further comprises: determining a plurality of intrinsic reflectivities for a plurality of volume elements within the environment from the plurality of respective measured reflectivities along the respective path of each of the respective radar beams based on a solution to a problem z j m (x,y) = ∥Φ x,y Z∥ 1 + ε j (x,y), j=1,2,..., N, wherein Φ x,y is obtained from a range variant kernel formed by the effective antenna pattern and relating the reflectivity and the mean received power, Z is intrinsic two dimensional distribution of reflectivity, and ε(x,y) is noise in the measurement for the jth radar, wherein the respective intrinsic reflectivities for the respective volume elements each have a respective resolution greater than the respective cross-azimuthal resolutions of the measured reflectivities. generating respective radar beams with each of a plurality of N radars (110, 210, 310, 610) disposed at different positions within an environment; determining a plurality of respective measured reflectivities z j m (x,y) of the environment along a respective path of each of the respective radar beams from the generated respective radar beams, wherein each respective measured reflectivity has a respective cross-azimuthal resolution, characterised in that the method further comprises: determining a plurality of intrinsic reflectivities for a plurality of volume elements within the environment from the plurality of respective measured reflectivities along the respective path of each of the respective radar beams based on a solution to a problem z j m (x,y) = ∥Φ x,y Z∥ 1 + ε j (x,y), j=1,2,..., N, wherein Φ x,y is obtained from a range variant kernel formed by the effective antenna pattern and relating the reflectivity and the mean received power, Z is intrinsic two dimensional distribution of reflectivity, and ε(x,y) is noise in the measurement for the jth radar, wherein the respective intrinsic reflectivities for the respective volume elements each have a respective resolution greater than the respective cross-azimuthal resolutions of the measured reflectivities. determining a plurality of intrinsic reflectivities for a plurality of volume elements within the environment from the plurality of respective measured reflectivities along the respective path of each of the respective radar beams based on a solution to a problem z j m (x,y) = ∥Φ x,y Z∥ 1 + ε j (x,y), j=1,2,..., N, wherein Φ x,y is obtained from a range variant kernel formed by the effective antenna pattern and relating the reflectivity and the mean received power, Z is intrinsic two dimensional distribution of reflectivity, and ε(x,y) is noise in the measurement for the jth radar, wherein the respective intrinsic reflectivities for the respective volume elements each have a respective resolution greater than the respective cross-azimuthal resolutions of the measured reflectivities.