Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1361 1738 (N47W44 L=131)  1746    308        6 (1321/1645-1405/1816)  52.4 W38-W52 N39-N52
      2  945 1560 (N41W02 L= 89)   303     73        1 ( 926/1527- 958/1592)  68.2 W00-W04 N38-N43
      3 1140 1495 (N34W15 L=102)  1015    224        2 (1054/1406-1224/1567) -54.0 W10-W19 N26-N39
      4 1707 1259 (N10W51 L=138)  1324    329        2 (1588/1170-1823/1341) -37.9 W42-W60 N02-N16
      5  212 1268 (N28E60 L= 27)   451    112        2 ( 186/1248- 241/1289)  -8.1 E63-E56 N26-N29
      6 1292 1206 (N13W19 L=106)   290    115        1 (1228/1174-1338/1235)   7.1 W14-W22 N11-N14
      7  794 1108 (N14E13 L= 74)   208     87        0 ( 774/1062- 810/1160) -80.5 E14-E12 N10-N16
      8  369 1044 (N14E44 L= 43)   670    166        0 ( 310/1000- 423/1102) -41.2 E49-E39 N10-N17
      9 1570  827 (S14W33 L=120)   283    108        1 (1533/ 805-1614/ 854) -33.7 W30-W36 S15-S11
     10  419  706 (S07E45 L= 42)   335     64        1 ( 398/ 684- 445/ 731)  26.6 E47-E42 S08-S06
     11 1570  581 (S31W33 L=120)  2245    316        4 (1396/ 492-1717/ 710)  19.1 W19-W45 S34-S24
     12  795  491 (S25E23 L= 64)   779    142        1 ( 734/ 457- 865/ 535)  11.3 E27-E17 S27-S22
  

    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