Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1235 1749 (N41W44 L=260)   760    120        1 (1186/1701-1273/1790)  12.1 W37-W51 N39-N42
      2 1696 1536 (N13W63 L=278)  2129    190        3 (1645/1459-1744/1614)  73.3 W56-W67 N08-N18
      3 1602 1306 (N05W42 L=258)   850    182        0 (1516/1244-1667/1351) -56.3 W38-W45 N00-N09
      4 1472 1219 (N05W31 L=246)   237     81        0 (1452/1173-1492/1257)  45.0 W28-W33 N02-N06
      5 1819 1035 (S16W52 L=267)   885    175        2 (1755/ 967-1864/1064) -48.4 W45-W56 S21-S12
      6 1813  919 (S23W49 L=264)   409     98        0 (1772/ 892-1842/ 959)   7.1 W44-W52 S24-S21
      7  319  892 (N16E49 L=167)   316     75        0 ( 288/ 876- 346/ 923)   9.5 E51-E45 N14-N17
      8  659  881 (N08E24 L=191)   223     69        0 ( 621/ 864- 690/ 904)   0.0 E26-E21 N06-N08
      9 1060  789 (S07E04 L=211)   677    223        1 ( 998/ 698-1130/ 884) -85.6 E05-E02 S13-N00
     10  457  779 (N07E41 L=175)   396    102        0 ( 399/ 767- 508/ 791)  -7.1 E44-E36 N05-N08
     11  705  653 (S06E28 L=187)   249    114        0 ( 651/ 636- 754/ 679)   0.0 E31-E24 S06-S05
     12  621  554 (S10E37 L=178)  2256    342        7 ( 515/ 435- 717/ 651) -65.2 E45-E31 S18-S04
     13  767  515 (S15E29 L=186)   273     80        0 ( 724/ 499- 796/ 537) -14.0 E32-E27 S16-S14
     14  930  396 (S27E25 L=191)  1630    248        3 ( 865/ 309-1001/ 512)  29.4 E32-E16 S32-S21
     15  603  405 (S18E46 L=170)  3113    558       10 ( 471/ 276- 792/ 491)   9.5 E70-E28 S25-S12
     16 1088  349 (S35E18 L=198)  1195    211        2 ( 986/ 307-1173/ 401)   3.6 E26-E10 S37-S33
     17  719  209 (S35E54 L=161)   467     14        0 ( 699/ 191- 739/ 229)  23.2 E57-E50 S36-S33
  

    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