Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1170 1764 (N43W40 L= 31)   349     37        0 (1147/1738-1189/1786)  63.4 W38-W42 N41-N45
      2 1420 1449 (N17W36 L= 27)   444    144        0 (1364/1423-1496/1464) -16.7 W31-W41 N14-N18
      3 1786 1144 (S10W51 L= 42)  1788    283        3 (1718/1004-1842/1256)  47.3 W45-W60 S18-S05
      4 1850 1152 (S12W58 L= 49)   363     47        3 (1834/1124-1865/1184)  76.0 W56-W60 S13-S09
      5 1412 1144 (N01W25 L= 16)   285     87        1 (1378/1119-1455/1168)   0.0 W22-W27 N00-N02
      6  391 1125 (N27E38 L=313)   315     47        1 ( 361/1105- 417/1142) -31.0 E40-E35 N25-N28
      7 1680 1009 (S14W39 L= 30)  1030    237        1 (1600/ 947-1726/1097)  54.5 W32-W43 S16-S09
      8 1301  950 (S06W13 L=  4)  1159    324        2 (1251/ 793-1343/1090)  45.0 W06-W19 S13-N00
      9  208  928 (N19E59 L=292)   500    104        0 ( 172/ 906- 250/ 963) -18.4 E63-E54 N17-N21
     10  390  840 (N11E44 L=307)   674    220        0 ( 328/ 761- 463/ 916) -63.4 E47-E40 N04-N16
     11 1716  610 (S40W37 L= 28)   491     65        1 (1683/ 582-1747/ 631) -40.6 W33-W40 S43-S37
     12 1679  524 (S45W33 L= 24)   936     84        0 (1650/ 479-1711/ 566)  71.6 W30-W37 S48-S42
     13 1544  422 (S48W19 L= 10)  1843    230        3 (1469/ 363-1637/ 478)   3.2 W10-W28 S51-S43
     14  424  400 (S16E62 L=289)   876    100        4 ( 402/ 356- 449/ 456)  35.0 E67-E57 S18-S11
     15  847  303 (S32E36 L=315)  1861    367        4 ( 795/ 128- 896/ 439)  33.7 E58-E25 S46-S23
     16  493  327 (S22E61 L=289)  2193    149        3 ( 461/ 241- 540/ 380)  22.6 E79-E53 S28-S18
  

    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