Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1608 (N46W08 L= 71)   306     29        0 ( 842/1580- 872/1633)  53.1 W06-W10 N44-N48
      2  942 1606 (N43W14 L= 77)   323     40        2 ( 913/1580- 969/1628) -36.9 W12-W16 N41-N45
      3 1564 1479 (N16W49 L=112)  1063    239        5 (1514/1384-1618/1558)  31.6 W44-W58 N10-N19
      4 1245 1461 (N25W27 L= 90)   225     57        2 (1216/1441-1284/1477) -11.3 W24-W29 N24-N26
      5 1450 1347 (N13W35 L= 98)   602    156        1 (1382/1294-1499/1377) -54.5 W30-W37 N08-N15
      6 1744 1099 (S10W48 L=111)   418     69        0 (1712/1071-1770/1129)   9.5 W44-W50 S10-S08
      7  467  987 (N18E36 L= 27)   600    144        1 ( 422/ 943- 526/1046)  24.0 E39-E30 N15-N20
      8  254  933 (N19E56 L=  7)   404     76        1 ( 226/ 905- 287/ 970) -45.0 E58-E52 N16-N21
      9  765  924 (N08E17 L= 46)   244     59        1 ( 751/ 891- 790/ 958)  56.3 E18-E15 N06-N09
     10  150  882 (N17E72 L=351)   763     86        0 ( 134/ 842- 167/ 922) -78.7 E74-E69 N14-N19
     11 1301  834 (S11W11 L= 74)   374    106        1 (1281/ 787-1327/ 883)  80.5 W09-W12 S13-S07
     12 1588  772 (S24W28 L= 91)   616    149        1 (1542/ 708-1631/ 827)  29.1 W23-W32 S26-S21
     13 1674  729 (S29W35 L= 98)   292     46        0 (1658/ 703-1689/ 760)  45.0 W33-W36 S31-S27
     14 1165  650 (S18E02 L= 61)  1094    371        5 (1005/ 624-1346/ 683) -22.8 E10-W09 S22-S13
     15 1703  589 (S41W38 L=101)   519     79        0 (1664/ 560-1728/ 619)  15.9 W33-W41 S42-S39
     16  399  453 (S13E62 L=  1)  1407    234        4 ( 377/ 364- 424/ 570)  49.1 E71-E52 S19-S04
     17  679  421 (S20E40 L= 23)   326     78        1 ( 661/ 382- 698/ 452)  71.6 E41-E37 S23-S17
  

    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