Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1759 (N30E52 L= 31)  1197    269        3 ( 611/1688- 672/1837) -53.1 E57-E46 N24-N36
      2  631 1644 (N22E44 L= 39)   587    144        1 ( 590/1615- 693/1670)  -6.3 E47-E38 N19-N23
      3  391 1610 (N14E63 L= 20)   834    111        3 ( 361/1563- 430/1659) -81.9 E66-E60 N11-N18
      4  469 1526 (N11E50 L= 33)   453    151        2 ( 432/1455- 499/1588) -71.6 E53-E48 N05-N14
      5  278 1316 (S05E58 L= 25)   597    144        4 ( 248/1268- 309/1369) -35.5 E62-E53 S07-S02
      6 1367 1229 (N16W14 L= 97)   855    275        2 (1293/1110-1399/1349) -45.0 W06-W18 N09-N22
      7  416 1189 (S10E42 L= 41)   421    145        1 ( 353/1145- 475/1221)  45.0 E45-E37 S13-S06
      8  735  990 (S14E16 L= 67)   232     81        0 ( 714/ 940- 762/1025)  90.0 E16-E14 S16-S11
      9 1484  841 (S02W32 L=115)   429    170        2 (1450/ 758-1503/ 926) -77.5 W29-W33 S07-N02
     10  564  881 (S24E24 L= 59)   293    143        0 ( 502/ 853- 624/ 898)  36.9 E28-E20 S27-S21
     11 1814  853 (N09W57 L=141)   373    123        0 (1777/ 827-1849/ 878)  39.8 W54-W60 N06-N11
     12  215  792 (S38E54 L= 29)   742    192        3 ( 172/ 764- 257/ 816)  -6.3 E59-E49 S39-S37
     13 1119  624 (S24W17 L=100)   487    178        4 (1033/ 603-1207/ 643)  20.0 W11-W22 S27-S21
     14 1401  569 (S19W37 L=120)   346    171        2 (1318/ 549-1479/ 595)  33.7 W31-W40 S22-S15
     15  340  540 (S54E44 L= 39)   409     60        0 ( 311/ 521- 366/ 561)  -9.5 E47-E40 S54-S52
     16 1490  363 (S28W54 L=138)  3692    707        8 (1238/ 320-1685/ 429)  31.5 W37-W69 S36-S17
  

    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