Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  887 1697 (N49E15 L=121)  3298    533        5 ( 787/1564-1004/1783)  56.3 E22-E04 N38-N57
      2 1456 1664 (N47W38 L=174)   690    131        1 (1412/1637-1498/1684)  16.7 W32-W42 N44-N48
      3 1317 1604 (N42W23 L=159)   815    156        3 (1276/1546-1350/1660)  56.3 W19-W27 N37-N46
      4 1593 1605 (N42W49 L=184)   495     54        1 (1568/1579-1619/1625) -33.7 W46-W50 N39-N43
      5  301 1548 (N34E63 L= 73)   703    105        3 ( 278/1521- 329/1569)  26.6 E65-E59 N31-N35
      6 1673 1549 (N37W54 L=190)   442     74        0 (1647/1535-1701/1562)   9.5 W50-W56 N36-N38
      7  467 1454 (N29E43 L= 93)   439    149        0 ( 456/1385- 488/1509) -82.9 E44-E41 N24-N32
      8 1334 1430 (N29W21 L=157)   293     50        1 (1313/1402-1355/1464) -76.0 W20-W22 N27-N31
      9 1021 1306 (N21E01 L=135)   429     79        1 ( 984/1277-1077/1337) -45.0 E03-W02 N18-N22
     10 1014 1228 (N16E01 L=135)   367    157        1 ( 924/1213-1077/1259)   6.3 E06-W03 N14-N17
     11 1612 1071 (N06W39 L=175)   415     86        1 (1580/1038-1651/1109) -38.7 W36-W42 N04-N08
     12 1312  893 (S04W18 L=154)   294     88        0 (1297/ 845-1330/ 943)  80.5 W17-W19 S07-S01
     13  296  855 (S10E52 L= 84)   317    127        0 ( 259/ 831- 337/ 894) -20.6 E55-E47 S11-S07
     14  357  765 (S15E47 L= 89)   510     83        1 ( 324/ 736- 391/ 789)  18.4 E49-E43 S17-S13
     15 1040  730 (S15W02 L=137)   276     62        0 (1021/ 695-1057/ 765) -78.7 W00-W02 S17-S12
     16 1017  645 (S21W00 L=136)   265     36        1 ( 992/ 620-1038/ 670) -26.6 E01-W01 S22-S19
     17  507  395 (S41E43 L= 93)  1507    221        0 ( 450/ 354- 557/ 456)  31.0 E51-E35 S45-S35
     18 1336  373 (S41W27 L=163)  1586    352        3 (1192/ 328-1508/ 428)  13.5 W15-W41 S44-S36
  

    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