Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  645 1787 (N32E55 L=261)  1251    155        1 ( 597/1747- 711/1824)  -3.8 E62-E45 N29-N34
      2  625 1611 (N20E43 L=274)   922    122        2 ( 594/1549- 654/1675) -60.9 E46-E39 N15-N24
      3  307 1549 (N09E68 L=248)  3520    266        4 ( 219/1480- 393/1626)  42.3 E79-E59 N03-N14
      4 1426 1441 (N31W14 L=330)  2091    520        5 (1310/1211-1520/1612) -41.0 W01-W24 N19-N40
      5  878 1396 (N13E18 L=298)   251     60        0 ( 858/1360- 892/1424) -45.0 E19-E16 N10-N14
      6 1664 1296 (N29W34 L=351)   832    129        1 (1605/1246-1713/1331) -11.3 W29-W39 N27-N30
      7 1110 1156 (N05W01 L=318)   903    273        1 (1010/1084-1160/1271)  80.5 E02-W03 N00-N12
      8 1479 1101 (N12W24 L=340)   263     99        0 (1449/1059-1495/1140) -53.1 W21-W25 N09-N13
      9 1744 1005 (N14W46 L=  2)  1420    249        2 (1708/ 921-1784/1101) -45.0 W40-W51 N09-N18
     10 1671 1058 (N15W39 L=355)   280     46        1 (1646/1038-1699/1080)  45.0 W36-W40 N13-N16
     11  634  986 (S16E22 L=295)   244     61        1 ( 612/ 955- 656/1016)  90.0 E23-E20 S18-S14
     12 1610  843 (N01W40 L=356)   834    154        2 (1538/ 813-1677/ 885)   4.8 W33-W45 N00-N03
     13 1222  732 (S15W19 L=335)   217     53        0 (1198/ 706-1243/ 756)  76.0 W17-W20 S17-S13
     14 1465  714 (S10W34 L=351)   268     49        1 (1443/ 691-1484/ 733)  90.0 W32-W35 S11-S08
     15 1297  687 (S16W25 L=341)   593    158        0 (1228/ 654-1382/ 708)  36.9 W21-W30 S19-S13
     16  254  617 (S49E59 L=257)  2384    205        5 ( 166/ 571- 338/ 686)  -3.6 E78-E44 S51-S45
     17 1678  604 (S09W55 L= 11)   357     96        0 (1655/ 571-1697/ 648) -53.1 W51-W56 S11-S07
     18 1067  359 (S42W27 L=343)   812    124        2 (1014/ 328-1128/ 396)  11.3 W21-W32 S43-S40
  

    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