Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  463 1648 (N18E59 L=272)  1812    201        4 ( 420/1539- 519/1713) -61.7 E64-E53 N10-N23
      2 1262 1460 (N27W02 L=333)  1215    417        1 (1200/1300-1327/1623) -50.2 E06-W10 N18-N36
      3 1489 1385 (N29W19 L=350)   653    117        1 (1429/1348-1543/1418)   6.3 W14-W23 N27-N30
      4  866 1238 (N03E14 L=317)   789    195        2 ( 770/1189- 932/1302)  58.0 E18-E10 S01-N08
      5 1257 1183 (N10W09 L=340)   434    190        0 (1192/1139-1290/1243)  66.8 W05-W10 N06-N13
      6 1583 1122 (N16W31 L=  2)   721    169        1 (1537/1049-1635/1188) -36.9 W26-W35 N12-N18
      7  274 1127 (S16E53 L=278)   407     97        1 ( 247/1089- 293/1169) -45.0 E56-E50 S18-S14
      8  450 1055 (S17E36 L=295)   224     46        1 ( 434/1025- 470/1078) -71.6 E37-E34 S18-S15
      9 1416  945 (N01W24 L=355)   883    214        3 (1322/ 915-1510/ 992)   8.7 W17-W30 N00-N04
     10 1875  828 (N09W66 L= 37)  1890    144        5 (1824/ 756-1921/ 891)  66.0 W60-W73 N04-N13
     11 1094  806 (S15W10 L=341)   628    140        1 (1038/ 767-1143/ 861)  69.4 W07-W13 S18-S10
     12  405  840 (S30E36 L=295)   306    102        0 ( 370/ 806- 455/ 856)  35.5 E39-E32 S33-S28
     13 1517  672 (S10W40 L= 11)   363    143        0 (1482/ 621-1534/ 726) -71.6 W36-W41 S13-S07
     14 1749  645 (S05W60 L= 31)   997    126        2 (1721/ 607-1798/ 699) -10.3 W55-W66 S07-S02
     15  216  633 (S49E67 L=265)   871     61        1 ( 179/ 618- 251/ 651)   7.4 E80-E57 S50-S47
     16 1159  529 (S29W23 L=354)  1124    378        1 (1025/ 438-1245/ 628) -10.3 W11-W33 S32-S25
     17  917  440 (S42W12 L=343)   595    117        1 ( 864/ 417- 982/ 471)  11.3 W07-W17 S43-S40
     18  634  396 (S55E10 L=321)   332     36        1 ( 617/ 370- 658/ 415)  76.0 E11-E07 S57-S53
  

    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