Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1131 1819 (N47E22 L=323)   365     93        0 (1113/1787-1147/1844) -56.3 E24-E19 N43-N49
      2 1298 1762 (N48E06 L=339)   450     80        1 (1275/1736-1323/1789)  90.0 E07-E03 N44-N51
      3  401 1643 (N15E62 L=283)  1221    206        2 ( 340/1606- 445/1697)  32.0 E67-E58 N11-N18
      4  962 1341 (N09E12 L=333)   537    229        1 ( 893/1248-1001/1437) -36.0 E18-E07 N04-N13
      5  303 1249 (S11E53 L=293)   827    174        2 ( 246/1214- 359/1278)  16.7 E57-E47 S13-S08
      6 1373 1167 (N10W16 L=  1)   395    257        2 (1358/1087-1391/1258) -60.9 W13-W18 N06-N15
      7 1210 1115 (N03W08 L=353)   287     85        0 (1177/1065-1244/1155) -38.7 W05-W10 N00-N04
      8  527  880 (S27E27 L=318)   607    291        1 ( 412/ 858- 663/ 901)  12.5 E35-E17 S30-S24
      9 1280  793 (S13W21 L=  6)   827    214        1 (1229/ 706-1356/ 888) -30.3 W15-W27 S15-S08
     10 1453  459 (S25W46 L= 31)   370     76        0 (1415/ 444-1494/ 473)  18.4 W42-W48 S26-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