Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1259 1592 (N30W32 L=234)   487    145        0 (1200/1569-1321/1610) -39.8 W27-W34 N27-N32
      2  871 1538 (N38W04 L=206)   856    242        0 ( 773/1497- 987/1579)   0.0 E03-W13 N35-N39
      3 1595 1425 (N12W47 L=249)   949    192        0 (1533/1341-1653/1498) -85.2 W44-W48 N05-N17
      4  157 1192 (N32E67 L=135)   725     96        0 ( 118/1164- 204/1219) -14.7 E77-E58 N29-N34
      5 1260 1015 (S02W13 L=215)   232     83        0 (1219/ 993-1292/1035) -36.9 W10-W14 S04-S01
      6 1528  807 (S21W26 L=228)   526    136        0 (1492/ 756-1553/ 874)  45.0 W22-W28 S23-S17
      7 1792  625 (S42W52 L=254)  1891    310        1 (1689/ 541-1874/ 757)  22.2 W39-W62 S45-S34
      8  827  513 (S22E25 L=177)   386    133        0 ( 794/ 448- 860/ 562)  36.9 E29-E21 S25-S19
      9  638  454 (S21E40 L=162)   913    203        1 ( 608/ 369- 671/ 544)  54.0 E44-E35 S27-S15
     10  390  508 (S13E57 L=144)   349     54        0 ( 375/ 479- 411/ 540)  23.2 E60-E53 S14-S11
     11 1356  338 (S50W05 L=206)  1647    331        3 (1210/ 284-1502/ 404)  -2.0 E09-W19 S51-S47
  

    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