Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  504 1240 (N31E28 L= 37)  1865    364        3 ( 361/1186- 672/1276) -17.7 E40-E15 N24-N34
      2  431 1021 (N20E37 L= 27)  1302    337        1 ( 303/ 958- 570/1096) -32.3 E47-E27 N13-N26
      3  192  951 (N20E62 L=  3)  1233    186        2 ( 148/ 871- 252/1026) -50.7 E67-E55 N14-N25
      4 1831  903 (S25W52 L=116)   388     50        1 (1805/ 878-1851/ 932)   9.5 W48-W54 S26-S24
      5 1494  835 (S18W22 L= 86)   381    153        0 (1477/ 764-1509/ 908)  54.5 W19-W24 S21-S14
      6 1692  572 (S43W36 L=100)  1197    227        3 (1592/ 516-1787/ 629)  -2.5 W24-W48 S44-S40
      7  391  442 (S13E63 L=  2)  2537    279        8 ( 315/ 379- 475/ 511) -68.2 E72-E55 S18-S07
      8  513  381 (S19E55 L=  9)   357     61        1 ( 485/ 363- 536/ 402)   8.1 E58-E51 S20-S17
      9  934  287 (S37E31 L= 33)  1966    229        5 ( 844/ 204-1010/ 368)  19.2 E44-E21 S41-S32
     10  523  287 (S26E64 L=  0)  1227     84        2 ( 480/ 255- 549/ 318) -78.7 E69-E58 S28-S23
     11  784  198 (S39E51 L= 13)   414     75        0 ( 769/ 177- 812/ 235)  39.8 E54-E46 S41-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