Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  547 1746 (N58E37 L=155)   603     45        0 ( 519/1722- 580/1767)   8.1 E41-E32 N55-N60
      2 1104 1689 (N41W18 L=210)  1264    192        1 (1007/1664-1201/1718) -10.6 W10-W26 N38-N43
      3 1281 1615 (N33W28 L=220)   826    160        0 (1229/1575-1362/1634)   4.8 W23-W35 N30-N34
      4 1645 1533 (N22W54 L=246)  1003    154        1 (1571/1497-1712/1568) -31.0 W48-W59 N18-N25
      5 1393  889 (S13W21 L=212)   242    119        0 (1336/ 875-1455/ 900)  -7.1 W16-W24 S13-S12
      6 1623  765 (S24W37 L=228)   442    131        0 (1583/ 706-1647/ 827)  50.2 W33-W38 S26-S20
      7 1012  675 (S20E06 L=186)   709    226        1 ( 971/ 582-1039/ 777)  81.3 E09-E03 S26-S13
      8 1425  604 (S31W21 L=213)   451     96        0 (1374/ 582-1467/ 635)  15.9 W17-W24 S32-S29
      9 1468  384 (S49W27 L=219)  1084    175        1 (1400/ 352-1543/ 420)   0.0 W20-W35 S52-S44
  

    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