Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1811 (N46E43 L=262)  1620    141        3 ( 626/1780- 756/1845) -14.0 E53-E33 N42-N49
      2 1173 1701 (N41W04 L=309)   613    145        2 (1093/1684-1240/1717)   0.0 E02-W09 N39-N42
      3  618 1658 (N32E38 L=268)   362    112        0 ( 571/1618- 651/1675) -18.4 E42-E33 N29-N33
      4 1532 1643 (N42W35 L=340)   518    110        0 (1481/1608-1573/1667) -23.2 W30-W39 N39-N43
      5  790 1600 (N29E22 L=283)   451    106        2 ( 763/1556- 819/1653) -66.8 E24-E19 N25-N32
      6 1757 1564 (N40W61 L=  7)   663     81        1 (1714/1538-1790/1584) -12.5 W55-W67 N38-N41
      7  316 1473 (N18E57 L=248)  1333    172        3 ( 277/1408- 374/1541) -25.3 E67-E48 N13-N22
      8  568 1477 (N19E35 L=270)   250     40        2 ( 548/1451- 586/1500) -56.3 E36-E33 N17-N20
      9  927 1423 (N18E10 L=296)   302    122        0 ( 888/1376- 980/1475)  63.4 E11-E07 N14-N22
     10 1122 1328 (N14W03 L=309)   307    172        0 (1078/1260-1162/1409) -58.0 W00-W06 N10-N18
     11 1562 1313 (N18W32 L=338)   419    182        0 (1493/1279-1646/1358)  -9.5 W27-W39 N16-N20
     12 1624 1249 (N15W37 L=343)   277     54        1 (1595/1229-1648/1273)   0.0 W34-W39 N13-N16
     13 1495 1180 (N09W28 L=333)   282    122        0 (1465/1123-1534/1225) -56.3 W25-W30 N05-N11
     14 1503 1048 (N01W29 L=335)   428    102        0 (1461/1022-1552/1071)  18.4 W26-W32 N00-N03
     15 1794  935 (S02W54 L=359)  1377    215        5 (1708/ 870-1886/ 979)  22.4 W46-W64 S06-N02
     16  448  629 (S34E40 L=265)   351     60        0 ( 432/ 589- 466/ 663)  45.0 E42-E37 S36-S31
     17 1656  639 (S21W49 L=355)   368    107        0 (1614/ 615-1713/ 659)  26.6 W46-W54 S22-S18
     18  866  387 (S48E06 L=300)   492     70        1 ( 824/ 353- 902/ 413)  15.9 E09-E02 S51-S45
  

    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