Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1743 (N39E23 L= 96)  2186    443        4 ( 933/1611-1163/1833)  -7.4 E35-E11 N30-N43
      2  520 1660 (N18E54 L= 64)  2238    426        6 ( 434/1596- 577/1748)  85.2 E59-E47 N13-N25
      3  609 1541 (N13E41 L= 78)   374     95        1 ( 572/1511- 653/1566)  59.0 E43-E37 N10-N15
      4  673 1451 (N09E33 L= 85)   668    234        1 ( 648/1364- 704/1547) -45.0 E37-E27 N04-N14
      5  719 1477 (N12E31 L= 87)   438    125        1 ( 693/1421- 754/1532) -81.9 E31-E29 N08-N15
      6  297 1484 (N03E65 L= 53)   850     76        0 ( 271/1446- 324/1524)  90.0 E68-E62 N00-N05
      7 1496 1248 (N20W21 L=139)   634    122        1 (1455/1208-1540/1309) -20.6 W17-W25 N17-N23
      8  252 1173 (S16E56 L= 63)   401     87        1 ( 223/1129- 275/1213)  71.6 E57-E54 S19-S13
      9 1602  446 (S19W58 L=176)   578     98        1 (1558/ 426-1645/ 465)  39.8 W55-W61 S21-S16
     10 1201  251 (S44W45 L=163)  1634    297        4 (1075/ 223-1319/ 287)  22.4 W34-W52 S48-S39
  

    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