Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1509 1722 (N31W60 L=232)  1256    197        2 (1402/1697-1584/1746) -19.4 W50-W68 N27-N34
      2 1136 1668 (N39W27 L=199)  1008    232        0 (1022/1621-1235/1705)  -3.0 W17-W36 N37-N39
      3 1342 1450 (N20W30 L=202)   291     90        0 (1300/1427-1375/1476) -59.0 W28-W32 N17-N22
      4  371 1288 (N33E40 L=132)   918    170        1 ( 291/1271- 456/1309) -19.4 E48-E31 N30-N36
      5  909 1189 (N15E03 L=169)   411    146        0 ( 836/1161- 976/1210) -39.8 E06-W00 N12-N17
      6  537 1161 (N22E27 L=145)   394    135        0 ( 494/1108- 590/1184)  15.9 E30-E23 N19-N22
      7 1668 1020 (S12W40 L=212)   703    276        1 (1592/ 947-1713/1143)  36.9 W33-W45 S15-S05
      8 1781 1048 (S13W50 L=222)   942    221        3 (1720/ 956-1847/1128)  25.0 W43-W58 S17-S10
      9 1548 1054 (S07W32 L=204)   264     60        1 (1511/1033-1573/1074) -14.0 W29-W33 S08-S05
     10 1791 1011 (S16W51 L=223)   490     73        0 (1778/ 971-1804/1051)  68.2 W48-W52 S18-S13
     11  180  971 (N17E62 L=110)  2297    242        2 (  85/ 916- 252/1045)   6.1 E82-E52 N14-N21
     12 1845  892 (S25W56 L=228)   896     91        2 (1807/ 858-1879/ 926)  14.0 W51-W60 S26-S22
     13 1277  866 (S11W11 L=183)   205      3        0 (1259/ 848-1294/ 886) -63.4 W09-W11 S12-S10
     14 1414  738 (S22W17 L=189)   283     65        0 (1399/ 702-1430/ 774)  53.1 W15-W18 S24-S20
     15 1218  632 (S24W02 L=174)   796    260        2 (1162/ 529-1279/ 745)  47.7 E03-W08 S29-S18
     16 1024  615 (S20E10 L=162)   838    262        4 ( 979/ 498-1080/ 730)  45.0 E15-E03 S26-S14
     17  705  484 (S21E34 L=138)  1326    400        0 ( 650/ 321- 768/ 657)  39.6 E47-E23 S30-S11
     18 1275  543 (S31W04 L=176)   406    117        0 (1220/ 514-1325/ 580)   0.0 W00-W08 S32-S29
     19 1526  407 (S50W23 L=195)  1028    179        0 (1438/ 377-1593/ 446)  -3.4 W13-W30 S52-S47
  

    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