Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  502 1383 (N12E41 L= 87)   271    104        1 ( 485/1347- 519/1429) -80.5 E43-E39 N09-N15
      2 1045 1350 (N19E04 L=124)   310     66        1 (1014/1331-1086/1378) -11.3 E06-E01 N17-N20
      3 1646 1256 (N24W39 L=167)   373     40        2 (1620/1234-1677/1275)  11.3 W36-W41 N22-N25
      4 1893 1120 (N20W68 L=196)   550     46        1 (1872/1101-1917/1140) -15.9 W64-W72 N19-N21
      5 1234 1055 (N05W12 L=140)   986    287        2 (1146/ 966-1333/1133) -23.2 W05-W19 N00-N08
      6 1877  961 (N10W66 L=194)   729     86        2 (1851/ 922-1901/1002) -38.7 W61-W70 N08-N12
      7 1831  865 (N04W61 L=190)   737    114        3 (1780/ 842-1869/ 900)  24.0 W55-W66 N01-N06
      8  868  711 (S22E04 L=124)   591    271        2 ( 825/ 627- 963/ 810) -45.0 E07-W03 S26-S16
      9  623  628 (S32E21 L=107)   271     60        1 ( 608/ 597- 639/ 666) -68.2 E21-E19 S34-S29
     10 1675  552 (S18W58 L=186)   802    121        0 (1637/ 505-1711/ 603)  81.9 W56-W60 S21-S14
     11 1164  547 (S27W18 L=146)   394    101        0 (1137/ 502-1190/ 596)  90.0 W17-W19 S30-S23
     12  883  430 (S41W02 L=130)  1119    236        5 ( 786/ 401-1003/ 458)  20.6 E06-W10 S44-S36
  

    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