Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  342 1111 (S15E46 L=207)  1202    229        2 ( 280/1022- 406/1204)  65.2 E50-E42 S21-S08
      2  242 1139 (S16E56 L=197)   427     60        0 ( 228/1101- 262/1170)  90.0 E57-E54 S18-S13
      3 1856 1059 (N21W56 L=310)   472    115        0 (1822/1037-1877/1098)  59.0 W52-W58 N19-N24
      4  424  988 (S21E37 L=217)   291     54        1 ( 396/ 973- 451/1003)   0.0 E39-E34 S21-S19
      5  175  819 (S37E59 L=195)  1417    307        1 ( 129/ 784- 211/ 866)  35.0 E64-E54 S39-S32
      6  390  787 (S34E37 L=216)   792    159        1 ( 322/ 747- 452/ 830)  -9.5 E43-E31 S35-S32
      7  787  771 (S25E07 L=246)  1763    452        1 ( 547/ 693- 977/ 822)   9.5 E24-W06 S28-S21
      8 1351  675 (S15W29 L=282)  1860    476        3 (1141/ 633-1566/ 731)  17.8 W14-W43 S18-S07
      9  546  549 (S46E21 L=232)   754    122        0 ( 510/ 510- 590/ 595) -14.0 E24-E16 S48-S42
     10  674  528 (S43E10 L=244)   304     79        0 ( 638/ 507- 710/ 551)  45.0 E12-E07 S45-S40
     11 1068  485 (S34W19 L=273)   230     63        0 (1037/ 461-1091/ 506)  -9.5 W16-W22 S35-S33
     12 1403  288 (S36W54 L=307)   573     80        1 (1361/ 271-1449/ 310)  51.3 W50-W55 S37-S32
  

    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