Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1344 1756 (N56W21 L= 86)   773    122        2 (1277/1735-1403/1776)   3.8 W12-W28 N54-N58
      2  745 1541 (N30E27 L= 38)  1095    214        1 ( 703/1435- 811/1599)  85.2 E28-E23 N22-N34
      3  418 1528 (N25E53 L= 12)   695    115        0 ( 389/1473- 458/1563)  90.0 E55-E50 N21-N27
      4  301 1478 (N20E64 L=  1)   494     62        0 ( 269/1461- 335/1496)  -5.2 E69-E58 N19-N21
      5  539 1441 (N21E39 L= 26)   361     65        0 ( 509/1411- 572/1470)  63.4 E40-E36 N18-N23
      6  484 1154 (N03E36 L= 29)   799    204        0 ( 404/1084- 532/1197)  -5.2 E42-E31 N00-N05
      7 1743 1134 (N15W49 L=114)   346     57        0 (1728/1099-1757/1165) -76.0 W47-W50 N12-N16
      8 1559  953 (N02W35 L=100)   606     83        2 (1521/ 907-1584/ 995)  90.0 W31-W36 N00-N04
      9  934  908 (S07E04 L= 61)   417    170        0 ( 884/ 841- 987/ 968)  60.9 E06-E01 S11-S02
     10 1772  617 (S16W63 L=128)  2532    358        8 (1672/ 465-1846/ 748)  75.3 W59-W71 S27-S07
     11 1672  586 (S19W53 L=118)   626    118        3 (1623/ 554-1718/ 623)  39.8 W48-W56 S21-S16
     12 1410  305 (S42W46 L=111)   857     75        0 (1368/ 276-1441/ 345)  78.7 W41-W50 S44-S38
  

    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