Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  814 1867 (N68W21 L= 27)  1159    105        0 ( 792/1818- 849/1910)  19.2 W12-W35 N62-N73
      2  907 1872 (N64W31 L= 37)   477     59        0 ( 892/1839- 923/1905)  33.7 W26-W39 N59-N67
      3 1157 1748 (N44W35 L= 41)   756    133        2 (1124/1703-1189/1802)  28.6 W29-W40 N40-N47
      4  919 1741 (N51W17 L= 23)   489     77        1 ( 877/1721- 961/1757) -23.2 W13-W21 N49-N53
      5 1442 1635 (N28W47 L= 53)   727    130        0 (1392/1606-1505/1671) -49.4 W44-W51 N24-N31
      6 1341 1533 (N25W34 L= 40)   481    152        1 (1268/1510-1411/1547) -16.7 W29-W39 N22-N26
      7 1556 1444 (N14W45 L= 50)   891    204        1 (1516/1355-1589/1543)  57.5 W42-W50 N08-N19
      8 1612 1427 (N11W48 L= 54)   273     35        0 (1595/1406-1634/1445)  14.0 W46-W50 N10-N12
      9  910 1416 (N29W03 L=  9)   504    162        1 ( 835/1357- 963/1441) -58.0 E01-W05 N24-N32
     10 1171 1418 (N22W19 L= 25)   494    217        1 (1066/1402-1273/1432) -17.1 W12-W25 N19-N25
     11 1110 1358 (N20W13 L= 19)   274    105        0 (1062/1343-1163/1379) -38.7 W11-W16 N18-N22
     12  339 1119 (N24E43 L=322)   480     85        0 ( 306/1082- 379/1157) -50.2 E46-E40 N21-N27
     13  302  997 (N18E48 L=317)   478     70        2 ( 275/ 966- 334/1026) -53.1 E51-E45 N15-N19
     14  321  864 (N10E49 L=316)   368    112        2 ( 302/ 812- 341/ 926)  50.2 E52-E46 N07-N13
     15  204  881 (N13E61 L=304)   422     58        0 ( 188/ 844- 218/ 918)  45.0 E63-E58 N10-N14
     16  808  642 (S13E22 L=344)   205     77        1 ( 796/ 597- 826/ 684)  59.0 E23-E19 S15-S10
     17 1282  441 (S38W01 L=  7)  1626    366        2 (1159/ 326-1383/ 548)  23.2 E11-W10 S44-S32
     18 1454  270 (S61W14 L= 19)   767     71        1 (1412/ 246-1498/ 290) -14.0 W08-W20 S63-S58
     19 1319  186 (S65E07 L=359)   566     71        1 (1285/ 167-1350/ 208)  -6.3 E10-E01 S67-S61
  

    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