Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  775 1810 (N63W11 L= 41)   510     40        0 ( 749/1781- 794/1837)  35.5 W07-W15 N59-N65
      2  608 1735 (N63E13 L= 17)  2037    277        6 ( 489/1668- 724/1794)   3.4 E30-W04 N58-N65
      3  449 1440 (N43E31 L=359)   526    150        0 ( 372/1421- 518/1457) -11.3 E39-E24 N40-N44
      4 1515 1299 (N07W36 L= 65)   448     78        2 (1488/1259-1540/1337) -90.0 W33-W37 N04-N09
      5  691 1260 (N25E14 L= 16)   515    132        2 ( 645/1218- 741/1300)   8.1 E17-E10 N23-N27
      6 1828 1266 (S04W60 L= 90)   354     48        0 (1806/1244-1848/1283) -36.9 W57-W62 S05-S02
      7 1833 1125 (S11W56 L= 85)   687    154        0 (1803/1055-1861/1198)  33.7 W51-W60 S14-S08
      8 1418 1104 (S01W24 L= 54)   236     24        1 (1396/1081-1436/1129)  45.0 W22-W25 S02-N00
      9  946  652 (S14E14 L= 16)   280    101        0 ( 926/ 612- 967/ 700)  63.4 E14-E11 S17-S11
     10  604  481 (S17E41 L=349)  1018    293        2 ( 549/ 367- 644/ 595)  35.2 E51-E33 S22-S10
     11 1643  515 (S44W32 L= 61)   578    124        1 (1582/ 499-1703/ 533) -13.0 W25-W38 S46-S41
     12 1373  420 (S42W06 L= 36)  1131    257        2 (1261/ 369-1489/ 448)  -6.0 E03-W16 S44-S38
  

    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