Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1806 (N44E21 L=353)   544     80        1 (1029/1784-1104/1829)   0.0 E24-E16 N42-N45
      2 1323 1677 (N41W03 L= 18)   643    123        2 (1276/1631-1371/1703) -15.9 W00-W07 N38-N43
      3  572 1574 (N16E42 L=332)   681    174        1 ( 550/1501- 591/1640) -52.1 E45-E38 N11-N20
      4 1167 1346 (N14W01 L= 15)  1067    295        1 (1111/1229-1252/1482) -49.8 E05-W08 N08-N21
      5 1878 1164 (N24W56 L= 70)   741    121        1 (1845/1125-1903/1225) -26.6 W51-W59 N21-N26
      6  246 1081 (S18E54 L=321)   793    167        0 ( 215/ 996- 265/1154) -56.3 E58-E51 S23-S14
      7  717 1052 (S12E18 L=356)   830    188        0 ( 643/ 965- 773/1140)  71.6 E21-E15 S18-S06
      8 1366  954 (S03W21 L= 36)  1518    416        2 (1253/ 787-1454/1134)  84.8 W18-W24 S14-N08
      9 1013  741 (S23W06 L= 21)   721    147        1 ( 936/ 714-1097/ 776)  16.7 W01-W11 S25-S21
     10 1777  646 (S08W62 L= 76)  1695    185        1 (1701/ 579-1838/ 721)  57.5 W56-W66 S13-S02
  

    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