Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  538 1630 (N58E20 L= 98)   427     54        1 ( 509/1612- 573/1651) -14.0 E24-E16 N55-N60
      2 1146 1613 (N35W28 L=147)  1040    186        1 (1055/1584-1240/1636) -32.5 W21-W34 N31-N38
      3  298 1325 (N42E48 L= 71)   506     82        1 ( 257/1293- 330/1357) -31.0 E54-E43 N39-N45
      4 1376 1029 (S04W20 L=138)   624    154        0 (1337/ 956-1417/1096) -90.0 W18-W20 S08-N00
      5 1700  966 (S17W40 L=158)  1166    324        1 (1593/ 860-1786/1063)  14.0 W29-W49 S19-S14
      6 1735  864 (S24W41 L=160)   297     86        2 (1716/ 821-1751/ 909)  51.3 W38-W43 S26-S21
      7  245  774 (N10E60 L= 59)  1122    205        0 ( 210/ 690- 302/ 856) -90.0 E63-E56 N04-N14
      8 1453  574 (S33W14 L=132)   856    167        2 (1413/ 510-1482/ 661)  85.2 W10-W16 S37-S25
      9 1416  475 (S39W09 L=127)   428     84        0 (1370/ 452-1443/ 501)   0.0 W05-W11 S40-S36
     10  636  429 (S19E42 L= 76)   419    100        0 ( 617/ 383- 653/ 485)  54.5 E46-E39 S21-S14
     11  564  327 (S24E55 L= 63)   402     60        0 ( 531/ 313- 598/ 343)  -8.1 E59-E52 S24-S22
  

    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