Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1041 1875 (N57W48 L=316)   919     85        2 (1001/1839-1073/1915) -90.0 W42-W59 N53-N59
      2 1709 1568 (N14W68 L=336)   506     60        0 (1689/1532-1724/1594)  36.9 W65-W71 N12-N15
      3 1833  816 (S31W52 L=320)   479     59        0 (1820/ 782-1848/ 849)  76.0 W50-W54 S32-S28
      4 1558  753 (S25W25 L=293)   882    263        1 (1452/ 652-1654/ 823)   9.5 W15-W33 S28-S22
      5 1282  676 (S21W05 L=274)   319    125        0 (1255/ 627-1330/ 708)  33.7 W02-W08 S23-S19
      6  411  673 (N00E48 L=220)   386     56        0 ( 383/ 651- 441/ 693) -14.0 E50-E45 N00-N01
      7 1345  608 (S27W08 L=276)   310     97        0 (1323/ 563-1368/ 658)  45.0 W05-W10 S29-S24
      8 1422  550 (S33W11 L=280)   465    118        0 (1365/ 535-1474/ 574) -14.0 W07-W15 S35-S31
      9 1691  496 (S49W37 L=306)   935    121        1 (1648/ 461-1717/ 554) -51.3 W31-W43 S53-S43
     10 1067  488 (S27E13 L=255)   553    123        1 (1041/ 431-1099/ 541)  40.6 E16-E09 S29-S23
     11  492  346 (S21E61 L=208)   496     61        0 ( 470/ 318- 515/ 374)  63.4 E62-E59 S23-S19
  

    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