Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1318 1537 (N35W05 L= 76)  1394    270        3 (1240/1417-1369/1646) -47.3 E02-W10 N28-N41
      2  400 1489 (N10E55 L= 16)   866    151        2 ( 374/1434- 428/1551) -66.8 E57-E50 N06-N13
      3  704 1451 (N14E30 L= 41)   269     75        0 ( 684/1415- 734/1484)  78.7 E31-E29 N11-N16
      4 1628 1165 (N20W35 L=106)  1341    251        3 (1567/1059-1689/1266)  90.0 W31-W39 N13-N27
      5 1844 1033 (N18W58 L=129)  1545    167        1 (1807/ 960-1873/1118) -48.8 W52-W62 N14-N22
      6 1327 1017 (N04W17 L= 89)   438     99        1 (1289/ 983-1367/1048)  51.3 W15-W19 N01-N06
      7  196 1009 (S23E61 L= 10)   428     60        2 ( 180/ 977- 218/1042)  90.0 E63-E58 S24-S20
      8  549  700 (S34E24 L= 47)  1029    210        1 ( 466/ 672- 650/ 744)   7.1 E31-E15 S35-S31
      9  911  706 (S24W01 L= 72)   282     68        0 ( 875/ 687- 942/ 727)  56.3 W00-W03 S25-S22
     10 1616  558 (S14W53 L=124)   618    122        0 (1576/ 512-1661/ 599)  80.5 W50-W55 S17-S11
     11 1676  546 (S13W59 L=130)   399     66        1 (1648/ 522-1698/ 578)  78.7 W57-W61 S15-S10
     12 1569  397 (S25W61 L=132)   649    123        2 (1518/ 369-1609/ 440)  56.3 W57-W65 S27-S21
  

    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