A method, implemented by a system coupled with a wireless telecommunication network, of estimating flows of persons that gathered at an Area of Interest for attending a public happening during a time interval on a day, wherein said Area of Interest is defined by an Area of Interest center and an Area of Interest radius and is covered by a mobile telecommunication network having a plurality of communication stations each of which is adapted to manage communications of user equipment in one or more served areas in which the mobile telecommunication network is subdivided, the method comprising the steps of: a) defining a plurality of calculated radius values of the Area of Interest radius, and, for each calculated radius value: b) identifying (1210,1212) a first number of user equipment associated with at least one event record of a corresponding event of interaction occurred between the user equipment and the mobile communication network during the time interval on the day within the Area of Interest; c) identifying (1214,1216) a second number of user equipment associated with at least one event record of a corresponding event of interaction occurred between the user equipment and the mobile communication network during the time interval for each day of a predetermined number of previous days preceding the day within the Area of Interest; d) combining (1218,1220) the first number of user equipment and the second numbers of user equipment for obtaining a statistical quantity; e) detecting an occurrence of the public happening if the statistical quantity reaches a certain threshold; f) computing an optimum radius value of the Area of Interest radius as the average of the calculated radius values within which the public happening is detected; g) identifying persons (1238 - 1268), owners of user equipment, that gathered for attending at the public happening within the Area of Interest having the Area of Interest radius equal to the optimum radius values during the time interval on the day within the Area of Interest based on a first time fraction indicating a probability that a user equipment has been in the Area of Interest during the time interval on the day and on a second time fraction indicating that the user equipment has been in the Area of Interest during the previous days for each user equipment identified at step b); h) computing at least one matrix (700m',n) accounting for movements of persons identified at step g) within a Region of Interest comprising the Area of Interest to the Area of Interest during at least one observation time period comprising the time interval, and i) computing at least one matrix (800m",n) accounting for movements of persons identified at step g) within a Region of Interest comprising the Area of Interest from the Area of Interest during at least one observation time period comprising the time interval. a) defining a plurality of calculated radius values of the Area of Interest radius, and, for each calculated radius value: b) identifying (1210,1212) a first number of user equipment associated with at least one event record of a corresponding event of interaction occurred between the user equipment and the mobile communication network during the time interval on the day within the Area of Interest; c) identifying (1214,1216) a second number of user equipment associated with at least one event record of a corresponding event of interaction occurred between the user equipment and the mobile communication network during the time interval for each day of a predetermined number of previous days preceding the day within the Area of Interest; d) combining (1218,1220) the first number of user equipment and the second numbers of user equipment for obtaining a statistical quantity; e) detecting an occurrence of the public happening if the statistical quantity reaches a certain threshold; f) computing an optimum radius value of the Area of Interest radius as the average of the calculated radius values within which the public happening is detected; g) identifying persons (1238 - 1268), owners of user equipment, that gathered for attending at the public happening within the Area of Interest having the Area of Interest radius equal to the optimum radius values during the time interval on the day within the Area of Interest based on a first time fraction indicating a probability that a user equipment has been in the Area of Interest during the time interval on the day and on a second time fraction indicating that the user equipment has been in the Area of Interest during the previous days for each user equipment identified at step b); h) computing at least one matrix (700m',n) accounting for movements of persons identified at step g) within a Region of Interest comprising the Area of Interest to the Area of Interest during at least one observation time period comprising the time interval, and i) computing at least one matrix (800m",n) accounting for movements of persons identified at step g) within a Region of Interest comprising the Area of Interest from the Area of Interest during at least one observation time period comprising the time interval.