Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  785 1705 (N57W05 L= 42)   428     53        1 ( 762/1677- 807/1734)  51.3 W02-W07 N54-N59
      2  439 1542 (N54E38 L=359)   403     45        1 ( 408/1520- 461/1563)   8.1 E43-E35 N52-N55
      3 1368 1534 (N27W38 L= 76)   277     97        0 (1334/1517-1415/1552) -26.6 W34-W40 N24-N28
      4 1206 1332 (N20W19 L= 57)   226    137        2 (1160/1269-1239/1360) -74.1 W16-W20 N15-N22
      5  945 1158 (N17E01 L= 36)   883    225        2 ( 858/1092-1010/1223) -56.3 E05-W02 N12-N22
      6 1763  926 (S19W46 L= 84)  1688    362        1 (1703/ 779-1814/1086)  47.0 W39-W54 S25-S10
      7  731 1053 (N15E17 L= 21)   245     88        2 ( 700/1017- 755/1085) -90.0 E18-E14 N12-N17
      8 1006  997 (N06E02 L= 36)   446    139        0 ( 970/ 936-1029/1058)  50.2 E05-W00 N03-N09
      9  292  970 (N19E52 L=346)   377     79        0 ( 282/ 925- 302/1012) -90.0 E52-E50 N16-N21
     10 1863  893 (S25W58 L= 96)   847    108        0 (1822/ 848-1902/ 939)   9.5 W52-W64 S26-S22
     11  907  814 (S02E12 L= 26)   355    104        2 ( 887/ 769- 927/ 867)  78.7 E12-E10 S04-N01
     12 1763  754 (S30W45 L= 83)   545     96        0 (1723/ 738-1798/ 782) -12.5 W40-W49 S31-S27
     13 1140  607 (S20E03 L= 34)   279     92        0 (1092/ 579-1179/ 626) -38.7 E06-E01 S22-S18
     14 1070  443 (S29E13 L= 25)   640    169        1 (1007/ 402-1140/ 500)  20.0 E17-E06 S30-S26
     15 1328  416 (S39W04 L= 42)  1162    209        4 (1257/ 348-1390/ 495)  36.0 E03-W09 S42-S34
     16  550  379 (S22E52 L=346)   373     94        0 ( 524/ 355- 588/ 394)   0.0 E55-E48 S23-S20
     17  619  323 (S27E50 L=348)   381     56        0 ( 595/ 299- 642/ 348)  78.7 E52-E48 S29-S24
     18 1130  291 (S43E16 L= 22)   408     56        1 (1109/ 270-1157/ 321) -90.0 E17-E14 S45-S39
  

    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