Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  926 1877 (N70W32 L=184)   471     50        0 ( 909/1851- 951/1894)  10.6 W24-W40 N67-N72
      2 1516 1521 (N25W47 L=199)  1637    260        1 (1411/1433-1596/1605) -70.3 W43-W50 N18-N32
      3  392 1433 (N44E45 L=107)   386     77        1 ( 366/1406- 432/1455) -12.5 E48-E39 N41-N45
      4 1140 1297 (N21W13 L=164)   302    145        0 (1098/1241-1193/1348) -69.4 W11-W14 N16-N24
      5 1252 1073 (N05W14 L=166)   229     81        0 (1229/1039-1288/1107)  31.0 W12-W17 N03-N06
      6 1553  994 (S06W32 L=184)   410    108        1 (1497/ 967-1599/1022) -35.5 W28-W35 S08-S02
      7 1646  926 (S12W38 L=190)   264     60        1 (1631/ 893-1659/ 960)  56.3 W36-W39 S14-S10
      8  599  905 (N07E28 L=123)   257     44        0 ( 572/ 888- 625/ 924)  14.0 E30-E26 N06-N08
      9  270  700 (S02E61 L= 91)   874     94        2 ( 239/ 651- 298/ 736) -63.4 E63-E57 S04-N00
     10 1682  682 (S29W40 L=191)   925    130        1 (1628/ 640-1727/ 732)  16.7 W34-W44 S31-S26
     11  267  553 (S11E71 L= 81)   756     77        3 ( 239/ 519- 296/ 580) -63.4 E74-E67 S13-S09
     12  442  505 (S16E52 L= 99)   455     62        1 ( 428/ 468- 459/ 540)  59.0 E54-E49 S18-S13
  

    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