Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  844 1594 (N35E10 L=324)   374    105        1 ( 821/1542- 860/1636)  81.9 E11-E08 N30-N37
      2 1618 1525 (N27W45 L= 20)   430     54        1 (1592/1494-1645/1560) -63.4 W43-W47 N24-N29
      3 1652 1295 (N12W43 L= 17)   321     58        1 (1633/1259-1667/1331) -90.0 W41-W44 N09-N13
      4 1272 1153 (N04W15 L=350)  1692    439        5 (1108/ 973-1347/1274) -54.8 W06-W19 S06-N12
      5  989 1169 (N06E02 L=333)   209     86        0 ( 937/1156-1032/1182)   0.0 E04-W01 N05-N07
      6  842 1111 (N04E10 L=324)   644    227        2 ( 806/1049- 917/1162) -50.2 E12-E05 N00-N06
      7  287  791 (S12E52 L=282)  2246    452        4 ( 225/ 646- 368/ 942)  43.4 E62-E43 S20-S03
      8  412  914 (S06E40 L=294)   234     75        0 ( 383/ 886- 438/ 931) -45.0 E42-E38 S07-S04
      9 1796  825 (S17W54 L= 28)   499    111        1 (1763/ 774-1833/ 861)  38.7 W50-W57 S19-S14
     10  957  611 (S27E07 L=328)   231     68        0 ( 932/ 583- 976/ 641)  45.0 E08-E04 S29-S25
     11 1523  436 (S42W38 L= 12)  1800    343        1 (1376/ 379-1681/ 522)  18.4 W25-W50 S46-S36
  

    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