Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1029 1679 (N45W25 L= 75)   887    152        1 ( 966/1643-1105/1703) -14.0 W18-W30 N42-N47
      2 1659 1538 (N16W61 L=111)   615    120        3 (1644/1486-1677/1604)  33.7 W57-W67 N12-N18
      3 1366 1545 (N26W39 L= 88)   523    108        2 (1313/1521-1414/1563) -14.0 W34-W42 N24-N27
      4 1561 1441 (N14W47 L= 96)   553     96        0 (1502/1423-1609/1459) -14.0 W41-W50 N12-N15
      5  915 1263 (N23W00 L= 50)   336     87        2 ( 876/1231- 942/1296)  18.4 E02-W01 N21-N25
      6 1634 1223 (N00W43 L= 92)   408    141        0 (1608/1162-1664/1292)  81.9 W41-W44 S03-N04
      7 1833 1143 (S11W58 L=108)  1326    208        3 (1775/1065-1901/1201)   5.7 W49-W70 S13-S08
      8  628 1057 (N18E23 L= 27)   585    186        1 ( 563/1004- 699/1124)  18.4 E28-E16 N16-N20
      9  170  986 (N23E66 L=343)  1906    172        1 ( 143/ 906- 217/1068)  78.7 E71-E60 N18-N28
     10  324  989 (N21E48 L=  2)   979    271        1 ( 232/ 899- 380/1049) -34.7 E57-E43 N15-N25
     11  962 1004 (N07E04 L= 46)   242     63        1 ( 933/ 968- 984/1030)  45.0 E05-E02 N05-N09
     12  258  869 (N15E57 L=353)   905    145        3 ( 229/ 812- 289/ 947) -76.0 E60-E53 N11-N19
     13 1492  907 (S12W24 L= 74)   296    107        3 (1466/ 865-1529/ 944)  53.1 W22-W26 S14-S10
     14 1740  857 (S24W42 L= 92)   914    189        1 (1684/ 772-1786/ 932)  35.0 W36-W47 S27-S20
     15 1795  804 (S29W48 L= 97)   363     55        1 (1772/ 780-1819/ 828)  11.3 W45-W50 S30-S27
     16 1404  748 (S19W15 L= 65)   760    297        1 (1284/ 730-1538/ 781) -19.4 W07-W24 S22-S15
     17  841  472 (S21E27 L= 23)   322     76        2 ( 820/ 444- 873/ 509) -78.7 E28-E25 S22-S17
     18  497  414 (S17E55 L=355)  1268    242        4 ( 441/ 351- 546/ 506)  23.2 E66-E45 S20-S11
  

    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