Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1123 1831 (N59W36 L= 56)  1504    263        3 (1006/1778-1234/1876)   0.0 W19-W55 N56-N60
      2 1457 1598 (N33W45 L= 65)   469    150        0 (1389/1559-1516/1618) -40.6 W40-W48 N29-N35
      3 1358 1500 (N29W32 L= 52)   301     66        0 (1323/1481-1390/1519) -33.7 W30-W34 N27-N30
      4  481 1370 (N35E35 L=345)  1545    295        1 ( 446/1231- 529/1499)  83.3 E37-E31 N26-N43
      5  265 1244 (N28E56 L=324)  2654    308        5 ( 171/1131- 357/1347) -27.3 E74-E44 N20-N35
      6  761 1233 (N23E14 L=  7)   423    153        0 ( 697/1202- 829/1270)   5.7 E18-E08 N21-N24
      7 1594 1059 (S01W36 L= 57)   656    124        1 (1527/1039-1652/1085) -12.5 W31-W40 S02-N00
      8  505  889 (N05E35 L=345)   343    115        1 ( 480/ 845- 528/ 936)  78.7 E36-E33 N02-N07
      9 1713  868 (S15W44 L= 64)  1305    195        4 (1636/ 817-1790/ 929)  11.3 W36-W51 S17-S12
     10 1595  852 (S13W34 L= 54)   277     30        0 (1574/ 830-1616/ 876) -56.3 W32-W35 S14-S11
     11  890  640 (S15E14 L=  6)  1191    410        0 ( 766/ 557- 972/ 793)  45.0 E23-E08 S19-S05
     12  566  643 (S11E35 L=345)  1115    190        3 ( 546/ 542- 602/ 718)  65.6 E39-E31 S17-S06
     13  719  648 (S12E25 L=355)   263     55        2 ( 685/ 631- 753/ 663) -26.6 E26-E22 S13-S10
     14  289  617 (S10E61 L=319)   496     91        2 ( 264/ 592- 329/ 651) -38.7 E64-E57 S12-S08
     15  286  536 (S16E67 L=313)   810    109        1 ( 259/ 502- 302/ 589)  33.7 E73-E61 S18-S12
     16 1507  368 (S47W29 L= 49)   947     86        2 (1463/ 345-1551/ 401)   6.3 W24-W34 S49-S44
     17 1365  306 (S49W15 L= 35)   729     80        1 (1339/ 260-1402/ 343)  90.0 W12-W19 S54-S45
     18 1291  252 (S53W08 L= 28)  1246    116        4 (1245/ 216-1331/ 296) -45.0 W02-W12 S58-S48
     19 1287  190 (S61W07 L= 27)   691     64        1 (1252/ 176-1320/ 213) -24.0 W02-W11 S64-S58
  

    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