Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1449 1630 (N26W47 L=137)   821    146        0 (1407/1564-1506/1667)  13.0 W40-W53 N23-N28
      2 1616 1524 (N16W53 L=143)   601    117        0 (1596/1483-1645/1585)  50.2 W50-W57 N13-N19
      3 1221 1535 (N27W26 L=117)   414     91        0 (1168/1511-1268/1552) -38.7 W23-W28 N24-N29
      4 1097 1338 (N19W12 L=102)  1004    338        1 ( 992/1270-1234/1441) -56.3 W08-W18 N12-N27
      5  771 1302 (N25E08 L= 82)  1816    574        3 ( 530/1193-1000/1366) -29.1 E24-W03 N14-N32
      6 1863 1217 (S07W62 L=153)  2406    387        2 (1805/1026-1898/1357)  52.6 W53-W70 S17-N01
      7  642 1157 (N20E19 L= 71)   212    101        0 ( 591/1146- 691/1166) -26.6 E22-E16 N18-N21
      8 1718  923 (S19W41 L=131)  1426    274        2 (1606/ 829-1823/ 980)   5.2 W30-W52 S22-S17
      9 1912  876 (S29W64 L=154)  1163     98        2 (1875/ 848-1964/ 896)  -3.0 W57-W76 S30-S26
     10  970  784 (S08E09 L= 81)   532    196        2 ( 917/ 725-1006/ 859)  45.0 E12-E05 S11-S03
     11  586  544 (S13E39 L= 51)  4386    851       11 ( 363/ 363- 859/ 756) -54.5 E48-E25 S27-N02
     12  162  694 (N03E74 L= 16)  1087    136        1 ( 130/ 664- 214/ 734)  -7.6 E83-E66 N01-N05
     13  442  584 (S08E48 L= 42)   694    135        2 ( 408/ 537- 474/ 650)  82.9 E51-E45 S11-S03
     14  368  449 (S15E63 L= 27)   944    124        2 ( 308/ 431- 431/ 477) -22.6 E68-E56 S16-S11
     15  571  270 (S30E59 L= 31)  3170    261        7 ( 498/ 168- 652/ 368)  39.0 E76-E51 S39-S22
     16 1497  280 (S61W18 L=108)   481     44        1 (1474/ 254-1519/ 309) -52.1 W14-W22 S65-S56
  

    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