Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1734 (N39E16 L= 62)   633     77        2 (1065/1709-1159/1754)   7.1 E20-E12 N36-N40
      2 1300 1663 (N39E01 L= 77)  1310    183        4 (1233/1610-1363/1712) -15.3 E06-W05 N35-N43
      3 1532 1624 (N45W19 L= 96)   416     53        0 (1507/1603-1563/1649) -11.3 W16-W21 N43-N47
      4  919 1527 (N19E19 L= 58)   806    205        0 ( 862/1455- 989/1613) -30.3 E25-E13 N15-N23
      5  258 1480 (N04E68 L=  9)  2072    353       10 ( 202/1421- 303/1546) -53.1 E76-E61 N00-N08
      6 1381 1400 (N24W11 L= 89)  1073    332        4 (1334/1239-1424/1529) -55.5 W05-W17 N14-N31
      7 1789 1204 (N25W45 L=122)  1515    246        0 (1706/1120-1860/1290) -17.4 W36-W53 N22-N27
      8 1542 1031 (N07W30 L=107)   241     83        0 (1514/1001-1577/1053)  53.1 W28-W31 N04-N08
      9  262  925 (S30E50 L= 28)   675    145        2 ( 215/ 896- 312/ 951)   6.3 E54-E45 S31-S27
     10  742  446 (S50E02 L= 75)   683    136        1 ( 685/ 416- 817/ 486)   0.0 E07-W05 S51-S46
     11  861  211 (S65W26 L=103)   432     20        0 ( 836/ 196- 886/ 225)   0.0 W21-W29 S66-S62
  

    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