Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1542 1629 (N32W46 L=291)   361     44        0 (1517/1613-1569/1644)  11.3 W43-W48 N31-N32
      2 1656 1569 (N27W53 L=298)   439     80        0 (1632/1549-1702/1592) -31.0 W50-W57 N25-N28
      3  296 1408 (N30E51 L=195)  1209    285        4 ( 236/1365- 368/1461) -28.1 E58-E42 N26-N34
      4 1137 1424 (N22W11 L=256)   435    132        0 (1072/1406-1204/1449)   0.0 W07-W15 N20-N23
      5 1759 1262 (N07W53 L=298)   777    167        1 (1697/1229-1793/1305)  24.0 W46-W57 N05-N09
      6  339 1238 (N18E45 L=201)   314     54        1 ( 320/1210- 359/1269)  76.0 E46-E42 N16-N20
      7  158 1110 (N13E63 L=182)  2841    462       10 (  92/1036- 211/1195)  23.6 E74-E56 N08-N17
      8  998 1098 (N03E01 L=244)   176     42        1 ( 981/1070-1011/1122) -90.0 E01-W00 N01-N04
      9 1325 1059 (S02W18 L=264)   447    105        0 (1278/1036-1382/1098) -26.6 W15-W21 S03-N00
     10  155  994 (N06E65 L=180)   550     77        2 ( 133/ 964- 178/1021) -90.0 E68-E61 N03-N07
     11 1405  976 (S07W23 L=268)   369     82        2 (1382/ 934-1428/1012)  53.1 W21-W24 S09-S05
     12 1155  853 (S13W06 L=252)   357     53        1 (1117/ 829-1180/ 881)  26.6 W03-W07 S14-S10
     13 1807  785 (S22W54 L=300)  1443    196        4 (1767/ 711-1838/ 866)  66.0 W50-W57 S27-S17
     14 1169  695 (S22W06 L=251)   503    154        1 (1120/ 623-1236/ 772)  45.0 W02-W10 S26-S18
     15  981  709 (S20E06 L=239)   786    224        5 ( 931/ 661-1066/ 770) -35.5 E08-W00 S22-S15
     16  332  604 (S19E54 L=191)   583    139        3 ( 310/ 569- 355/ 640) -76.0 E57-E51 S21-S17
     17  353  510 (S26E56 L=189)   714    173        0 ( 323/ 472- 391/ 542)  18.4 E60-E51 S27-S23
     18  954  343 (S45E15 L=230)   520     84        2 ( 920/ 308- 985/ 382)  56.3 E18-E11 S47-S41
  

    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