Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1764 (N42E16 L=173)   817    194        1 (1088/1701-1184/1829) -28.3 E23-E10 N38-N47
      2  578 1578 (N14E43 L=146)   667    159        3 ( 561/1518- 603/1638) -53.1 E47-E39 N10-N18
      3 1314 1084 (N04W15 L=204)   492    144        1 (1245/1056-1400/1106)  18.4 W10-W19 N00-N06
      4 1797 1032 (N16W49 L=238)   362     71        0 (1784/ 992-1809/1073) -78.7 W47-W50 N12-N17
      5  261 1021 (S24E51 L=138)   804    277        1 ( 183/ 959- 342/1069)  31.6 E58-E44 S28-S20
      6 1444 1021 (N04W24 L=213)   225     28        0 (1428/ 997-1460/1044) -45.0 W22-W24 N02-N05
      7 1465  829 (S06W30 L=220)   689    135        0 (1443/ 759-1483/ 896) -54.5 W27-W33 S09-S02
      8 1196  768 (S17W17 L=206)  1027    224        2 (1071/ 732-1302/ 816)  42.3 W10-W21 S21-S11
      9  909  510 (S40W09 L=198)  1665    306        1 ( 775/ 475-1081/ 548)  15.3 E01-W21 S44-S35
     10 1502  453 (S24W49 L=238)  1733    391        0 (1351/ 423-1642/ 501)  24.4 W38-W60 S29-S18
  

    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