Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1528 1531 (N22W48 L=183)   888    294        3 (1494/1420-1548/1643)  55.3 W41-W53 N15-N28
      2 1780 1491 (N11W70 L=205)   745     96        0 (1749/1441-1806/1536)  80.5 W68-W71 N08-N14
      3 1695 1084 (S07W43 L=178)   730    240        2 (1631/ 967-1729/1167) -81.9 W40-W46 S14-N00
      4 1874 1076 (S14W61 L=195)  1339    284        3 (1838/1016-1910/1143)  20.0 W55-W67 S16-S10
      5 1603  670 (S29W29 L=164)  1133    186        3 (1497/ 571-1713/ 790)  19.3 W19-W39 S32-S24
      6  404  564 (S08E53 L= 82)   427     81        1 ( 387/ 525- 426/ 603)  80.5 E54-E51 S10-S04
      7  707  488 (S18E33 L=102)  1162    348        5 ( 634/ 415- 814/ 560) -63.4 E37-E29 S24-S12
      8  559  416 (S20E47 L= 88)   753    177        2 ( 497/ 351- 605/ 487) -84.8 E49-E44 S24-S13
      9  426  398 (S19E60 L= 74)  1517    192        4 ( 365/ 350- 483/ 464) -83.7 E64-E55 S22-S13
  

    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