Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1682 (N30W56 L=190)  1938    315        3 (1333/1652-1587/1709) -20.0 W43-W67 N24-N35
      2 1757 1531 (N12W70 L=205)   851     98        1 (1723/1498-1792/1566) -78.7 W67-W73 N09-N14
      3 1326 1467 (N23W31 L=166)   681    128        1 (1254/1440-1394/1489) -32.0 W27-W35 N19-N25
      4  783 1161 (N20E11 L=124)  1237    318        3 ( 639/1095- 930/1234) -40.9 E18-E03 N13-N28
      5 1896 1128 (S13W65 L=199)  4654    313       12 (1822/1052-1955/1233) -61.7 W55-W78 S19-S06
      6 1690 1131 (S06W43 L=178)   256     99        0 (1662/1086-1724/1162)  15.9 W39-W46 S07-S04
      7  229 1026 (N24E57 L= 78)   811    185        3 ( 182/ 958- 275/1111) -47.7 E62-E51 N19-N30
      8 1243  983 (S01W11 L=146)   259     82        1 (1210/ 955-1277/1010) -53.1 W09-W12 S03-N01
      9 1666  841 (S21W35 L=169)   983    150        4 (1632/ 773-1703/ 914)  63.4 W31-W38 S25-S17
     10  502  716 (N02E39 L= 95)   755    189        0 ( 481/ 634- 529/ 800)  58.0 E41-E35 S03-N05
     11  950  710 (S09E12 L=122)   327    164        0 ( 932/ 651- 988/ 782)  82.9 E13-E09 S12-S04
     12 1131  677 (S15E03 L=132)   767    178        3 (1079/ 610-1198/ 744)  50.2 E05-W02 S19-S11
     13 1558  658 (S29W24 L=158)   916    169        3 (1499/ 600-1599/ 735)  29.1 W18-W27 S32-S25
     14  854  467 (S21E26 L=109)   229     81        0 ( 843/ 427- 867/ 514)  63.4 E27-E23 S23-S17
     15 1450  299 (S55W09 L=143)   402     54        1 (1418/ 283-1482/ 317)  -8.1 W05-W12 S56-S52
     16  960  243 (S40E32 L=103)  1323    160        4 ( 886/ 202-1021/ 292)  10.6 E40-E24 S42-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