Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1618 (N22E64 L=198)   609     94        2 ( 381/1591- 443/1652) -13.0 E70-E57 N20-N23
      2  359 1473 (N12E58 L=204)  1085    157        1 ( 338/1412- 389/1546) -66.0 E61-E52 N07-N16
      3 1634 1288 (N27W36 L=298)   566     82        2 (1607/1253-1662/1336) -59.0 W33-W38 N25-N30
      4  549  830 (S23E27 L=235)   323    139        0 ( 501/ 783- 583/ 889)  63.4 E29-E24 S26-S18
      5  243  803 (S31E56 L=206)   951    175        2 ( 189/ 747- 295/ 880) -30.3 E62-E49 S33-S26
      6 1046  696 (S21W08 L=270)   458    133        3 ( 988/ 668-1106/ 732)  35.5 W04-W11 S23-S17
      7  703  573 (S36E12 L=250)   593    144        1 ( 646/ 541- 772/ 617) -11.3 E16-E06 S37-S34
      8 1574  487 (S22W53 L=315)   787    154        1 (1528/ 437-1631/ 527)  66.8 W51-W57 S25-S18
      9 1657  443 (S22W66 L=327)  2105    235        7 (1578/ 381-1729/ 513)  70.0 W60-W74 S27-S16
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y