Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1761 1512 (N27W65 L= 77)   819    237        2 (1714/1486-1798/1529)  -5.7 W59-W69 N25-N27
      2 1531 1349 (N20W38 L= 50)   525    123        1 (1480/1316-1587/1383) -39.8 W34-W41 N17-N22
      3  429 1157 (N17E40 L=333)   635    115        4 ( 389/1121- 473/1195) -26.6 E42-E36 N14-N19
      4  240 1076 (N12E57 L=315)   531    192        0 ( 203/ 999- 268/1148) -52.1 E61-E54 N07-N16
      5  461  897 (N01E38 L=334)   419    100        2 ( 433/ 846- 481/ 941)  51.3 E40-E36 S02-N03
      6  371  863 (S01E46 L=326)   338    144        0 ( 336/ 810- 387/ 923)  50.2 E49-E43 S04-N02
      7 1566  788 (S15W34 L= 46)   390    127        0 (1513/ 757-1619/ 815)  14.0 W30-W38 S16-S13
      8  625  538 (S23E32 L=340)   665    175        1 ( 595/ 472- 653/ 604)  48.4 E36-E28 S27-S18
  

    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