Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1512 1764 (N49W66 L=169)   602     39        1 (1490/1735-1531/1788)  15.9 W60-W75 N47-N51
      2  531 1546 (N41E39 L= 64)   448     59        2 ( 514/1514- 555/1578) -51.3 E41-E35 N38-N43
      3 1041 1558 (N39W05 L=108)   529    108        1 ( 980/1533-1097/1575)  -6.3 W00-W09 N37-N40
      4 1237 1459 (N31W18 L=121)   319    108        0 (1201/1414-1264/1500) -71.6 W16-W19 N27-N33
      5  514 1306 (N24E34 L= 68)   568     96        1 ( 495/1251- 535/1355) -71.6 E36-E32 N21-N27
      6 1826 1243 (N11W63 L=165)   566     97        1 (1807/1196-1851/1285) -71.6 W60-W64 N07-N13
      7 1784 1242 (N11W58 L=160)   555     75        1 (1758/1203-1805/1278) -78.7 W55-W59 N08-N13
      8 1738 1132 (N05W51 L=153)   330     42        0 (1713/1115-1765/1149) -21.8 W48-W53 N04-N06
      9 1865  973 (S06W63 L=166)  1167    163        5 (1812/ 914-1914/1036) -36.9 W56-W71 S10-S01
     10 1779  783 (S17W54 L=157)  4572    670        8 (1704/ 602-1861/ 997) -85.6 W46-W64 S29-S03
     11  198  887 (S01E63 L= 39)  1438    254        0 ( 174/ 775- 234/ 991)  58.4 E66-E57 S08-N05
     12 1818  869 (S12W58 L=160)   426     60        1 (1785/ 848-1850/ 895)   0.0 W53-W61 S12-S10
     13 1395  797 (S12W22 L=125)   280    149        0 (1339/ 736-1431/ 842)  45.0 W18-W24 S15-S09
     14  900  764 (S11E09 L= 93)   266     92        2 ( 871/ 731- 948/ 792)  31.0 E11-E06 S12-S09
     15  406  653 (S16E46 L= 56)   365     86        0 ( 379/ 613- 425/ 692) -90.0 E47-E44 S18-S13
     16  866  462 (S31E15 L= 87)  1872    480        3 ( 688/ 391- 996/ 601)  23.7 E29-E04 S36-S21
     17  426  495 (S27E51 L= 52)   500     78        2 ( 388/ 477- 468/ 517)  -7.1 E54-E46 S28-S25
  

    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