Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1833 (N57E17 L=354)  1394    151        3 (1172/1790-1317/1868)   3.4 E25-E07 N54-N59
      2  602 1709 (N26E53 L=318)  2930    341        7 ( 565/1565- 653/1830) -25.9 E75-E39 N16-N34
      3 1655 1629 (N57W39 L= 49)   832     69        2 (1634/1590-1675/1665)  40.2 W33-W49 N52-N63
      4  496  857 (S27E29 L=341)   338    108        0 ( 483/ 799- 507/ 905) -63.4 E31-E27 S30-S24
      5  868  784 (S22E03 L=  8)   193     83        1 ( 855/ 749- 880/ 825) -63.4 E04-E01 S23-S19
      6  649  742 (S30E16 L=355)   262    119        0 ( 592/ 726- 700/ 765)  35.5 E19-E12 S32-S27
      7  784  666 (S31E05 L=  6)   399     61        0 ( 751/ 641- 824/ 690)  59.0 E06-E02 S33-S28
      8  619  329 (S63E10 L=  1)   501    101        0 ( 576/ 311- 662/ 346)  22.6 E16-E04 S66-S59
      9 1395  294 (S35W57 L= 68)  1202    143        1 (1329/ 272-1456/ 328)  11.3 W48-W63 S37-S32
  

    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