Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1521 1577 (N21W51 L=209)  1122    224        0 (1444/1524-1628/1608) -17.4 W42-W58 N16-N22
      2 1104 1551 (N33W22 L=180)   507    115        0 (1052/1523-1160/1572) -14.0 W18-W26 N31-N34
      3  646 1388 (N37E14 L=144)   265     70        0 ( 600/1375- 685/1399) -23.2 E17-E10 N35-N38
      4 1745 1330 (N01W55 L=213)   839    134        0 (1714/1252-1774/1391)  60.3 W52-W58 S03-N04
      5 1285  765 (S15W08 L=165)   454    216        2 (1210/ 693-1351/ 848)  22.6 W01-W13 S18-S12
      6  827  671 (S08E21 L=137)   958    292        0 ( 746/ 548- 893/ 780) -90.0 E22-E18 S17-N00
      7 1615  582 (S38W27 L=185)   305     30        0 (1593/ 564-1635/ 603)  14.0 W24-W28 S39-S36
      8 1314  549 (S30W04 L=161)   373    102        1 (1266/ 534-1366/ 566) -15.9 W00-W07 S31-S27
      9  974  345 (S33E25 L=132)   724    122        1 ( 930/ 290-1000/ 389)  71.6 E28-E22 S37-S28
  

    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