Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  691 1640 (N46E28 L=168)   445     87        0 ( 659/1614- 737/1655)  15.9 E31-E23 N43-N47
      2 1357 1546 (N36W28 L=224)   420     66        1 (1321/1520-1389/1568) -11.3 W24-W30 N34-N38
      3 1756 1514 (N31W69 L=266)   696     55        0 (1728/1497-1788/1530) -14.0 W64-W75 N29-N32
      4  357 1298 (N22E49 L=148)   445    125        0 ( 324/1242- 380/1324) -29.7 E52-E45 N18-N23
      5  236 1201 (N15E60 L=137)  1336    184        2 ( 219/1105- 259/1281) -73.3 E63-E56 N09-N19
      6  696  950 (N00E21 L=176)   280     81        3 ( 654/ 924- 729/ 980)  11.3 E23-E18 S01-N01
      7  841  837 (S07E12 L=185)   807    237        0 ( 757/ 767- 930/ 917) -42.0 E16-E06 S11-S01
      8 1466  866 (S07W28 L=224)   315    104        0 (1434/ 832-1484/ 896)  53.1 W25-W29 S09-S05
      9 1313  768 (S13W18 L=214)   515    171        3 (1263/ 716-1350/ 836)  50.2 W14-W20 S15-S08
     10  537  659 (S18E34 L=162)  2571    453       10 ( 396/ 526- 695/ 774) -41.6 E44-E24 S27-S10
     11  885  593 (S23E10 L=186)  1093    253        1 ( 808/ 492- 976/ 676)  45.0 E16-E04 S30-S17
     12 1494  536 (S29W33 L=230)  1538    306        3 (1424/ 448-1604/ 661)  54.5 W28-W40 S35-S21
     13 1452  588 (S25W29 L=225)   527    142        0 (1385/ 541-1502/ 623)  32.0 W24-W32 S28-S23
     14  828  387 (S38E18 L=179)  1074    251        2 ( 713/ 348- 941/ 414)   0.0 E26-E08 S41-S36
     15 1177  327 (S45W11 L=207)   863    123        0 (1142/ 284-1219/ 386) -56.3 W07-W14 S48-S39
     16 1229  330 (S45W15 L=212)   332     42        1 (1210/ 308-1256/ 354) -45.0 W13-W17 S46-S42
     17  835  219 (S55E24 L=172)   406     35        1 ( 815/ 197- 854/ 244)  39.8 E27-E21 S57-S51
     18 1093  178 (S61W04 L=201)   468     57        1 (1076/ 157-1125/ 202)  38.7 W02-W08 S64-S57
  

    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