Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  523 1672 (N18E54 L= 43)   952    189        1 ( 475/1582- 581/1738) -81.3 E56-E50 N11-N24
      2  831 1676 (N26E32 L= 65)   382     92        1 ( 781/1654- 868/1695)  51.3 E34-E29 N23-N28
      3  419 1627 (N14E60 L= 36)  1418    134        2 ( 376/1578- 469/1687) -82.9 E64-E56 N10-N18
      4 1336 1586 (N35W03 L= 99)   667    158        1 (1308/1511-1387/1646) -79.7 W01-W05 N29-N40
      5  950 1538 (N20E19 L= 78)   329    103        1 ( 931/1496- 985/1574) -51.3 E20-E16 N17-N22
      6  405 1353 (S03E47 L= 49)   377     78        1 ( 388/1313- 435/1378)  76.0 E48-E45 S05-S01
      7  156 1220 (S15E69 L= 27)  2933    333        6 (  79/1099- 214/1361)  72.6 E85-E63 S21-S05
      8  255 1285 (S10E59 L= 38)   422     86        2 ( 235/1244- 276/1322) -53.1 E61-E55 S11-S07
      9  226 1240 (S13E60 L= 36)   638     73        2 ( 203/1205- 262/1274) -45.0 E63-E56 S14-S10
     10 1756 1162 (N22W42 L=139)   425     86        0 (1738/1121-1776/1205)  71.6 W40-W44 N19-N25
     11  271 1127 (S18E53 L= 44)  1479    263        4 ( 165/1088- 375/1188)  21.8 E65-E43 S22-S13
     12 1398 1108 (N08W18 L=115)   224     55        0 (1369/1087-1419/1130)  45.0 W16-W19 N06-N09
     13 1635  975 (N08W37 L=133)   315     60        0 (1613/ 945-1654/1012) -36.9 W34-W38 N06-N09
     14  955  923 (S14E01 L= 95)   327     89        0 ( 935/ 877- 973/ 972) -59.0 E03-W00 S16-S11
     15 1920  893 (N14W68 L=165)  1214     91        3 (1897/ 846-1944/ 932) -38.7 W63-W73 N11-N16
     16 1709  543 (S12W62 L=158)   418    117        0 (1668/ 507-1742/ 587)  80.5 W59-W63 S14-S08
  

    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