Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  663 1703 (N61E06 L=245)   446     85        0 ( 644/1663- 695/1736)  45.0 E09-E01 N57-N63
      2 1279 1648 (N35W39 L=290)  1407    251        4 (1171/1588-1390/1690) -58.4 W33-W43 N28-N41
      3  998 1407 (N29W09 L=260)   259     92        0 ( 960/1393-1047/1417) -18.4 W06-W12 N28-N30
      4  322 1155 (N31E46 L=205)   846    162        2 ( 246/1125- 407/1183) -19.4 E55-E38 N27-N33
      5 1777  853 (S25W46 L=297)  2809    457        3 (1660/ 709-1890/1010)  19.7 W33-W62 S31-S19
      6  931  588 (S16E17 L=234)   217     61        0 ( 905/ 572- 965/ 611) -56.3 E18-E15 S17-S13
      7  905  438 (S25E24 L=227)   650    253        0 ( 848/ 400- 960/ 493) -15.9 E27-E20 S27-S20
      8 1560  437 (S47W22 L=273)   483    106        0 (1525/ 396-1604/ 470)  12.5 W18-W27 S49-S44
      9  534  349 (S22E56 L=195)  2863    218        6 ( 453/ 297- 635/ 411) -45.0 E64-E49 S26-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