Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1161 1588 (N33W26 L=147)   556    155        1 (1099/1535-1215/1612) -60.3 W23-W28 N28-N36
      2  699 1577 (N47E06 L=115)   242     57        1 ( 669/1563- 736/1589) -18.4 E08-E02 N45-N47
      3 1376 1514 (N22W36 L=157)   296    111        0 (1340/1475-1425/1550)   6.3 W32-W41 N20-N23
      4  647 1201 (N24E18 L=103)   247    108        1 ( 627/1144- 671/1248) -90.0 E19-E17 N20-N26
      5  205 1100 (N28E58 L= 63)  1591    236        6 ( 137/1003- 273/1177)  29.4 E68-E49 N22-N31
      6 1858  973 (S22W57 L=177)   916    115        2 (1839/ 922-1885/1027) -69.4 W54-W60 S25-S17
      7  784  677 (S09E23 L= 98)   255     54        0 ( 765/ 645- 806/ 699)  26.6 E24-E20 S10-S07
      8 1623  587 (S39W29 L=150)  1297    311        2 (1488/ 521-1730/ 674)   5.0 W17-W40 S40-S35
      9  852  403 (S27E29 L= 92)   459     78        1 ( 807/ 383- 896/ 427) -31.0 E31-E26 S29-S25
     10  632  347 (S26E48 L= 73)   582    133        0 ( 596/ 306- 693/ 378)   8.7 E53-E40 S27-S24
     11 1031  152 (S54E38 L= 82)   854     73        2 ( 995/ 122-1068/ 184)   0.0 E44-E31 S57-S51
  

    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