Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1251 1631 (N44W19 L=128)   991    232        1 (1226/1547-1274/1708)  61.7 W15-W23 N37-N50
      2 1562 1614 (N42W47 L=156)  1369    168        4 (1509/1556-1605/1663) -69.4 W41-W50 N37-N45
      3  484 1532 (N35E44 L= 66)   507     79        1 ( 464/1494- 512/1563) -51.3 E46-E41 N32-N37
      4  771 1458 (N31E19 L= 91)   418    135        2 ( 755/1399- 786/1519) -82.9 E19-E16 N26-N34
      5 1631 1443 (N30W46 L=156)   408     93        1 (1600/1413-1675/1475) -38.7 W43-W49 N27-N31
      6 1397 1308 (N21W25 L=134)   459    108        0 (1348/1261-1453/1350) -39.8 W21-W28 N18-N23
      7 1116 1270 (N19W06 L=115)   226     54        0 (1093/1253-1148/1290) -33.7 W04-W07 N17-N19
      8  902 1196 (N14E08 L=102)   314    139        0 ( 878/1131- 929/1266) -71.6 E09-E06 N09-N18
      9 1379 1217 (N15W23 L=132)   228     91        1 (1341/1206-1431/1232)   9.5 W20-W26 N14-N16
     10  681 1178 (N12E22 L= 87)   334    111        0 ( 659/1123- 716/1231)  81.9 E23-E20 N08-N15
     11 1339 1122 (N09W20 L=129)   514    230        1 (1269/1072-1459/1187) -30.3 W15-W27 N06-N13
     12  351 1024 (N02E46 L= 63)   522    149        2 ( 299/ 978- 391/1086)  40.6 E50-E42 S01-N05
     13 1695  863 (S07W46 L=156)   309     52        3 (1678/ 833-1713/ 894) -56.3 W44-W48 S08-S05
     14  609  831 (S09E27 L= 83)   217    104        0 ( 562/ 812- 666/ 857) -15.9 E29-E22 S10-S07
     15  542  643 (S22E33 L= 76)   478     98        0 ( 500/ 617- 591/ 681) -26.6 E35-E29 S23-S19
     16  370  610 (S25E48 L= 61)   952    215        1 ( 300/ 583- 434/ 643)  13.0 E55-E42 S26-S22
     17  827  359 (S42E16 L= 94)   804    152        1 ( 783/ 314- 867/ 424)  63.4 E20-E12 S46-S36
     18  706  398 (S39E25 L= 84)   306     65        1 ( 678/ 370- 733/ 422) -45.0 E27-E23 S41-S37
     19 1539  365 (S42W46 L=156)  1549    288        1 (1449/ 334-1637/ 413)   6.7 W37-W56 S44-S38
  

    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