Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  499 1664 (N58E45 L= 28)   603     74        1 ( 473/1637- 536/1691)  23.2 E51-E41 N55-N60
      2 1709 1459 (N19W61 L=134)   484    131        0 (1649/1423-1747/1487) -32.0 W55-W63 N16-N21
      3  779 1438 (N36E10 L= 63)   229     74        0 ( 743/1423- 814/1449)   0.0 E13-E07 N35-N36
      4  773 1185 (N20E13 L= 60)   379    119        1 ( 730/1149- 819/1212) -26.6 E16-E10 N16-N21
      5  388 1066 (N17E43 L= 31)   544    117        1 ( 354/1016- 415/1113) -60.3 E45-E40 N13-N20
      6 1360  598 (S26W15 L= 89)  1933    664        3 (1103/ 433-1556/ 784)  23.2 E05-W30 S33-S17
      7 1205  621 (S21W05 L= 79)   851    322        0 (1103/ 535-1307/ 739)  33.7 E02-W13 S25-S15
      8  706  584 (S16E28 L= 46)   479    124        1 ( 682/ 527- 740/ 635)  50.2 E30-E24 S19-S13
      9 1670  546 (S37W43 L=117)   546     81        0 (1632/ 523-1707/ 579)  15.9 W38-W47 S38-S35
     10 1443  458 (S38W21 L= 95)   643    197        0 (1376/ 406-1524/ 516)  21.0 W15-W28 S41-S35
     11  808  446 (S27E25 L= 48)   674    126        1 ( 781/ 402- 847/ 492) -74.1 E27-E22 S30-S23
  

    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