Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1390 1715 (N35W45 L=175)  1888    302        3 (1298/1654-1490/1760)  -6.7 W37-W54 N31-N36
      2  537 1601 (N50E25 L=105)   313     85        0 ( 519/1566- 552/1632) -71.6 E26-E22 N47-N53
      3 1676 1193 (S02W44 L=174)  1551    369        4 (1622/1041-1757/1324)  32.3 W36-W55 S09-N03
      4 1380 1231 (N07W25 L=155)   419    102        1 (1330/1178-1423/1285) -66.8 W23-W27 N03-N10
      5  112  939 (N15E70 L= 59)   814    121        2 (  87/ 913- 138/ 964) -26.6 E74-E66 N13-N16
      6 1648  811 (S24W35 L=165)   440    108        1 (1601/ 788-1680/ 851)  -9.5 W31-W37 S24-S21
      7 1234  517 (S32W02 L=131)   936    186        2 (1182/ 447-1260/ 608)  55.0 E03-W04 S36-S26
      8 1489  479 (S42W20 L=150)   406     62        1 (1454/ 452-1513/ 500) -11.3 W17-W22 S44-S40
      9  788  419 (S28E30 L= 99)   561    137        1 ( 761/ 367- 815/ 476)  36.9 E34-E26 S31-S24
     10  366  446 (S18E62 L= 68)   529    101        1 ( 349/ 417- 387/ 472)  26.6 E65-E58 S19-S16
  

    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