Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  707 1766 (N33E49 L=266)   404     70        4 ( 683/1729- 727/1794) -45.0 E52-E45 N30-N35
      2 1100 1670 (N38E15 L=300)  1131    227        1 (1040/1571-1149/1756) -42.7 E22-E08 N31-N43
      3 1583 1663 (N57W29 L=344)  1120    125        2 (1525/1629-1642/1688)  27.6 W20-W43 N50-N63
      4  830 1557 (N22E27 L=288)   955    343        3 ( 742/1401- 909/1679) -32.5 E39-E16 N14-N28
      5  280  944 (S26E50 L=265)  1354    302        2 ( 209/ 827- 386/1015)  36.4 E59-E40 S35-S20
      6  678  683 (S33E13 L=302)   681    221        1 ( 578/ 665- 793/ 699)  25.0 E20-E05 S35-S28
      7 1204  517 (S28W27 L=342)   759    143        2 (1160/ 451-1235/ 587) -66.0 W24-W28 S32-S23
      8 1515  422 (S24W54 L=  9)   363     46        1 (1495/ 395-1535/ 451) -53.1 W51-W56 S26-S22
      9 1340  402 (S31W42 L=357)   929    280        1 (1222/ 358-1441/ 439)  11.9 W33-W52 S33-S28
     10  556  360 (S62E20 L=295)  1862    247        5 ( 474/ 296- 665/ 438) -16.3 E34-E04 S66-S56
     11 1272  175 (S48W65 L= 20)   603     47        1 (1252/ 147-1284/ 206) -14.0 W58-W74 S49-S45
  

    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