Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  608 1640 (N46E37 L=169)   326     68        1 ( 577/1612- 641/1659)  45.0 E39-E33 N43-N47
      2 1173 1597 (N41W14 L=220)   323     76        2 (1138/1570-1203/1627) -45.0 W11-W15 N38-N43
      3 1715 1488 (N30W61 L=266)  1289    155        1 (1646/1450-1784/1527) -19.7 W54-W69 N26-N32
      4  267 1312 (N22E59 L=147)   404     49        0 ( 245/1282- 285/1334) -18.4 E62-E56 N19-N23
      5  314 1242 (N18E52 L=154)   570    128        5 ( 293/1182- 348/1303) -82.9 E54-E47 N13-N21
      6  172 1233 (N16E70 L=136)  1510    117        1 ( 148/1179- 192/1290) -50.2 E76-E65 N13-N19
      7 1775 1274 (N16W57 L=263)   400     76        0 (1734/1255-1812/1287) -14.0 W53-W61 N14-N16
      8  656  862 (S05E24 L=182)  1147    368        1 ( 505/ 761- 775/ 959) -30.5 E33-E16 S11-N00
      9  501  774 (S11E35 L=170)   280     64        0 ( 475/ 748- 531/ 798)  31.0 E37-E32 S12-S09
     10 1153  742 (S14W07 L=213)   373    143        1 (1088/ 711-1202/ 794)  23.2 W03-W10 S15-S10
     11  417  691 (S17E44 L=162)  2113    283        6 ( 295/ 605- 492/ 749) -28.1 E54-E36 S22-S13
     12  539  699 (S16E33 L=172)   278     47        1 ( 512/ 670- 565/ 721)  36.9 E35-E31 S17-S14
     13 1247  451 (S34W15 L=221)  3705    758        8 (1010/ 301-1517/ 705)  34.4 E03-W32 S46-S17
     14  734  573 (S24E21 L=185)  1344    351        1 ( 648/ 433- 828/ 687)  45.0 E30-E13 S34-S16
     15 1346  608 (S23W21 L=227)   439    123        1 (1307/ 568-1391/ 646)  45.0 W18-W23 S26-S20
     16  575  600 (S22E32 L=173)   316     39        0 ( 555/ 571- 591/ 624) -63.4 E33-E31 S24-S20
     17  774  415 (S36E21 L=185)   480     86        0 ( 733/ 392- 803/ 447)  33.7 E24-E18 S37-S33
     18  654  388 (S38E32 L=173)   398     38        1 ( 631/ 364- 678/ 410) -76.0 E34-E29 S40-S36
     19  927  197 (S58E15 L=191)   450     45        0 ( 906/ 177- 951/ 224) -80.5 E17-E12 S60-S54
  

    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