A method for hydraulic fracture characterization using microseismic datapoints comprising: identifying a location of all microseismic datapoints; setting an initial orientation for a plane, the initial orientation comprising a strike angle and a dip angle; identifying the principal strike orientations of the average microseismic datapoints that relate to one or more sets of planes and their orientation uncertainty using a computer processor determining a dip of each set of planes considering a contact weight between each microseismic datapoint and each plane within each set of planes; determining a location of each plane within each set of planes; and performing a statistical analysis for each set of planes with a common identified principal strike orientation to determine an uncertainty related to an orientation, spacing and dip of a respective set of subsurface fractures, the method being characterized in that the principal strike orientations of the average microseismic datapoints that relate to one or more sets of planes and their orientation uncertainty are identified by: creating the plane that has a normal vector and passes through the average location of all microseismic datapoints for the orientation of the plane; calculating a perpendicular distance between each microseismic datapoint and the plane; creating a set of planes parallel to the plane wherein each plane in the set of planes intersects a volume represented by a location uncertainty of each microseismic datapoint; counting the number of planes in the set of planes and the number of microseismic datapoints that contact the planes in the set of planes; calculating the average microseismic datapoints that contact each plane in the set of planes for each orientation of the planes in the set of planes; incrementing the strike angle or the dip angle of the plane by 1 degree; and repeating each of the prior steps until the strike angle equals 0 degrees or 360 degrees and the dip angle equals 0 degrees or 90 degrees. identifying a location of all microseismic datapoints; setting an initial orientation for a plane, the initial orientation comprising a strike angle and a dip angle; identifying the principal strike orientations of the average microseismic datapoints that relate to one or more sets of planes and their orientation uncertainty using a computer processor determining a dip of each set of planes considering a contact weight between each microseismic datapoint and each plane within each set of planes; determining a location of each plane within each set of planes; and performing a statistical analysis for each set of planes with a common identified principal strike orientation to determine an uncertainty related to an orientation, spacing and dip of a respective set of subsurface fractures, the method being characterized in that the principal strike orientations of the average microseismic datapoints that relate to one or more sets of planes and their orientation uncertainty are identified by: creating the plane that has a normal vector and passes through the average location of all microseismic datapoints for the orientation of the plane; calculating a perpendicular distance between each microseismic datapoint and the plane; creating a set of planes parallel to the plane wherein each plane in the set of planes intersects a volume represented by a location uncertainty of each microseismic datapoint; counting the number of planes in the set of planes and the number of microseismic datapoints that contact the planes in the set of planes; calculating the average microseismic datapoints that contact each plane in the set of planes for each orientation of the planes in the set of planes; incrementing the strike angle or the dip angle of the plane by 1 degree; and repeating each of the prior steps until the strike angle equals 0 degrees or 360 degrees and the dip angle equals 0 degrees or 90 degrees. creating the plane that has a normal vector and passes through the average location of all microseismic datapoints for the orientation of the plane; calculating a perpendicular distance between each microseismic datapoint and the plane; creating a set of planes parallel to the plane wherein each plane in the set of planes intersects a volume represented by a location uncertainty of each microseismic datapoint; counting the number of planes in the set of planes and the number of microseismic datapoints that contact the planes in the set of planes; calculating the average microseismic datapoints that contact each plane in the set of planes for each orientation of the planes in the set of planes; incrementing the strike angle or the dip angle of the plane by 1 degree; and repeating each of the prior steps until the strike angle equals 0 degrees or 360 degrees and the dip angle equals 0 degrees or 90 degrees.