Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  890 1661 (N27E28 L=313)   630    162        0 ( 868/1584- 916/1726) -55.0 E31-E24 N21-N31
      2 1332 1637 (N39W00 L=341)   767    181        0 (1301/1535-1354/1704) -63.4 E03-W04 N32-N44
      3 1455 1350 (N24W16 L=357)   298    115        0 (1407/1311-1498/1386)  -8.1 W12-W19 N22-N24
      4 1650 1291 (N27W31 L= 12)   362     93        0 (1636/1251-1665/1347) -63.4 W29-W33 N24-N30
      5  236 1243 (S13E60 L=281)  1715    234        2 ( 194/1147- 261/1332) -61.4 E67-E55 S18-S07
      6 1724 1246 (N27W38 L= 19)   330     91        0 (1707/1194-1734/1292) -59.0 W36-W39 N24-N29
      7  309 1218 (S13E51 L=289)   388     94        1 ( 275/1185- 335/1259)  59.0 E54-E49 S14-S09
      8  443 1190 (S12E39 L=301)   374    109        1 ( 408/1153- 481/1236) -20.6 E43-E35 S13-S09
      9 1402 1187 (N13W17 L=357)   490     91        2 (1360/1163-1440/1218)  -9.5 W13-W19 N10-N14
     10  393 1090 (S18E41 L=300)   598     99        3 ( 341/1063- 438/1124)  20.6 E45-E37 S20-S15
     11 1295 1102 (N05W13 L=353)   404     93        0 (1247/1084-1352/1117)  33.7 W09-W15 N03-N07
     12  889  934 (S15E05 L=336)   261     54        2 ( 867/ 902- 916/ 962) -45.0 E06-E03 S16-S13
     13  801  841 (S23E08 L=333)   742    278        2 ( 771/ 709- 861/ 954) -79.4 E10-E04 S30-S14
     14 1777  880 (N08W51 L= 32)   524    120        2 (1756/ 839-1816/ 933) -45.0 W48-W55 N06-N11
     15 1738  795 (N03W51 L= 32)   466     81        1 (1706/ 770-1774/ 823)  -8.1 W47-W54 N01-N03
     16  261  667 (S48E55 L=286)   585     96        0 ( 217/ 644- 309/ 692)   0.0 E62-E47 S48-S46
     17 1289  538 (S25W31 L= 12)   243     37        0 (1266/ 514-1311/ 563) -14.0 W29-W33 S26-S23
     18 1543  487 (S20W51 L= 32)   979    218        4 (1505/ 413-1582/ 550) -63.4 W45-W57 S23-S15
     19 1082  461 (S36W23 L=  4)   311     87        0 (1043/ 443-1131/ 479)  45.0 W20-W25 S38-S33
  

    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