Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1070 1812 (N57E12 L= 95)  1318    312        0 ( 927/1775-1216/1851)  -1.7 E30-W04 N54-N57
      2  370 1628 (N30E69 L= 38)  1856    228        4 ( 310/1532- 452/1684)  71.6 E74-E62 N22-N34
      3 1046 1621 (N39E08 L= 99)  1126    318        1 ( 891/1584-1182/1680)  -2.4 E21-W03 N37-N42
      4  402 1498 (N22E53 L= 53)  1271    259        3 ( 360/1383- 443/1612) -68.2 E58-E48 N14-N29
      5 1613 1373 (N30W38 L=145)  2790    446        3 (1411/1309-1802/1445)  -4.8 W21-W58 N27-N32
      6 1791 1100 (N14W54 L=161)   543    124        0 (1774/1042-1825/1133) -80.5 W52-W57 N10-N16
      7  338 1066 (S05E47 L= 60)   539    116        0 ( 308/1010- 375/1118)  63.4 E48-E43 S09-S01
      8  599 1044 (S03E27 L= 80)   370    106        0 ( 568/1005- 617/1102)  71.6 E28-E25 S06-N00
      9  864 1031 (S01E10 L= 97)   556    116        0 ( 821/ 984- 903/1085)  63.4 E11-E07 S04-N02
     10  313  905 (S15E48 L= 59)   319     86        0 ( 297/ 863- 324/ 949) -78.7 E49-E46 S18-S12
     11 1160  716 (S16W12 L=119)   755    286        0 (1078/ 664-1216/ 778)  32.0 W07-W15 S20-S12
  

    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