Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  503 1666 (N40E56 L=169)  1105    137        1 ( 486/1593- 531/1727) -32.0 E66-E50 N34-N44
      2  855 1503 (N31E16 L=209)  1696    433        3 ( 744/1363-1002/1607) -26.6 E26-E04 N22-N38
      3  262 1379 (N17E64 L=161)  1163    172        3 ( 237/1318- 298/1462) -45.0 E73-E58 N13-N22
      4 1807 1358 (N28W64 L=289)   531     96        0 (1778/1313-1833/1391) -56.3 W60-W66 N24-N30
      5 1185 1300 (N20W08 L=234)   225     82        1 (1159/1269-1213/1338) -53.1 W06-W10 N18-N22
      6 1318 1298 (N21W17 L=243)   293     51        1 (1295/1272-1340/1322) -71.6 W15-W19 N19-N22
      7 1634 1177 (N15W41 L=266)   373    113        0 (1603/1137-1683/1216) -33.7 W38-W45 N13-N17
      8 1766 1128 (N13W53 L=278)   543     76        0 (1729/1107-1797/1150)  15.9 W49-W56 N11-N14
      9 1771  887 (S02W55 L=280)  1558    158        4 (1693/ 837-1840/ 925)  11.3 W47-W62 S05-N00
     10 1092  829 (S09W06 L=231)   228     78        0 (1070/ 793-1126/ 858)  53.1 W04-W07 S11-S07
     11  693  779 (S16E20 L=206)   615     93        1 ( 642/ 746- 743/ 821)  15.9 E23-E16 S18-S12
     12 1721  675 (S15W55 L=280)   628    117        1 (1690/ 637-1764/ 727) -26.6 W50-W61 S17-S12
     13 1135  649 (S21W10 L=235)   215     40        0 (1119/ 619-1149/ 672)  90.0 W08-W11 S22-S19
     14 1569  642 (S18W41 L=267)   304     21        1 (1546/ 617-1586/ 665) -76.0 W39-W43 S20-S16
     15  961  606 (S25E01 L=224)   372     90        0 ( 915/ 591-1009/ 620)   0.0 E04-W02 S25-S23
     16  858  450 (S37E08 L=218)   450    138        0 ( 792/ 433- 924/ 468)  16.7 E12-E02 S38-S35
  

    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