Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1603 1703 (N25W72 L= 40)   781     65        0 (1569/1677-1630/1719) -56.3 W68-W78 N23-N26
      2 1087 1656 (N41W27 L=355)   253     75        1 (1053/1638-1127/1677) -18.4 W24-W31 N40-N42
      3 1620 1603 (N20W62 L= 30)   388     40        3 (1601/1575-1640/1627)  18.4 W59-W65 N18-N20
      4 1440 1571 (N24W45 L= 13)   630    185        1 (1367/1506-1491/1609) -77.5 W41-W48 N19-N28
      5 1330 1485 (N23W33 L=  1)   800    268        0 (1219/1441-1425/1519) -33.7 W24-W37 N17-N26
      6 1068 1380 (N25W13 L=341)   605    191        2 (1017/1300-1111/1461) -90.0 W11-W14 N19-N31
      7  304 1432 (N50E52 L=276)   543     69        0 ( 280/1404- 332/1458)  31.0 E56-E47 N47-N51
      8  535 1313 (N36E25 L=303)  1128    186        4 ( 455/1266- 635/1354) -29.7 E31-E17 N31-N40
      9 1209 1235 (N13W16 L=344)   510    185        2 (1162/1172-1279/1292) -73.3 W14-W18 N07-N17
     10  728 1207 (N24E13 L=315)   269     97        0 ( 693/1176- 768/1243) -56.3 E14-E10 N21-N27
     11 1624  983 (S13W35 L=  3)  1075    277        0 (1562/ 869-1676/1107)  35.5 W27-W41 S17-S07
     12  865 1013 (N10E09 L=319)   767    247        1 ( 794/ 928- 955/1103) -77.0 E11-E05 N03-N16
     13  278  962 (N20E52 L=276)  1927    329        6 ( 213/ 862- 339/1084)  77.0 E59-E46 N14-N27
     14  573  943 (N13E29 L=299)   298     70        1 ( 544/ 921- 605/ 967) -45.0 E30-E27 N11-N14
     15 1433  724 (S22W16 L=344)   316    140        0 (1372/ 698-1488/ 765)  -6.3 W11-W20 S22-S20
     16  263  728 (N07E60 L=268)   408     71        0 ( 227/ 714- 298/ 744)   0.0 E65-E56 N06-N07
     17 1631  706 (S30W30 L=358)   862    113        2 (1594/ 655-1679/ 743)  14.0 W26-W34 S32-S26
     18  835  333 (S30E35 L=293)   598    102        1 ( 799/ 306- 878/ 367)  12.5 E38-E29 S31-S27
     19 1084  300 (S40E21 L=307)   582    128        1 (1053/ 246-1106/ 363)  48.4 E26-E17 S43-S34
  

    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