Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  492 1619 (N32E50 L=141)   402     57        1 ( 477/1592- 507/1646)  90.0 E52-E48 N29-N33
      2  772 1377 (N19E21 L=171)   616    249        2 ( 742/1263- 800/1490) -78.7 E23-E17 N12-N27
      3  352 1458 (N19E56 L=136)   448    104        2 ( 318/1438- 386/1483)  26.6 E59-E53 N17-N21
      4 1654 1408 (N33W43 L=234)   301     63        2 (1631/1380-1674/1430) -33.7 W40-W44 N30-N34
      5 1731 1380 (N32W49 L=241)   502     53        1 (1708/1352-1752/1407) -26.6 W47-W51 N29-N33
      6 1805 1095 (N14W55 L=247)  1077    218        9 (1760/ 986-1847/1212) -52.4 W50-W60 N07-N20
      7  200  915 (S16E60 L=131)   481     94        1 ( 189/ 878- 215/ 962)  80.5 E61-E58 S18-S12
      8  319  850 (S18E47 L=144)   319     56        1 ( 298/ 823- 337/ 881) -76.0 E49-E45 S20-S16
      9  341  815 (S20E45 L=146)   594    201        2 ( 280/ 774- 375/ 865)  35.5 E51-E42 S23-S16
     10 1765  659 (S13W60 L=251)  5231    599       15 (1652/ 536-1847/ 808) -66.0 W50-W68 S21-S03
     11  499  736 (S23E32 L=159)   276     61        0 ( 472/ 718- 531/ 764) -11.3 E34-E29 S24-S21
     12 1138  675 (S19W12 L=203)   516    190        2 (1042/ 653-1217/ 706)  15.3 W05-W16 S21-S16
     13 1669  673 (S13W50 L=241)   583    123        2 (1622/ 657-1717/ 699)  12.5 W45-W54 S14-S11
     14 1419  537 (S25W34 L=225)   429    131        2 (1378/ 482-1465/ 570)  60.3 W31-W36 S29-S22
     15  823  479 (S37E07 L=184)   510    106        0 ( 768/ 456- 863/ 511)  29.7 E11-E04 S39-S34
     16 1348  295 (S44W41 L=232)   371     77        1 (1324/ 277-1379/ 320)  11.3 W38-W43 S45-S41
  

    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