Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1633 (N37W07 L= 33)   281     29        1 (1029/1614-1075/1658) -56.3 W05-W08 N35-N38
      2 1737 1456 (N21W55 L= 82)   348     65        1 (1715/1433-1761/1475) -33.7 W53-W57 N19-N21
      3 1773 1405 (N17W57 L= 84)   852    176        4 (1740/1356-1810/1455)  45.0 W53-W61 N14-N20
      4 1003 1401 (N21W01 L= 28)   252     66        0 ( 968/1374-1044/1425) -45.0 W00-W03 N19-N22
      5  193 1207 (N15E59 L=327)   598    144        1 ( 167/1145- 214/1246) -63.4 E62-E56 N11-N17
      6 1244 1221 (N08W15 L= 41)   207     63        0 (1203/1210-1276/1231) -14.0 W12-W16 N07-N09
      7  438 1047 (N04E37 L=349)  1166    251        1 ( 375/ 922- 476/1152) -72.9 E41-E34 S03-N10
      8  194 1002 (N03E59 L=327)  1958    324        7 ( 130/ 957- 292/1050) -16.4 E67-E49 N00-N06
      9  288  779 (S11E53 L=334)   791    198        3 ( 238/ 727- 333/ 832)  24.4 E59-E47 S13-S08
     10  285  596 (S22E59 L=328)   402     89        1 ( 264/ 566- 307/ 631)  76.0 E60-E56 S24-S20
     11  337  581 (S24E54 L=333)   449     92        1 ( 313/ 531- 370/ 611) -80.5 E55-E51 S27-S21
     12  989  403 (S41E08 L= 19)   819    203        1 ( 905/ 361-1094/ 462)  18.4 E15-W00 S44-S37
  

    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