Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1290 1511 (N32W02 L=124)   383    104        0 (1257/1460-1317/1549) -38.7 W00-W05 N28-N33
      2 1702 1410 (N40W38 L=159)   430     68        2 (1673/1389-1736/1437)  -9.5 W34-W42 N37-N41
      3 1627 1277 (N27W32 L=153)  1339    241        3 (1554/1177-1701/1333) -19.7 W24-W39 N23-N30
      4  548 1293 (S01E36 L= 86)   733    138        1 ( 475/1243- 591/1330) -12.5 E41-E32 S03-N00
      5 1446 1225 (N18W19 L=141)   723    258        2 (1348/1148-1515/1325)  65.0 W14-W22 N11-N26
      6  327 1143 (S15E49 L= 73)   317    143        0 ( 282/1103- 389/1167)  29.1 E52-E43 S18-S12
      7 1802 1025 (N18W52 L=173)  1482    193        4 (1746/ 964-1860/1089) -35.5 W45-W58 N14-N22
      8  172  960 (S29E65 L= 57)  2186    211        2 ( 133/ 874- 209/1056)  40.6 E74-E59 S34-S22
      9 1648  801 (N01W45 L=166)   908    211        1 (1621/ 712-1679/ 903) -74.7 W42-W47 S04-N07
     10 1399  814 (S06W27 L=149)  1027    238        4 (1303/ 773-1483/ 860)  45.0 W22-W32 S10-S01
     11 1786  833 (N07W56 L=177)   371     63        0 (1751/ 819-1818/ 844)  15.9 W52-W59 N05-N08
     12  183  805 (S39E65 L= 56)   647     52        1 ( 156/ 785- 207/ 826)  21.8 E72-E60 S40-S36
     13 1232  772 (S13W19 L=140)   532    182        2 (1146/ 726-1297/ 796)   5.7 W13-W23 S14-S11
     14  698  701 (S32E11 L=110)   927    261        2 ( 617/ 602- 797/ 785)  66.4 E15-E06 S40-S24
     15 1712  564 (S10W62 L=184)   630     90        0 (1678/ 540-1756/ 594)  45.0 W59-W65 S11-S06
     16 1429  314 (S32W58 L=179)  2052    333        2 (1289/ 272-1562/ 341)  12.3 W44-W78 S37-S28
  

    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